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A study was made of the Raman spectrum of sulphuric 
acid at concentrations varying every five percent from 100 
to 5 percent by volume. A dispersion of 10A/mm at 4200A 
and 15A/mm at 4500A was used. Eight lines were found, 
which varied gradually in frequency, width, and intensity 
with changing concentration. Their frequency ranges are 


408-432, 561-596, 913-895, 970-973, 980-983, 1031-1051, 
1142-1174, 1216-1236, 1352-1332. A gradual change in 
frequency with concentration was found for a number of 
the lines. The measurements were checked with a micro- 
photometer. 





NUMBER of investigators have studied 

the Raman spectrum of sulphuric acid at 
varying concentrations. This substance is of par- 
ticular interest because of the presence of mole- 
cules and two types of ions. In addition to giving 
data on the molecular structure, Raman spectra 
give information regarding the relative abun- 
dance of the ions. Nisi! found that the lines 
change in intensity and width with increasing 
dilution. Ramakrishna Rao? and Woodward and 
Horner’ reported a frequency shift for some of the 
lines. The latter found that this shift could be 
interpreted on the basis of components of slightly 
different frequency which appeared and disap- 
peared at different concentrations. 

In a preliminary report‘ the writers found evi- 
dence of a gradual shift rather than components. 
It was not possible to decide fully whether com- 
ponents existed without a more exhaustive study. 
Consequently, the Raman spectrum was photo- 
graphed every five percent from 100 to 5 percent 
by volume. The conventional set-up with mer- 


—. 


' Nisi, Jap. J. Phys. 5, 119 (1929). 
* Ramakrishna Rao, Ind. J. Phys. 8, 123 (1933). 
* Woodward and Horner, Proc. Roy. Soc. A144, 129 


(1934), 
* Bell and Jeppesen, J. Chem. Phys. 2, 711 (1934). 


cury excitation was used. A prism spectrograph 
was used with a dispersion of 10 A/mm at 4200A 
and 15 A/mm at 4500A. Lines excited by both 
4047A and 4358A were measured as a check on 
the frequencies. The temperature was kept fairly 
constant during exposures which varied up to 
seventy-two hours. Among the difficulties en- 
countered was the continuous background which 
was especially troublesome for 100 percent acid 
(F.P. 10.45°C). The broad water band at low 
concentrations masked some of the lines excited 
by 4047A. The plates were measured on a com- 
parator and the frequencies calculated from an 
iron comparison spectrum. For the sharp lines 
they are accurate to about 3 cm™'; for the broad 
ones to about 5 cm. These frequencies with the 
widths and intensities are given in Table I. 
Measurements made with KHSQ, and Na.SO, 
solutions are given in Table II. 

It will be seen from Table I that all of the eight 
lines have varying characteristics at different 
concentrations. Line 1 is present everywhere 
except at very low concentrations. It could not 
be measured at 100 percent on account of the 
strong continuous background. The intensity of 
the nearby mercury lines made the measurements 
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Mole 
% 


Vol. 
% 





100 





50 








(f) too faint to be measured on comparator 
(w) masked by water band 

(m) masked by mercury line 

( ) uncertain 


at low concentrations uncertain. This line is 
broad and of constant intensity at all concentra- 
tions. It seems to increase in frequency as is 
shown in Fig. 1. Line 2 is similar to line 1. It has 
the same width and is a little more intense. It 
appears at all concentrations from 100 to 10 
percent and gradually increases in frequency. The 
excitation from 4047 coincides with the 4142 
mercury line at high concentrations, but is easily 
separated at low concentrations. The similarity 
of lines 1 and 2 may mean that they originate 
from the same scatterer. 

Line 3 is the most intense line at 100 percent 
and practically disappears for very dilute solu- 
tions. It broadens noticeably and decreases in 
frequency for low concentrations. Its great in- 
tensity at 100 percent seems to indicate a con- 


The columns refer to the several lines observed. For a given con- 
centration the first number is the Av from 22,938 cm~ (43584), the 
second the Ay from 24,705 cm=! (4047A), the third the width of the 
line in cm~!, and the fourth the intensity. 


nection with the molecule. Line 4 appears first at 
about 55 percent and increases in intensity with 
dilution. The frequency remains constant as far 
as can be determined, and the line is uniformly 
sharp wherever it appears. Because its frequency 
checks with the strongest line observed in sul- 
phate solutions (Table II), it has been ascribed 
to the SO, ion. This would indicate the presence 
of a large number of SO, ions from 35 to 5 per- 
cent. From 95 to 85 percent a weak 971 cm™ 


TABLE II. 








Mole % D 


2.2 21,955 
21,890 


0.4 21,958 
23,724 


983 (4) 
1048 (5) 


980 (5) 
981 (4) 
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Fic. 1. Average frequencies of the five most intense lines 
plotted against the concentration in volume percent. 


line is found. It is entirely absent at low concen- 
trations. It should not be confused with the 
981 cm line. The strong continuous background 
may mask it at 100 percent. Line 5 does not ap- 
pear at 100 percent, but shows up faintly in pure 
acid which has been exposed to the air for only a 


few hours. Since this is in agreement with the 
strongest line found in acid sulphate solutions 
(Table II), it has been ascribed to the HSO, ion. 
The intensity variation shows that this ion is 
present in solutions containing only a trace of 
water and that it is present in great abundance 
from 45 to 5 percent. Its frequency seems to 
decrease slightly and then increase. Sufficient 
data are not available to determine whether this 
frequency shift is also characteristic of the 1048 
acid sulphate line. 

Line 6 is never strong but appears at practically 
all concentrations. A shift of 1146 cm~ has been 
reported® for SO, solutions. Line 7 was not found 
at 100 percent, but at intermediate concentra- 
tions it is slightly more intense than line 6. 
Woodward and Horner found evidence for this 
line in their microphotometer curves. Line 8 ap- 
pears very faint at a number of concentrations. 
It is in approximate agreement with the 1340 
cm~' SO, line. The frequency shift of lines 6, 7 
and 8 is open to question on account of the diffi- 
culty in measuring such weak lines. They could 
be measured best in the region excited by 4047 


* Bhagavantam, Ind. J. Phys. 5, 59 (1930). 
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Fic. 2. Reproduction of microphotometer curves of lines 
excited by 4358A. 











because of the absence of background. The 3450 
cm water band first appears at 80 percent. It 
broadens and increases in intensity with dilution. 
Siva Rao® has found a gradual decrease in fre- 
quency for this band, going from 3468 cm at 
32 percent to 3424 cm for pure water. A careful 
search failed to show any lines corresponding to 
S—H or O—H bonds in sulphuric acid. 

If each line consists of components whose in- 
tensities depend upon the concentration there 
should be a particular concentration favorable to 
the simultaneous appearance of two components. 
Photographs taken every 5 percent have failed to 
reveal this. If they do exist their frequency dif- 
ferences must be much less than that previously 
reported. On the other hand a gradual shift is in- 
dicated for some of the lines. This is shown in 
Fig. 1. In most cases the change is most rapid 
around 85 percent. 

Microphotometer measurements were made 
over the range studied. Three representative 
curves are reproduced in Fig. 2. The symmetrical 
shape of all the lines makes the existence of any 
components doubtful. A study of all the micro- 
photometer data seems to verify the frequency, 
width and intensity variations reported in 
Table I. 

The writers wish to thank the Physics Depart- 
ment of Johns Hopkins University for the use of 
their microphotometer. 


6 Siva Rao, Ind. J. Phys. 9, 195 (1934). 
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Raman spectra of a number of organic compounds containing a neopentyl group have been 
obtained. Among these is tetramethyl methane-d,. The spectrum of this compound shows the 
various types of isotope effects to be expected for a polyatomic molecule. It is also apparent 
from the spectrum that the vibrations of tetramethylmethane-d, are less degenerate than the 
vibrations of tetramethylmethane. The spectrum of tetramethylsilicon has 5 lines in the internal 
vibration spectrum. The explanation for the existence of this extra line is probably the same 


as that given by Fermi in the case of CCl,. 





INTRODUCTION 


HE study of organic chemistry has furnished 
us with an enormous number of different 
species of molecules whose geometrical structure 
has been determined quite satisfactorily on the 
assumption of a static molecule. In dealing with 
the vibrations of the atoms in a molecule the 
chemist’s static model is invaluable. From this 
model we are able to predict the number of 
possible types of vibrations and the degree of 
degeneracy in the case of symmetrical molecules. 
The molecules which have been investigated in 
the present work have been chosen because of 
the various degrees of simplicity introduced into 
their vibration spectrum because of their sym- 
metry properties. 

The title of this paper is in one sense a mis- 
nomer, since no pentatomic molecules were in- 
vestigated. Cross and Van Vleck! have treated 
the ethyl halides as though they were triatomic 
molecules with a fair degree of success. If we 
consider a methyl group for example as a single 
atom, then the molecules in the present investiga- 
tion which will be discussed can be considered 
to be of the pentatomic type. 

The molecules whose Raman spectrum has 
been investigated were tetramethylmethane-d,; 
silicon tetramethyl; 2,2,4-trimethylpentene-3; 
2,2,4-trimethylpentene-4;  2,2,4-trimethylpen- 
tane; 2,2,4,4-tetramethylpentane; 3-hydroxypen- 
tane and 2-chloropentane. The structural formu- 
lae of these compounds and a schematic diagram 


*From a dissertation submitted by D. H. Rank in 
partial fulfilment of the requirements for the Ph.D. 
degree at The Pennsylvania State College. 

1 Paul C. Cross and J. H. Van Vleck, J. Chem. Phys. 1, 
350 (1933). 


of their Raman spectra are shown in Table I. 

The experimental technique employed in the 
present work has been described by one of us’ 
previously. The spectrograph had a dispersion of 
14A/mm at 4500. However, in the case of sili- 
contetramethyl and 2,2,4,4-tetramethylpentane 
the spectrograph was fitted with an F=9 lens 
(dispersion 25A/mm at A4500). The use of the 
faster spectrograph in the above two cases was 
necessary because of the small amounts of these 
compounds available and the comparative feeble- 
ness of the scattered light. 

The 2,2,4-trimethylpentene-3 and 2,2,4-tri- 
methylpentene-4 were samples of the compounds 
separated by Whitmore and Wrenn* and were 
considered to be in a high state of purity. The 
2,2,4-trimethylpentane, 3-hydroxypentane and 
2-chloropentane were furnished us by F. C. 
Whitmore from a stock of very carefully purified 
compounds. In the case of the compounds men- 
tioned above samples of approximately 15 cm’ 
size were used to obtain the Raman spectrum. 


Tetramethylmethane-d 

A 4 cm’ sample of material which consisted of 
approximately a 50-50 mixture of tetramethyl- 
methane and tetramethylmethane-d, was used in 
this case to obtain the Raman spectrum. The 
preparation} of this material has been described 
elsewhere.‘ 


2D. H. Rank, J. Chem. Phys. 1, 572 (1933). 

Whitmore and Wrenn, J. Am. Chem. Soc. 53, 31356 
(1931). 

+A preliminary report on the preparation and the 
spectrum of tetramethylmethane-d, was given by one of us 
(R) in conjunction with Mr. G. H. Fleming before the 
Florida meeting of the American Chemical Society. _ 

4 Whitmore, Fleming et al., J. Am. Chem. Soc. 56, 749 
(1934). 
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2,2,4,4-letramethyl pentane 

A 4 cm’ sample of this material’ which was 
prepared and purified by Dr. F. C. Whitmore and 
Dr. H. A. Southgate was used to obtain the 
Raman spectrum. 


Tetramethylsilicon 

This material was specially prepared for us by 
Dr. F. W. Breuer of the Chemistry Department 
and had the following physical constants: 


Nd = 1.37315 at —5.5°C, 
b.p. 26-26.5°C at 720 mm. 


We have tabulated the frequencies for tetra- 
methylmethane along with that of the mixture 
of tetramethylmethane and_ tetramethylmeth- 
ane-d; for comparison purposes. The frequency 
shifts are given in cm~!. The numbers in paren- 
theses are eye estimates of the relative intensity 
of the lines. 


Tetramethylmethane 


335 (10) 414 (3) 
1455 (15) 2712 (5) 
2892 (5) 2911 (10) 


Tetramethylmethane-d ,+ tetramethylmethane 


283 (3) 338 (6) 418 (4) 723 (20) 736 (20) 
930 (15) 1223 (4) 1251 (15) 1302 (3) 1377 (3) 1395 (3) 
1453 (15) =2179(5) = 2521 (2) 2714 (6) 2779 (3) 2796 (3) 
2870 (10) 2900(5) 2917 (15 V.B) 2959 (15 V.B) 

2916 (3) 2961 (5 V.B) 


819 (6) 


Silicon tetramethyl 


202 (20B) 239 (15B) 
1264 (10B) 1427 (15B) 


598 (20) 
2905 (10B) 


696 (15B) 
2963 (10B) 


863 (10B) 


2,2,4-trimethyl pentene-3 


133 (2B) 202 (2) 300 (2) 
568 (8) 760 (20) 820 (4) 
1027 (4) 1071 (6) 1159 (4) 
1357 (4) 1386 (10) 1395 (6) 
1667 (15) 2903 (10) 2948 (15B) 


327 (8) 
924 (10) 
1204 (6) 
1450 (15) 
3000 (15B) 


506 (2) 
940 (2) 
1257 (4) 
1466 (8) 


2,2,4-trimethyl pentene-4 


124 (4B) 190 (2) 294 (4) 
461 (2) 561 (10) 685 (15) 
894 (4) 910 (8) 936 (8) 
1101 (2) 1155 (8) 1207 (6) 
1297 (2) 1328 (2) 1413 (15) 
2906 (6) 2947 (10B) 2996 (10) 


315 (4) 
765 (20) 
999 (2) 
1240 (10) 
1450 (12B) 
3021 (6) 


387 (4) 
829 (6) 
1049 (2) 
1268 (2) 
1648 (15) 


2,2,4-trimethyl pentane 


302 (2B) 512 (4) 
902 (8) 929 (8) 
1250 (6B) 1358 (4) 
2998 (15B) 


745 (20) 
1101 (6) 
1458 (10B) 


788 (4) 
1167 (4B) 
2905 (10B) 


828 (4) 
1207 (6) 
2943 (15B) 


2,2,4,4-tetramethyl pentane 


274(2) 425 (2) 557 (2B) 


274 (2 736 (20) 
930(5B) 1250 (10B) 1450 (10) 


1480 (5) 


862 (2) 881 (5) 
2915 (10B) 2966 (10B) 


*Whitmore and Southgate (preparation of this com- 
pound to be published). 


PENTATOMIC MOLECULES 


3-Hydroxy pentane 


400 (3) 472 (2) 
965 (4B) 1013 (2) 
1307 (2) 1363 (2B) 

3005 (5) 


2-Chloropentane 


295 (2B) 342 (3) 
442 (3) 511 (2) 
846 (5) 873 (5) 
1087 (3) 1124 (2) 
2914 (5) 2967 (7B) 


829 (2) 
1128 (2) 
2918 (10) 


763 (3B) 
1038 (4) 
1454 (10B) 


857 (4) 
1276 (2) 
2974 (10B) 


399 (3) 
612 (10B) 
1029 (3) 
1309 (3) 


426 (3) 
672 (4B) 
1049 (3) 
1451 (10B) 


380 (3) 
536 (3) 
921 (2B) 
1275 (3) 
3013 (4) 


DISCUSSION OF RESULTS 


The vibrations of pentatomic tetrahedral mole- 
cule are discussed in detail in Born’s Optik.® 
As is well known a symmetrical pentatomic 
molecule has four distinct vibrational frequencies. 
In the case of CCl, for which complete data are 
available the multiplicities and depolarizations 
have been completely determined. The lowest 
frequency has a multiplicity of two and the next 
lowest a multiplicity of three. The totally sym- 
metric vibration is single while the highest fre- 
quency has been shown to be a triple vibration. 
In the infrared spectrum only the two triply de- 
generate frequencies make their appearance 
while in the Raman spectrum all four of these 
fundamental vibrations are observed. The doub- 
let and triplet lines in the Raman spectrum show 
a depolarization of ?. However, the totally sym- 
metric single vibration shows a depolarization 
of zero. 

The Raman spectrum of tetramethylmethane 
shows six lines in the frequency interval 0 to 
1500 cm~'. Of these the four lowest can be as- 
cribed to the four fundamental frequencies of a 
symmetrical pentatomic molecule if we are 
privileged to consider a CH; group as a single 
atom. The infrared spectrum of tetramethyl- 
methane is not available so there is no direct evi- 
dence to allow us to choose the triply degenerate 
frequencies with certainty. However, there seems 
to be no good reason why we should not assign 
the multiplicities of 2, 3, 1, 3 in this molecule for 
the four vibrations in the order of increasing 
frequency as has been shown to be true in the 
case of CCl,. There can hardly be any doubt as to 
the fact that this assignment is correct in the 
case of the higher two of the four frequencies. 

The frequencies 1252 and 1455 cm~ must un- 
doubtedly be ascribed to vibrations of deforma- 


tion of the methy! groups in which the hydrogens 


®’ Max Born, Ein 
Lichtheorie, p. 553. 
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are directly involved. Regarding the multiplicity 
of these vibrations we can say practically nothing 
except that they are probably highly degenerate, 
since in the case of derivatives of tetramethy]l- 
methane we obtain a large number of lines in this 
region which definitely cannot be ascribed to 
splitting up of the fundamental “framework” 
vibrations. 

The most striking feature seen on the plates 
showing the Raman spectrum of these molecules 
containing the “‘neopentyl” group is the high in- 
tensity and the extreme sharpness of the line 
which varies from 721 to 765 cm depending on 
the compound. The 598 cm line of silicontetra- 
methyl shows exactly the same characteristics. 
In the symmetrical compounds no lines of a 
sharpness comparable to this one are found. In 
the case of CCl, Langseth’ has found the 458 
cm line to be exceedingly sharp as is attested by 
his beautiful photographs in which he shows the 


7 A. Langseth, Zeits. f. Physik 72, 350 (1931). 


resolution of the lines which correspond to mole- 
cules containing 1 and 2 atoms, respectively, of 
chlorine 37. It is well known that the usually 
termed name Raman lines is a misnomer and that 
actually a Raman line is an unresolved band 
caused by the overlapping rotational structure. 
It appears from the plates that the lines corre- 
sponding to the totally symmetrical vibration 
consist only of the closely spaced rotational 
branches in contradistinction to the other broader 
lines where apparently the more widely spaced 
branches were also allowed. The existence of this 
single extremely sharp line in the spectrum of 
these compounds shows definitely that rotational 
structure rather than perterbation of the energy 
levels by surrounding molecules in the liquid is 
responsible for the breadth of the other lines. 
The spectrum of silicon tetramethyl shows the 
same degenerate character of the vibration 
spectrum as tetramethylmethane. The internal 
vibrations of this molecule lie at somewhat lower 
frequencies than is the case with tetramethy]- 
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methane. The frequencies 1264 and 1427 cm™ 
which are undoubtedly associated with the 
methyl group vibrations, remain substantially 
the same as in tetramethylmethane. In the case of 
lead tetramethyl apparently, according to Dun- 
can and Murray, ® this is no longer true and these 
vibrations are degraded to lower frequencies as 
well as the internal vibrations. In the case of 
silicontetramethyl we obtain a fairly intensé line 
at 863 cm~' which does not fit into the scheme of 
vibrations for a pentatomic molecule. It seems 
quite probable that the same explanation as that 
used by Fermi® for CCl, would hold in this case. 
If this explanation would be true we must as- 
sume that the highest triply degenerate fre- 
quency of the internal vibration spectrum is still 
unknown but very close to the sum of 202+ 598 
cm! which produces a splitting resulting in the 
production of the two lines, 696 cm™ and 863 
cm~! instead of the single line ordinarily to be 
expected. 


Tetramethylmethane-d , 

The results in the case of this molecule are not 
entirely unambiguous. About 50 percent of the 
molecules were ordinary tetramethylmethane 
and it is not possible to say whether all the lines 
found in tetramethylmethane are absent in the 
deutero derivative or not. The slight differences 
in the frequencies of lines known to be due to 
tetramethylmethane as measured on the two sets 
of plates are in some cases due to experimental 
error but in other cases may be due to superposi- 
tion of the two sets of lines from the two different 
compounds present. For this reason it seems diffi- 
cult to speculate with any degree of certainty 
concerning the splitting up of the internal vibra- 
tions of the molecule due to partial destruction 
of the symmetry by the substitution of a deuter- 
ium atom for a hydrogen atom. However, the 
spectrum clearly shows that the vibrational de- 
generacy of the tetramethylmethane-d; molecule 
is not as complete as that of the tetramethyl- 
methane molecule. The lines 1223, 1302, 1377, 
1395 cm are almost certainly due to the 
splitting up of the highly degenerate frequencies 
which lie at 1252 and 1455 cm~, respectively, in 





*A. B. F. Duncan and J. W. Murray, J. Chem. Phys. 2, 
636 (1934), 


*E. Fermi, Zeits. f. Physik 71, 259 (1931). 





the tetramethylmethane molecule. The totally 
symmetrical vibration exhibits an isotope effect 
of 13 cm™. The line appearing at 2179 cm“ is 
undoubtedly a vibration of the type commonly 
called C—-H ‘“‘stretching vibrations’”’ in which the 
deuterium atom forms part of a simple oscillator 
and is vibrating against part or all of the molecule 
with little or no interaction with other vibrations. 
The line at 2521 cm~ is probably of the same 
type but with the added restriction of interaction 
with other vibrations of the molecule. 

In the case of monochlorotetramethylmeth- 
ane, 2,2,4-trimethylpentene-3, 2,2,4-trimethyl- 
pentene-4, and 2,2,4 trimethylpentane the large 
number of lines shows that the degeneracy of the 
tetramethylmethane molecule has been practi- 
cally completely removed. An inspection of 
Table I seems to show a grouping of lines of 
varying separation which might be traced to the 
fundamental degenerate frequencies of tetra- 
methylmethane. 


CONCLUSION 


This series of compounds has been studied with 
the idea in mind to treat them as pentatomic 
molecules. The spectrum of the symmetrical 
molecules tetramethylmethane and silicon tetra- 
methyl has been satisfactorily explained by con- 
sidering the vibrations of a pentatomic molecule. 
The spectrum of tetramethylmethane-d; shows 
the various types of isotope effects expected for 
a polyatomic molecule. It also shows that the 
substitution of a deuterium atom for a hydrogen 
atom in tetramethylmethane partially destroys 
the dynamic symmetry and as a consequence 
partially removes the vibrational degeneracy. 
The large number of lines appearing in the spec- 
trum of the other tetramethylmethane deriva- 
tives shows practically complete removal of vi- 
brational degeneracy. , 
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Scattering of X-Rays by Polyatomic Liquids. n-Heptane 
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A Fourier integral analysis of Katzoff’s data for the 
x-ray scattering of n-heptane gives not only the arrange- 
ment of molecules in the liquid, but also the most probable 
distances within the single molecule. It is shown that the 
molecule consists of a zigzag carbon chain. Adjacent mole- 
cules are somewhat oriented and increase in temperature 


makes a more random condition. The Fourier integral 
theorem is shown to be valid even though the function for 
which it is used does not converge within the limits of the 
data. The distribution function is used to calculate the 
intensity, giving a fair agreement with the experimental 
values. 





ECENT publications! ? in the field of x-ray 
scattering by liquids have led the writer to 
the questions: 


(a) How much information concerning the structure of 
single molecules may be obtained from liquid scatter- 
ing data? 

(b) Is the use of the Fourier integral theorem valid for 
data extending over a small range of sin 0/d? 


The data of Katzoff* for n-heptane have per- 
mitted calculations that furnish definite answers 
to the questions. These data are particularly 
excellent for they are taken with monochromatic 
x-rays, the absorption in the sample is deter- 
mined, and they extend over a large angular 
range. 


THEORY 


Zernike and Prins‘ and Debye and Menke’® 
have developed equations for the x-ray scattering 
of atomic liquids in terms of a distribution func- 
tion. Warren' and Zachariasen® have applied 
these in slightly modified form to the scattering 
of molecular liquids having atoms of equivalent 
scattering power. The Fourier integral theorem 
is used to solve these equations for the distribu- 
tion function. 

The fundamental scattering equation 


Taw /NTP-1=(1/K) (a(o)/ sin Krdr (1) 


gives on transformation 


1B. E. Warren, Phys. Rev. 44, 969 (1933). 

2 W. H. Zachariasen, Phys. Rev. 47, 277 (1935). 

3S. Katzoff, J. Chem. Phys. 2, 841 (1934). 

a Zernike and J. A. Prins, Zeits. f. Physik 41, 184 
(1927). 

5 P. Debye and H. Menke, Physik. Zeits. 31, 797 (1930). 


a(r)=(2r/x) [ K (Ivovs.) / NT F? —1) 
7 Xsin KrdK, (2) 


where: 

Tops.) is the observed scattering, corrected for absorption, 
polarization and incoherent scattering; 

N is the number of irradiated atoms; 

T is the scattering per electron; 

F is the atomic scattering factor for the given value of K; 

K is 4x sin 0/X, where @ is half the scattering angle and } is 
the wavelength; 

q(r) is the number of atoms in excess of the average lying 
at distance ry Angstroms from a given atom; 

Twons,)/ NTF? is for brevity denoted by £, in keeping with 
common usage. 


Eq. (2) is identical in meaning with Eq. (2) of 
Katzoff’s paper, since the average number of 
atoms at 7 is 4zr’p (p is atoms per A*). 

In practice the data do not extend to infinity, 
nor does K(E—1) converge at the largest ob- 
tainable values of K. The results, however, 
show that no great error is introduced by inte- 
grating from K=0 to K =2r. 


METHOD OF CALCULATION 


The E-curve is determined as shown in Fig. 1. 
Katzoff’s intensity curves for n-heptane were 
replotted to.a convenient scale and the calculated 
F? and incoherent scattering curves drawn to the 
same scale. At the point s(2 sin 6/A) =1, it may 
be assumed that the atoms scatter independently 
and that the observed scattering is thus the sum 
of the independent coherent scattering plus the 
modified scattering. The point s=1 was chosen 
because on the theoretical gas scattering curve it 
lies half-way between a maximum and a mini- 
mum. The procedure is not very exact for large 
s values because of uncertainty both in the F? and 
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Fic. 1. Scattering curves for n-heptane (Katzoff); E curves 
and E curve calculated from g(r) curve. 


the experimental curves, but subsequent calcu- 
lations showed that the results were chiefly due 
to portions of the scattering curve at smaller s 
values where the accuracy is greater. James and 
Brindley’s® F values were used for carbon and 
hydrogen. In the F? curve the values are for C= 
(Z=8) and the hydrogen is neglected. The E 
curves are obtained by dividing the experimental 
curves by the F? curves. Since both are in the 
same units, £ is in absolute units. 

Integration of Eq. (2) was carried out in 0.25A 
steps for values of 7 from 0.25A to 15A. This 
laborious calculation was materially shortened by 
the use of a Coradi rolling sphere harmonic 
analyzer.” The function K(E—1) was plotted 
against K, (Fig. 2) and expressed as a Fourier 
series by tracing with the harmonic analyzer. The 
resulting expression for K(E—1) was multiplied 
by sin Kr and each term integrated from K =0 
to K=2z. The corresponding g(r) curves are 
shown in Fig. 3. 


asa W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1931), 
‘Available through the courtesy of Professor Henry 
ein of the Department of Economics, University of 
icago. 
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n-HEPTANE 


DISCUSSION 


The g(r) curves have unexpectedly sharp peaks 
at 1.4, 2.5, and 3.8A. These peaks are due to the 
outer parts of the scattering curve, where the 
molecules scatter quite independently, and they 
undoubtedly show within experimental limits the 
most probable distances within a single molecule. 
This is verified by the slight change with tem- 
perature and by calculations made with an arbi- 
trarily drawn K(E—1) curve which minimized 
the fluctuations at large K values. This curve 
(Fig. 2A) follows the —50° curve to K=z and 
diverges in the range where the experimental un- 
certainty is greatest. It gives the distribution 
function of Fig. 3A. The peaks appear at the 
same r-values, but are much less pronounced. 
Since the experimental uncertainty is probably 
not as great as shown in the arbitrary curve the 
true distribution function probably lies between 
the two curves. 

The g(r) peaks at larger distances show the 
concentration of scattering matter in adjacent 
molecules. The peak for the lower temperature 
curve is at r=4.5 to 6.5A, with a mean value of 
5.25 to 5.5A. This distance is somewhat greater 
than the mean diameter of the molecule because 
of the distribution of scattering atoms along the 
length of the adjacent molecule. At the higher 
temperature the peak extends outward, giving a 
mean value of about 5.75A. Elevation of the 
temperature has only increased the amplitude of 
the molecular vibration and given a more dis- 
orderly arrangement. 

The most probable distance (a) of scattering 
atoms may be more directly determined, as 
shown by both Zachariasen and Warren, by ap- 
plication of the Ehrenfest-Keesom equation 


d/2 sin 6(max) =d prage = 0.814. 


The value so calculated from the major peak of 
the scattering curve is the same as given by the 
g(r) curve. Application of this equation to the 
second and third peaks gives distances of 2.54 
and 1.43A, respectively. 

The radial distribution obtained does not imply 
the assumption of a structure for the liquid. Al- 
though they may be in random motion and orien- 
tation with respect to a given point, two adjacent 
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Fic. 2. K(E—1) curves for n-heptane 
smoothed data for —50° 


molecules must be in a somewhat parallel ar- 
rangement because of their shape and size. The 
free space in a liquid is probably small in com- 
parison to the volume occupied by the molecules, 


and there is not sufficient space for adjacent. 


molecules to turn endwise with respect to one 
another. Two widely separated molecules may, 
however, have any orientations. 

The packing of the molecules seems to be in a 
hexagonal arrangement for it appears from the 
q(r) curve that each molecule has on the average 
about six neighbors. This is deduced from the 
area under the peak, about four atoms. Since the 
average number of atoms at a distance of 4.5 to 
6.5A is about twenty-five (as computed from 
molecular volume) the total number at this dis- 
tance is about twenty-nine. This would require 
about six neighboring molecules. 

The results of the present calculation may also 
be interpreted by the cybotactic theory, for if the 


° 


; 0 0 calculated from q(r). 


arrangement were crystalline the distribution 
function would be about the same as found. 
Warren and Gingrich* have applied the Fourier 
integral analysis to the scattering from a solid, 
sulfur, and obtain a distribution similar to that 
given by a liquid. 

The question as to the validity of the use of 
the Fourier integral theorem was raised because 
the data do not extend beyond K = 27 while the 
integration must be performed in the limits 
K=0 to K=~. An attempt to find the magni- 
tude of the probable error was made by calculat- 
ing the intensity function from the g(r) curve. 
The portion of this from r=0 to r=4m was ex- 
pressed as a Fourier series by the harmonic 
analyzer and the integration of Eq. (1) was per- 
formed. The calculated K(E—1) curve is shown 
by the circles of Fig. 2 and by the upper E curve 


* B. E. Warren and N. S. Gingrich, Phys. Rev. 46, 368 
(1934). 
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Fic. 3. Distribution function for n-heptane. The g(r) values give the number 
of atoms in excess of the average at radial distance r. 


of Fig. 1. The calculated intensity is somewhat 
too high by a nearly constant factor but. the 
shape of the curve is the same as the original and 
the peaks are at the same positions. The small 
fluctuations at large r values in the g(r) curve may 
be due to the finite integration. (Another possi- 
bility is in the extrapolation of the scattering 
curve to zero scattering angle.) 


The success of the calculation in giving both 
the inner and outer structures of the liquid em- 
phasizes the importance of the Fourier integral 
method of analysis and shows the need for more 
quantitative data extending to large angles. 

The writer desires to express his thanks to 
Professor W. H. Zachariasen for his helpful 
criticisms. 
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of H.O and HS 


W. C. Price,* Physical Laboratory, The Johns Hopkins University 
(Received March 10, 1935) 


The absorption spectra of formaldehyde and of compounds of the type HOH, ROH, ROR and 
HSH, RSH, RSR have been observed in the region 2300-1000A. An analysis of the spectrum of 
formaldehyde yields a value of 10.83+0.01 volts for the first ionization potential of the C=O 
bond and about 164 Cal./mole for the strength of this bond. The spectra of the alkyl deriva- 
tives of H.S are all shown to be shifted to the long wavelength region relative to those of H,O. 
Rough estimates of their ionization potentials are given and the influence of dipole moment on 
the ultraviolet absorption spectra and ionization potentials of polyatomic molecules is discussed. 





HE author has for some time been investi- 
gating the absorption spectra of various 
organic molecules in the far ultraviolet. Up to 
the present, work of this nature has been done 
almost exclusively with fluorite spectrographs by 
using the H»2 continuum. The instruments used 
in this work are two vacuum grazing incidence 
spectrographs which have dispersions of 11- 
9A/mm and 2.7—2.0A/mm, respectively, in the 
region 2000-1000A. The source of the con- 
tinuum was a Lyman discharge tube of the type 
developed by Collins and Price.! The use of this 
continuum in preference to the He continuum 
enabled the absorption to be followed down to 
1000A. Both gratings were of glass and each was 
ruled with 30,000 lines to the inch. The radius 
of curvature of the smaller was 40 cm and that 
of the larger was two meters. The path lengths 
of the absorbing columns in the two spectro- 
graphs were about 35 and 150 cm, respectively. 
The gas to be investigated is led into the body 
of the spectrograph through a fine capillary and 
pumped out through the main evacuating pump. 
The pressure of the gas used in obtaining these 
absorption spectra is usually of the order of a 
tenth of a millimeter. The author has shown that 
in a large number of cases? the bands obtained 
under these conditions consist of various elec- 
tronic series with a superimposed vibrational 
pattern which go to a limit agreeing fairly well 
with the ionization potential of the molecule as 
obtained by electron impact. 


* Commonwealth Fellow. 

1G. B. Collins and W. C. Price, Rev. Sci. Inst. 5, 423 
(1934). 

2 W.C. Price, Phys. Rev. 45, 843 (1934); 46, 529 (1934); 
Phys. Rev. 47, 444, 510 (1935). 
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A photograph of the absorption spectrum of 
formaldehyde taken under high dispersion is 
given in Fig. 1. The gas was obtained by heating 
paraformaldehyde. The absorption is very strong 
and most of the bands appear at a pressure of 
0.1 mm in a path length of 1.5 meters. On 
account of the diffuseness of the bands the 
analysis is rather difficult. It is however plain 
that the bands converge to a limit around 
1150A. 


TABLE I. Comparison between calculated and observed values 
for the frequency differences between successive 
series members in formaldehyde. 





Difference between 
successive mem- 
bers (cm™!) 


n Intensity Wave number Obs. Cale. 
SERIES (i) 

2 10 71600 
7820 7766 

3 10 79420 
3300 3281 

4 8 82720 
1650 1687 

5 8 84370 
980 980 

6 7 85350 
600 619 

7 6 85950 

SERIES (ii) 

2 10D 64270 - 
13400 1066: 

3 10D 77670 
4080 4142 

4 9D 81750 
2060 2028 

5 8 83810 
1140 1142 
6 8 84950 ; 
700 705 
7 7 85650 
470 466 

8 5 86120 
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The strong bands can be classified into two 
electronic series. The frequencies of their mem- 
bers are given in Table I. As most of the members 
are a few Angstroms broad, the measurements 
can only be relied upon to within 50 wave 
numbers. Also, they all contain a considerable 
amount of fine structure due to the low moment 
of inertia of the molecule about the CO axis. 
Even with the high resolution used here,’ little 
could be done with this and the measurements 
given correspond to the center of the band. The 
higher members of (i) and (ii) conform fairly 
well to the following Rydberg series: 


(i) vo" = 87,830—R (n+0.60)?, 


(ii) ve"=87,710—R, (n+0.30)?,. n=2, 3,4 ete. 


An idea of the accuracy with which these 
formulas represent the series may be obtained 
from an inspection of Table I. 

The first member of series (ii) is somewhat off 
the Rydberg formula. This may probably be 
accounted for by the presence of a dipole in the 
CH. group which means that the electron is 
escaping in a field given approximately by 


v=e r1+(p cos 6) r,? 


bp” being the moment of the dipole, ‘‘7.”’ the 
distance of the electron from the center of the 
dipole and @ the angle between 72 and the dipole. 
It is obvious that the second term becomes 
negligible for the higher members of an electronic 
series, and thus they may be expected to conform 
to the Rydberg series. Further it can be shown 
that the values of the dipole moment usually 
associated with a CH, or CH; group can produce 
quite large effects on the ultraviolet absorption 
spectra and ionization potentials of polyatomic 
molecules. Thus the method used in getting 
these electronic series should be to try and fit 
the higher members into a Rydberg series and 
then to fit in the lower members by their in- 
tensities and the fact that they should fall fairly 
clos» to where the band of the Rydberg series 
would fall. Similarities in vibration pattern also 
hel) where these can be observed. 

The extremely diffuse band at 1745A (57,310 
cm ) can possibly be a lower member of Class I 
though it is very far off the series. If we assume 
the predissociation of this band to be due to 
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Fic. 1. The far ultraviolet absorption spectrum 
of formaldehyde. 


interaction with the ground state the process 
being H.CO—CH.+0O we find the heat of 
dissociation of the C=O bond to be about 
164 Cal. mole which is in fairly good agreement 
with the value ~172 Cal. mole obtained from 
thermochemical data. The author has also ob- 
served that all the bands in acetaldehyde and 
acetone which fall in the region below 1745A 
are predissociated, the strength of the pre- 
dissociation depending on how far the bands are 
removed from this wavelength. The bands at 
wavelengths longer than this are sharp. It is 
thus indicated that the C=O bond is roughly 
of the same strength in all three substances. 

The series (i) and (ii) in HeCO have approxi- 
mately the same limit which is 10.83+0.01 volts. 
This compares favorably with the value 11.3 
+0.5 volts determined by T. N. Jewitt® for the 
ionization potential of formaldehyde by the 
methods of electron impact. 

It is fairly certain that the ionization is from 
the C=O bond. If it were from the CH bond 
the ionization potential would be expected to be 
considerably higher and the bands would prob- 


37. N. Jewitt, Phys. Rev. 46, 616 (1934). 
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ably start at about the same wavelength as the 
corresponding ones in methane‘ (i.e., around 
1300A). Investigations of the spectra of acet- 
aldehyde and acetone’ also indicate that in these 
molecules the excitation is from the C=O bond. 

There are only two electronic transitions for 
which the vibrational pattern is at all well 
developed. These are the bands at 1556A and 
1397A. There are two frequency differences 
accompanying each. In the case of the first band 
the differences are 1250 cm™ and 1470 cm™ and 
for the second band they are 1120 cm™ and 
1260 cm~'. These vibration frequencies are pre- 
sumably types which involve vibration in the 
C=O bond. Herzberg*® has found a frequency of 
1187 cm“ in the near ultraviolet absorption of 
formaldehyde which he attributes to a vibration 
in this bond. A few frequency differences of 
~2900 cm involving vibrations in the CH 
bond were observed in the bands reported here 
but they are very weak. 

Although about 150 bands could be observed, 
the value of giving a complete table of their 
wavelengths is questionable because the latter, 
even under high dispersion, cannot be deter- 
mined very accurately as the bands do not 
possess sharp edges. A few things however can 
be said about particular cases: The band at 
64,270 cm™ is really a diffuse doublet with a 
separation of 80 cm~, the frequencies of the com- 
ponents being 64,230 cm™ and 64,310 cm". 
That at 71,600 cm" is likewise composed of two 
bands separated by 50 cm. Also all the upper 
bands of series (ii) have multiple heads which 
are possibly analogous to the 7R, etc., branches 
reported by Dieke and Kistiakowsky’ for the 
near ultraviolet bands. 

The spectra of the alkyl derivatives of H.S 
and H,0O are given in Fig. 2. They were taken on 
the smaller spectrograph. It is immediately 
evident that the spectra of the compounds con- 
taining sulphur are all shifted to the long wave- 
length region relative to those containing oxygen. 
This is in accord with the fact that the ionization 


4A. B. F. Duncan, J. Chem. Phys. 2, 851 (1934). 

5G. Scheibe, F. Povenz and C. F. Linstrom, Zeits. f. 
physik. Chemie B20, 292 (1933); W. A. Noyes, A. B. F. 
Duncan and W. M. Manning, J. Chem. Phys. 2, 717 (1934). 

6 G. Herzberg, Trans. Faraday Soc. 27, 378 (1931). 

7G. H. Dieke and G. B. Kistiakowsky, Phys. Rev. 45, 4 
(1934). 
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Fic. 2. The far ultraviolet absorption spectra of the alky] 
derivatives of H»S and H.O. 


potential of S(1) is ~10.3 volts and that of 
O(1) is ~ 13.6 volts. The spectra of H2S has been 
photographed under high dispersion. A complete 
discussion of it is reserved for publication else- 
where. The bands can be arranged into electronic 
series which go to several ionization potentials in 
the neighborhood of 10.5 volts. This agrees very 
well with the value 10.4 volts* obtained by 
electron impact. 

The bands in H2O have been followed to a 
limit around 980A by Rathenau.” Though no 
analysis has been made of the bands this limit 
is about 12.6 volts and it is obvious that they 
correspond to electronic series going to the 
ionization potential at 12.7 volts determined by 
Smyth and Mueller.'® Thus it is fairly sure that 
the bands of the other molecules are of the same 
type, ie., electronic series going to the first 
ionization potentials. 

There is a very considerable amount of evi- 
dence available for spectra of the type reported 
here, which indicates that the bands of any 
particular alkyl compound start around the 
same wavelength and go to about the same 
ionization potential whatever the nature of the 
alkyl group, e.g., the spectrum of methyl not 
ethyl alcohol is given in Fig. 2. However, the 
spectrum of ethyl alcohol starts at almost the 
same wavelength. Further, the author has photo- 
graphed the absorption spectra of ethyl, propyl, 


5C. A. MacKay, Phys. Rev. 24, 319 (1924). 

*G. Rathenau, Zeits. f. Physik 87, 32 (1933). 

” H. D. Smyth and D. W. Mueller, Phys. Rev. 43, 116 
(1933). 
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isopropyl and butyl mercaptans and sulphides 
and found that they all begin around 2100A and 
2300A, respectively. Also Scheibe, Herzberg" 
and others have obtained the spectra of methyl, 
ethyl and tertiary butyl iodide and found that 
they all start around 2010A. This must be 
related to the fact that the excitation is from 
the non-alkyl part of the molecule, and also that 
the dipole moments of the various alkyl groups 
are not sufficiently different to affect the term 
value of the excited electron.* 

While it is still too early in the development of 
the spectra of polyatomic molecules to make 
many definite statements, some tentative re- 
marks on these spectra may however be in order. 
The excitation in the case of the alcohols and 
the mercaptans is surely from the CO and the 
CS bonds, respectively, as the other bonds have 
higher ionization potentials. That of the sul- 
phides and the ethers is also from the CO and 
CS bonds but the substitution of an alkyl group 
for a hydrogen atom has introduced a dipole 


1G. Herzberg and G. Scheibe, Zeits. f. physik. Chemie 
B7, 390 (1930). 

* A dipole of 0.3 10~"* e.s.u., such as is usually associ- 
ated with a methyl group, will change the term value of an 
electron situated 1 Angstrom away by almost one electron 
volt and can thus change the position of the absorption 
bands by more than 100A. 
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which reduces the term values in the adjacent 
bonds (the negative end of the dipole is towards 
the excited electron). This diminishes the ioniza- 
tion potential and shifts the absorption bands 
towards the red. This successive shift towards 
the red is also found in the spectra of HeCO, 
CH;3HCO and (CHs3)2CO. Their estimated ioniza- 
tion potentials are 10.83, 10.5 and 10.1 volts, 
respectively. 

It is possible to predict very roughly that for 
alcohols and ethers the first ionization potentials 
are around 10.8 and 10 volts, respectively, and 
for mercaptans and sulphides around 9 and 8.5 
volts, respectively. 

Preliminary measurements with fluorite spec- 
trographs have been made on methyl and ethyl 
alcohol by Herzberg and Scheibe" and on methyl 
and ethyl ether by Scheibe and Grieneisen.”” 
A wavelength table is not given here as it is 
intended in the near future to photograph the 
bands with higher dispersion and at a greater 
variety of pressures. 

In conclusion the author wishes to express his 
gratitude to Dr. G. H. Dieke and Dr. K. F. 
Herzfeld for many discussions on the material of 
this article. 


12 Scheibe and Grieneisen, Zeits. f. physik. Chemie B25, 
52 (1934). 
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The Calculation of the Free Energy of Polyatomic Molecules from Spectroscopic 
Data. II 


A. R. Gorpvon, Chemistry Department, University of Toronto 
(Received January 30, 1935) 


The methods already developed for computing the free energy from spectroscopic data for 
diatomic molecules and for triatomic unsymmetrical rotators are extended to handle the case 


where one or more of the fundamental frequencies is degenerate. New tables, necessary for such 
calculations, have been computed and the procedure is illustrated by sample calculations for 


nitrous oxide, hydrogen cyanide and acetylene. 


URING recent years, the calculation of 
thermodynamic quantities from spectro- 
scopic data has proved to be a most useful method 
of solving problems of chemical equilibrium. At 
the moment of writing, with the exception of 


steam,' carbon dioxide* and nitrous oxide,’ 
precise calculations, taking account of the second 
and third order terms in the expression for the 


1 Gordon, J. Chem. Phys. 2, 65 (1934). 
2 Kassel, J. Am. Chem. Soc. 56, 1838 (1934). 
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energy levels, have been carried out only for 
diatomic molecules; consequently the extension 
of approximation methods to handle the more 
complicated polyatomic problem is clearly a 
matter of some importance. Steam (a typical 
example of a triatomic unsymmetrical rotator) 
is in some respects a specially simple case, since 
the question of vibrational degeneracy does not 
arise, and as a result the calculation can be 
brought under the diatomic form. For almost all 
other types of polyatomic molecules, however, 
one or more of the fundamental frequencies will 
be degenerate, and, as will be shown below, the 
rough assumption that the contribution of a 
doubly degenerate mode of vibration will be twice 
that of a singly degenerate mode of vibration 
with the same spectroscopic constants, may lead 
to appreciable errors in the calculated thermo- 
dynamic quantity. In the present paper, the 
methods previously developed for diatomic 
molecules‘ and for the triatomic unsymmetrical 
rotator! are extended to deal with the case where 
one (or more) of the frequencies is doubly de- 
generate; this extension has required the calcula- 
tion of new tables which are appended. 

The examples given below to illustrate the 
procedure are all for the collinear type of mole- 
cule, since it is only for these that a reasonably 
complete analysis of the energy levels has been 
made; the quantities tabulated in this paper, 
however, may be used in calculations for any 
type of molecule possessing a doubly degenerate 
mode of vibration. In the examples below, in 
order to illustrate the precision of the approxima- 
tions, the calculation is sometimes carried to a 
greater number of significant figures than is 
justified by the accuracy of the experimental 
data. 


THE TRIATOMIC COLLINEAR TYPE, 
HCN anp N,O 


The evaluation of the State Sum for this kind 
of molecule has recently been considered by 
Kassel,? who has computed the free energy, etc., 
of carbon dioxide and nitrous oxide by an exten- 


’ For an excellent general survey of the present state of 
spectroscopic thermodynamics, see a recent series of 
review articles by Zeise, Zeits. f. Elektrochemie 39, 758 
(1933); 39, 895 (1933); 40, 662 (1934); 40, 885 (1934). 

* Gordon and Barnes, J. Chem. Phys. 1, 297 (1933). 
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sion of the methods he has previously used. Car- 
bon dioxide is very much a special case owing to 
the resonance between two of the fundamentals, 
and since Kassel has dealt with it in detail, it 
will not be considered here. In general the vibra- 
tional levels for such molecules are given by an 
expression of the type 


ev /he = 0101 +017X11+V2W2+V2"Xo0+/?x, 
+U3w3t+3"X33+V2X12 
Hv 3Xis +2303. (1) 


Here 7; and v; are the quantum numbers for the 
two nondegenerate “‘valence’’ frequencies v; and 
v3, while ve is for the doubly degenerate ‘‘defor- 
mation” frequency »,. Such a degenerate fre- 
quency requires a subsidiary quantum number / 
to define it, / taking (for a given v) the v2+1 
values —vo, —ve+2, «++, Y2—2, vo; V1, Ve and 2; 
take all values from zero to their respective series 
limits. 

The rotational levels are in general given by an 
expression of the ordinary diatomic form® 


«,/hc=B,J(J+1) +DI2°(T+1)? (2) 


where B,=Bo(1—aywvi—aede—azv3+---) and 
J=|1l|. When no question of symmetry nor of 
nuclear spin arises, the summation over J for 
moderate and high temperatures may be re- 
placed without sensible error by an integration: 


> (2J+1) exp [— {B,J (J+1) 


J=1u 
+DJ*>(J+1)*}he/kT } 
we?) aye el 4ay-g,(14+1/12¢,+d,), (3) 


5 Neglecting the D term, the rotational energy is 
(strictly speaking) Bhc(J?+J—F); see Dennison, Rev. 
Mod. Phys. 3, 280 (1931). For the purposes of computation, 
however, the term scheme can be represented most 
conveniently by Eqs. (1) and (2). See, for example, Choi 
and Barker, reference 12, pp. 780-783 and Adel and Barker, 
reference 11 p. 278. In order to represent the rotational 
energy accurately, it may be necessary to have terms 
involving powers of J other than those of Eq. (2); such 
terms will introduce in Eq. (3) items similar to those 
discussed in Gordon and Barnes’ example 6, reference 4. 

6 The approximation of Eq. (3) decreases in percentage 
accuracy with increasing |/| and decreasing g. For ex- 
ample, for g=100 and D =0, the values of (3) are 


|| 0 1 2 3 4 
Direct summation 100.334 99.334 96.393 91.685 85.476 
Eq. (3) (approx.) 100.334 99.336 96.400 91.698 85.499 


For g=500 on the other hand, and |/| =4, the exact value 
is 484.52 while approximation (3) gives 484.58. : 
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where q,=kT/B,hc, dy= —(2D/By)qv. As usual, 
for each atom possessing j, units of nuclear spin, 
Eq. (3) must be multiplied through by a factor 
(2j,+1) if “absolute” thermodynamic quantities 
are desired. 

In the rotational-vibrational state sum, the 
/-summation occurs as the factor 


x exp [—/?(y:+1/qp) ], (4) 


where y,=x,hc/kT; this can be expanded and the 
resulting series in /*, /*, etc., can be summed as 
shown by Kassel.? The exponentials of the cross 
terms involving x12, etc., can be expanded simi- 
larly, and by an obvious rearrangement the 
logarithm of the state sum can be written to a 
sufficiently good approximation 


3 
In Qrv+Eo/RT = Yi (In Qvy + €0/RT); 
jon 


+1n got+1/4qot+ln (14+f), (5) 
where 
tet!*?. Ov} ;= {Do p.-exp [(eo— ev) /RT J}; 
= {>'p,-exp [—vw+(v?+v)wx]};, (6) 
f=1/12go+do—(yr+1/qo) (v2?+ 202) /3 
+ (yi +1/qo)?(302'+ 1202+ 8022 — 802) /30 
+ 120102 +130 103+ Ye30003 
+aibit+ari2+asb3, (7) 
vi" =(Lpsvr-e vl} /{ Spee vit}, (8) 
w=hcw/kT, (9) 
for i=1, 3 | 
=v+1 
Vij = —xijhc/kT. 


for 1=2, (10) 


(11) 


The replacement in f of g, and d, by gq and dy 
(the values of these quantities in the ground 
vibrational state) introduces only negligible er- 
rors in the final result. The neglect of some of the 
higher order term not included in Eq. (7) may be 
a more serious matter; such items, however, can 
be written out in any actual case where there is 
any doubt, e.g., tefms such as — yi2(y,;+1/qo) 
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X (6:29 + 25,022) /3, —-yr2"710?/2, —-Yr2V1901 Dads, 
123010203, a1V127;2b. to mention only a few, and 
their values estimated. It will usually be found 
that such terms are relatively unimportant, and 
also that there is some mutual cancelation. 

Of the quantities on the right of (5), (In Qy 
+¢€/kT), 6 and v? can be found at once in the 
case of »; and v3 from Gordon and Barnes’ 
tables for w=0.6 and x=0.02. For , the same 
three quantities for the same ranges of w and x 
may be obtained from Tables I, II and III of this 
paper. It should be noted that in entering the 
tables, the w defining w is not the w, of Eq. (1) © 
but is the difference (w;—x,,) of that equation, 
and that the x of the tables is the ratio —x;;/ 
(w.—xii) of the equation (cf. Eq. (6) above). 
The only remaining untabulated quantities in 
Eq. (5) are v° and 24 in the fourth term of 
Eq. (7). The contribution of this term is in gen- 
eral small and in most cases negligible; hence no 
serious error is introduced in estimating it by 
assuming that the vibrational levels of v2 are 
those of a harmonic oscillator with frequency 
(we+%22). On this assumption it is possible to 
write down closed expressions for 2‘, etc., and 
the term can be replaced by 


(yr+1/qo)?(t+130?+ 2403+ 1244), 


where t=e~—”/(1—e7*). 


(12) 


Tables I, II and III were computed by summa- 
tion (up to the series limit where necessary) for 
the first three entries in each column and then 
for every alternate entry in the column; the re- 
maining numbers were obtained by interpolation, 
using the method of mean second differences. 
All entries were originally computed to an extra 
decimal place to ensure accuracy in the last place 
as printed. The intervals are sufficiently fine to 
permit interpolation using not more than second 
differences. For values of w and x outside the 
range of the tables, the integral approximations’ 
discussed by Gordon and Barnes may be used in 
the case of v; and v3; for ve, similar integral ap- 
proximations can readily be set up. For example, 
the integral required to evaluate Qy is in their 
notation J,+J,; similarly, in the evaluation of 


7 There is an obvious misprint in their expression for J2 
at the bottom of p. 304. The equation should read 


T2=(s?+ 3) J, —(1/2r)-{(s—Mr)-EEm]—(s—Nr)-(En]}. 
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TABLE I. In Qy+e0/RkT as a function of wand x when p,=v+1. 








0.0025 0.0050 0.0075 


0.0100 0.0125 0.0150 0.0175 0.0200 





1.6130 
1.5645 
1.5182 
1.4738 
1.4312 


1.6358 
1.5864 
1.5392 
1.4939 
1.4505 


1.6604 
1.6099 
1.5616 
1.5155 
1.4712 


1.3903 
1.3510 
1.3132 
1.2768 
1.2418 


1.4089 
1.3689 
1.3305 
1.2936 
1.2580 


1.4287 
1.3880 
1.3489 
1.3113 
1.2751 


1.2081 
1.1756 
1.1442 
1.1139 
1.0847 


1.2237 
1.1907 
1.1589 
1.1282 
1.0985 


1.2402 
1.2066 
1.1743 
1.1431 
1.1129 


1.0564 
1.0291 
1.0027 
0.9772 

-9524 
0.9285 


1.0698 
1.0421 
1.0153 
0.9894 

-9643 
0.9400 


1.0838 
1.0557 
1.0285 
1.0022 
0.9767 
0.9521 


1.6875 
1.6357 
1.5862 
1.5389 
1.4936 


1.7182 
1.6647 
1.6136 
1.5648 
1.5182 


1.7538 
1.6980 
1.6449 
1.5943 
1.5461 


1.7953 
1.7368 
1.6813 
1.6285 
1.5782 


1.5303 
1.4845 
1.4408 
1.3989 
1.3588 


1.4502 
1.4086 
1.3686 
1.3302 
1.2933 


1.4737 
1.4310 
1.3900 
1.3507 
1.3130 


1.5000 
1.4560 
1.4138 
1.3734 
1.3346 


1.2578 
1.2236 
1.1906 
1.1589 
1.1282 


1.2767 
1.2418 
1.2081 
1.1757 
1.1444 


1.2973 
1.2615 
1.2271 
1.1939 
1.1619 


1.3204 
1.2835 
1.2480 
1.2139 
1.1811 


1.0986 
1.0700 
1.0424 
1.0156 
0.9898 
0.9648 


1.1142 
1.0851 
1.0570 
1.0298 
1.0035 
0.9781 


1.1311 
1.1014 
1.0727 
1.0449 
1.0181 
0.9922 


1.1495 
1.1190 
1.0896 
1.0612 
1.0338 
1.0073 


1.1696 
1.1383 
1.1081 
1.0790 
1.0509 
1.0238 





x= 0.000 0.005 0.010 0.015 

0.9310 
-8746 
-8224 
-7739 
-7288 


0.020 


0.8614 
-8095 
-7615 
-7168 
-6752 


0.8826 
-8294 
-7801 
.7344 
-6918 


0.9056 
-8510 
-8003 
-7533 


0.9598 
9011 
-8468 
-7966 
-7500 


-6364 
-6002 
-5663 
-5346 
-5050 


-6521 
-6150 
.5804 
-5480 
-5176 


-7067 
-6664 
-6288 
-5937 


-6869 
-6479 
6114 
-5773 
-5454 


4771 
-4510 
-4265 
-4035 
-3818 


-4892 
-4625 
4374 
-4139 
3917 


-5156 
-4875 
-4612 
-4365 


SUID 


3614 
3421 3512 
-3240 3327 
-3069 «3152 
0.2908 0.2987 


-3708 


5 
0 
5 
0 
5 
0 
5 


Seses 





x= 0.000 0.005 0.010 0.015 0.020 


~ 


NNNNhY 
keane 


0.2613 
-2349 
-2113 
-1902 
-1713 


0.2685 
-2415 
.2174 
-1958 
-1765 


0.2761 
-2485 
-2238 
-2017 
-1819 


0.2841 
-2558 
-2305 
-2079 
-1876 


0.2925 


-1543 -1591 
-1391 -1435 
-1255 1295 
-1132 -1169 
-1021 -1056 


-1693 
-1529 
-1382 
-1250 
-1130 


WQNNNN 
Senna 


-1022 1059 
0925 0959 
-0837 .O869 
0758 .O787 
-0687 0713 


.0922 
-0832 

0752 
-0679 
0613 


-0954 
-0862 

0779 
-0704 
-0637 


OwWwaw 
ad 


-0661 


-0576 
-0521 
-0471 
.0426 
0.0385 


-0598 
-0541 
-0490 
-0444 
0.0402 


.0622 .0646 
0563 0586 
.0510 0531 
-0462 0482 
0.0419 0.0437 


-0554 
0501 
.0453 
-0409 
0.0370 


wm Gd Go Go Go 
Soe ND 








> (v+1)-v"- exp [—vw+ (0? +0) wx], 


where n=1 or 2, required in computing % and 
v*, the integrals involved are their Z,+J; and 
(2s+r)I;/r—TI;/r?, respectively. 

The calculation of the entropy, heat content 
and heat capacity for molecules with a degenerate 
fundamental frequency would require additional 
tables similar to those constructed by Gordon 
and Barnes for the diatomic case, and the calcu- 
lation of such tables is already under way. It 
should be noted in passing, however, that a table 
of F°/T as a function of temperature, accurate 
to 0.0005 cal./deg. will yield on tabular differen- 
tiation® values of the derivative (and conse- 
quently of S° and H/T) numerically accurate to 


* See any standard text on tabular methods, e.g., Rice, 
Theory and Practice of Interpolation, pp. 107-113 (1899). 


0.005 cal./deg., i.e., a value of the heat content 
even at 2000° reliable to 10 cal.; for most prac- 
tical purposes, such accuracy is all that is 
required. 

It is important to note that the principal dif- 
ference between the method outlined here and 
that used by Kassel lies in the fact that he takes as 
a first approximation the vibrational energy 
levels to be those of a series of uncoupled har- 
monic oscillators, and in evaluating the state 
sum expands the exponentials of all anharmonic 
and cross terms as a correcting factor. In the 
present paper the first approximation consists in 
assuming that the vibrational levels are those of 
uncoupled anharmonic oscillators, and the ex- 
ponential expansion applies only to the cross 
terms and the term involving x,, The great ad- 
vantage of Kassel’s method lies in the fact that 
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TABLE II. 0 as a function of w and x when p,=v+1. 








x =0.0000 0.0025 0.0050 0.0075 0.0100 0.0125 0.0150 0.0175 0.0200 x =0.000 0.005 0.010 O15 0.020 
0.60 43 2.513 2.604 711 3.026 3.265 3.570 3.848 ‘ 0.998 1.042 1.093 85: 1.233 
62 od 2.403 2.487 586 3.082 3.358 3.622 ‘ -862 0.898 0.939 .987 1.047 
64 my 2.300 2.379 470 2.915 3.165 3.414 i .749 779 813 8! 0.899 
66 . 2.205 2.278 362 2.763 2.989 3.222 ‘ 655 680 .708 ° .778 
2.115 2.183 261 2.624 2.828 3.045 oa .574 -596 620 647 -678 
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1% 
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wn 


NRNwNN Nh 
NNwhK 


2.030 2.094 .167 2 ax} 497 2.680 2.882 
2.010 O78 85 obs od 2.545 2.731 
1.932 995 , a8 2. 2.420 2.592 
1.858 917 -985 -067 : 2.464 
1.788 843 .907 98. J 2.345 


.506 .525 545 aa 594 
463 481 5 .522 
410 426 44. 461 
.364 .378 39. 409 
.324 .336 ol .363 


Ne ee ee ee 
Ccaunn 
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722 774 83: 1.903 ‘ 2. 2.235 
-600 -646 ‘ 1.759 8: 925 2.038 


1 258 .267 ; 288 
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1.490 531 e 1.631 
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1 


.207 .214 : .231 
.166 -172 ‘ .186 
.134 -139 14 .150 
.109 -113 r .122 


391 428 ‘ 1.516 
301 334 37 1.413 
1.218 .249 : 1.320 


= Nw 
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WNNNLK 
SCRA SN 
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.072 O75 .07 081 
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048 050 52 054 
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TABLE III. v as a function of w and x when p,=v+1. 








x= 0.0000 0.0025 0.0050 = 0.0075 0.0100 = 0.0125 0.0150 0.0175 0.0200 x =0.000 0.005 0.010 0.015 0.020 


Ww 





0.60 11.31 12.14 13.14 A 16.31 19.33 23.52 28.92 32.86 , 2.49 : 2.99 3.i 3.96 
62 10.46 11.20 12. a 14.83 17.37 20.98 25.76 29.54 , 1.98 . 2.34 ‘ 2.98 
64 9.70 10.36 11. : 13.54 15.69 18.79 23.01 26.58 ‘ 1.59 ° 1.86 .0: 2.30 
.66 9.01 9.60 od - 12.40 14.23 16.89 20.61 23.94 1.30 a a d 1.82 
68 8.38 8.92 es ¥ 11.40 12.95 15.26 18.51 21.58 os 1.07 . By 5 mF 1.46 


8.30 8. os 11.84 13.84 16.68 19.49 . 0.89 95 02 ‘ 1.20 
7.74 . , 10.88 12.59 15.07 17.63 . 75 85 0.99 
7.23 7 2 < 10.01 11.49 3. 15.99 d 63 67 ° . 82 
6.76 . ° “ 9.24 10.53 ‘ 14.53 : .54 oS ‘ ‘ 69 

P , 8.56 9.68 32 13.22 . -46 é 52 os 59 


7.94 5 ee 12.06 . 34 ‘ 38 40 43 
6.88 d ° ; ’ -26 : 

6.01 d ¢ De d -20 

5.29 5.7 43 iP ‘ 15 

4.69 5. Des ¥ 3. 0.12 





7.25 
6.33 
5.57 
4.93 
4.38 
3.91 


he is able to write closed expressions for all might be necessary to apply correcting terms to 

items and that as a result the use of tables is his expressions. 

avoided ; the disadvantage lies in the very com- Example 1: N;O at 1000°K 

plicated form of his correcting factor, the analog : x od ; 

of the (1+) of Eq. (5) above (see, for example, The vibrational energy” - rie by co he 

his last equation on p. 1838 -of reference 2). 7 rai — 3.301" + 593.00 — raed + oni 
There is one other difference of minor impor- __ 13.803° + 1.03/° — 9.80102 — 26.701: os papewing 

tance between the two methods. Kassel appar- For he and v3, w of the tables nig 1.8504, 

0.85387 and 3.2254, respectively, and x is 0.00255, 


ently takes all summations over a given vibra- 4 ancon 3 0.00613; hence from Table 1 of this 


tional number to run from zero to : , 

same ere 7 dB ‘eables and inthe POPet and Gordon and Barnes’ Tables I and II, 
; : ; = 
Se ee See | Me a sNEY Eat ny, 1 tel my, to CLT, 1.1901 


tables of this paper, the summations run only and 0.0423, respectively. Similarly, from Table 


from zero to the series limit. While it is very II of this paper and Gordon and Barnes’ Table V 
doubtful that-in any actual case the use of sum- = ———____ 


mations to infinity will introduce any appreciable * Plyler and Barker, Phys. Rev. 38, 1827 (1931) as 
. y y : PP : quoted by Kassel, reference 2. For revised values of the 
error in the calculated thermodynamic quantity, constants, see Barker, Phys. Rev. 41, 369 (1932). 


ca ee ee e a - : , 10 In this and the other examples, the universal constants 
It is r n s mples, “ 
is conceivable that with a combination of a used were h=6.554X10-7, c=2.9986X10", &=1.372 


large anharmonic term and high temperature it x 10-", 
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d is 0.201, 1.563 and 0.043, while from Table III 
of this paper, v2? is 8.35. The rotational-vibra- 
tional interaction terms and the second order ro- 
tational terms involving D, etc., are not known, 
but the moment of inertia is 66.1 10-”. Hence 
the approximate value of (In Qry+Eo/RT) from 
Eq. (5) is 8.7733, while the exact value found by 
summation over 7,023 and / is 8.7731. In the cal- 
culation of f, the cross terms X12, x13 and X23 
contribute +0.00441, +0.00033 and +0.00128, 
respectively ; the term — (y:+1/qo)(v2?+2 2) /3 is 
— (0.00579 and the fourth term in (7) is +0.00010, 
so that f.is only +0.00038. The assumption that 
the levels could be treated as those of three 
harmonic oscillators (two singly degenerate with 
v=1285 and 2224 and one doubly degenerate 
with v=590) plus a rigid rotator with /=66.1 
<10-* would have given the value (In Qrv 
+ E,)/RT)=8.757; the error here (0.016) comes 
almost entirely from the neglect of the an- 
harmonic term in 1. 


Example 2: HCN at 1500°K 


The vibrational levels" are given by 2037.0 2, 
+52.0 v,°+712.1 V2—2.9 Vo? + 3364.2 V3—51.3 V3" 
+3.1 1?—0.2 vye—15.1 v1v3 —19.6 vQv3. The rota- 
tional constants!” are B,= 1.47890(1 —0.00625 a 
— 0.00726 v;), D= —3.63 X 10-*; hence go= 708.08, 
dy= 0.00348. The positive sign for the anhar- 
’ monic term of 7; is a little surprising and leads to 
untabulated quantities; (In Qy+¢/k7) and 6 for 
this fundamental, however, can be evaluated by 
the Kassel method" and are 0.1440 and 0.167, 
respectively. For v, and v3, w is 0.68279 and 
3.2616, and x is 0.00406 and 0.01502, respec- 
tively; hence, from the tables, for v, (In Qy 
4+ €/kT)= 1.4372, 5=2.144 and v?=9.22. For v3 
from Gordon and Barnes’ tables (In Qy+e/kT) 
= 0.0434 and 3=0.045. Hence the approximate 
value of (In Qrv+Eo/RT) is 8.1758, as compared 
with the exact value obtained by summation, 
8.1757. Here the terms x12, X13 and x23 contribute 
+0.00007, +0.00011 and +0.00181 to f; 
—(yr+1/o)(02?+2i2)/3 is —0.01967, 1/12g0 
= 0.00012, a1, 1 is +0.00104, a303 is +0.00033 and 
the fourth term in (7) is +0.00105; thus f is 


1 Adel and Barker, Phys. Rev. 45, 277 (1934). 

12 Badger and Binder, Phys. Rev. 37, 800 (1931); Choi 
and Barker, Phys. Rev. 42, 777 (1932); Herzberg and 
Spinks, Proc. Roy. Soc. A147, 434 (1934). 

13 Kassel,” Phys. Rev. 43, 364 (1933). 
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— 0.01166. Here the harmonic oscillator-rigid ro- 
tator approximation would give (In Qrv+Ey/RT) 
= 8.170, with a resulting error of only 0.01 cal./ 
deg. in F°/7; the reason for such close agreement 
in this case is that there is a fortuitous cancela- 
tion of the second and third order terms. 


THE TETRA-ATOMIC COLLINEAR TYPE 


Example 3: C,H; at 1000°K 


In spite of the very considerable amount of 
work that has been done on the acetylene spec- 
trum, no entirely satisfactory analysis of the 
vibrational levels has been accomplished as yet. 
For the purpose of a numerical example il- 
lustrating the use of the method, however, either 
the Sutherland'* or the Herzberg-Spinks" as- 
signments may be used. Sutherland gives 


ey /hc = 19740, + 342102 —49.2v0? + 3293v3 
— 5v3?+ 73Cv4+ 6055 — 30172 — 90103 — 20105 
— 146.8v2v3 — 10v0v4+ 2720; — 1337; — 9 47;,. 





Here »;, ve and v3 are the three nondegenerate 
valence frequencies and » and »; are the two 
doubly degenerate deformation frequencies. 
These latter are characterized by the subsidiary 
quantum numbers /, and /; which play the same 
role’® as the / of Examples 1 and 2. Since the co- 
efficients of /,? and /;? in the expression for ey are 
not known, the analog of Eq. (5) above becomes 


In Qev+E,/RT= E(In Ov+eo/kT): 
+1/4qo-+1n goin (1+f)—In 2, (13) 
where 
f=yidwot +++ ty Dstyi2%1? 172/24- + 
+-yi2V13012 Bads+ + ++ — (V2 +24)/3q0 
— (052+ 285)/3go+ait:+++asi5 
+1/12qo+do, 


(In Qy+¢0/kT), 5, v?, go, yi; and dy are as defined 
in Eqs. (6)—(11), p, is 1 for i=1, 2 and 3, and is 
v+1 for i=4 or 5, and —In 2 is the symmetry cor- 
rection. From Herzberg and Spinks’ rotational 


14 Sutherland, Phys. Rev. 43, 883 (1933). 


16 Herzberg and Spinks, Zeits. f. Physik 91, 386 (1934). 
16 Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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analysis B,=1.17692(1-0.005240, —0.00550v. 
—0.00521v3 —0.00051v,+0.00263v;), D=—1.83 
x10; hence go= 593.18, and d>=0.00184. The 
arguments for entering the tables and the corre- 
sponding tabulated quantities are 


V1 v2 v3 v4 V5 


w 2.8276 
x 0 
In Ov te/kT 


4.9711 

0.0142 

0.0610 0.0081 0.0090 0.8661 
0.063 0.008 0.009 1.084 


v? oe oe am 2.85 

Hence, from Eq. (15), (In Qrv+E)/RT) 
=7.7416, as compared with the exact value 
found by direct summation, 7.7412. The cross 
terms X12 to x3; make a relatively slight contribu- 
tion (+0.00056) to f, but ysis; and (yas)? 
-04°052/2 are +0.02027 and +0.00109, respec- 
tively; the terms involving @ altogether amount 
to —0.00288, while —(v4?+2%4)/3q. and —(v;? 
+245)/3q are —0.00282 and —0.00421. A 
calculation assuming a rigid rotator (B=1.177) 
plus five harmonic oscillators (v=1974, 3372, 
3288, 730 and 605, the last two degenerate) gives 
(In Orv+E,o/RT)=7.728 with a resulting error 
in F°/T of 0.03 cal./deg. 

Herzberg and Spinks" have recently given a 
different interpretation to the overtone and 
combination bands in the spectrum; their assign- 
ment is in some respects more satisfactory than 
Sutherland’s (particularly as regards the valence 
overtones and combinations) but leads to sur- 
prisingly large values for the anharmonic term 
for vy; (—26.3 cm) and for the cross term x4; 
(—18.8 cm=). If, however, Sutherland’s assign- 
ment of the band at 2683 cm™ as v3—y»; be re- 
tained (which involves taking the band at 2670 
cm as vy ++v;—v;) and the remainder of 
Herzberg and Spink’s Table V be accepted un- 
altered, the expression for the vibrational energy 
becomes 


ey /hc= 19749, + 3358.402+ 13.622 + 3318.403 
— 30.40.?+ 729.3v,+ 609.70; — 4.705? — 4.9u;02 


4.7241 


1.04567 
0.0015 0 


0.86661 
0 


1.0907 
1.451 


' 4.60 


— 4.9003 — 1.30404 — 26.2010; — 112.1 2903 
— 9. 3v0v4 — 1.20905 + 5.0v90; — 5.80405. 


This gives for v5, w=614.4 hc/kT=0.88008 and 
x= 4.7/614.4=0.00765. The tabulated quantities 
are 


v1 v2 v3 v4 v5 
0.0610 0.0080 0.0090 0.8671 1.1137 
v 0.063 0.008 0.009 1.086 1.534 
v 0.07 0.01 0.01 2.85 5.24 


In Qy +e0/kT 
v 


With the same rotational constants as before, Eq. 
(13) gives (In Qrv+E)/RT)=7.7615, as com- 
pared with the exact value found by summation, 
7.7613. The cross terms x;; and their powers alto- 
gether contribute +0.01853 to f; of this terms of 
the type yi0j; amount to +0.01764; (445)? 
v4°052/2= +0.0052, (y15)?-v12v52/2=0.00026 and 
¥45¥15010s052= 0.00011. The rotational-vibrational 
interaction terms a0,- « -a;5;amountto — 0.00306, 
so that f is +0.00996. 

In conclusion, it should be emphasized that for 
the five gases so far considered (HzO, COz, N2O, 
HCN, C;H2) the rigid rotator-harmonic oscillator 
approximation will yield values of (F°—E,°)/T 
even for temperatures as high as 1000—-1500°K 
accurate to better than 0.1 cal./deg., and in four 
of the five cases to better than 0.05 cal./deg. 
This agreement is to some extent due to a partial 
cancelation of the terms neglected, and while one 
cannot assert that such cancelation will necessa- 
rily occur, one would expect it to be the rule 
rather than the exception. This assumption, if 
verified by future calculations for other mole- 
cules, should be a matter of some importance to 
those interested in problems of chemical equilib- 
rium; it implies that a knowledge of the moment 
(or moments) of inertia and of the fundamental 
frequencies is sufficient to permit the calculation 
of equilibrium constants for the temperature 
range of greatest experimental interest to better 
than +5 percent. 
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(1) The reasons for various failures to observe the 
absorption spectrum of free hydroxyl radicals are dis- 
cussed. (2) The absorption spectrum of OH in the electric 
discharge through water vapor is observed with a spec- 
trograph of high resolving power. After interrupting the 
discharge, the presence of OH is traced by its absorption 
spectrum over time intervals up to 1/8 sec. (3) With alter- 
nating half-periods of 60 cycle a.c. cut out by a switch, the 
emission of OH vanishes almost entirely with the current 
whereas the absorption persists rather uniformly. It is 
concluded that the persistence of the absorption spectrum 
is to be attributed to a long life of OH radicals rather than a 
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long life of O and H atoms furnishing a fresh supply of short 
life OH radicals. (4) The rotational intensity distribution 
within the OH band at strong currents (being anomalous 
to a lesser extent than at faint currents) indicates the 
presence of free OH in the discharge. (5) For a comparison 
of the present result with previous estimates, leading to a 
much shorter life of free hydroxyl, the number per sec. of 
emission processes of the OH band is measured, assuming 
that each of these processes leads to the production of one 
free OH radical. (6) The process by which OH radicals 
are consumed is discussed. 








I. PROBLEM AND PREVIOUS WORK 


HE fact has been emphasized by H. N. 
Russell,! investigating absorption spectra of 
stellar atmospheres, that in the visible and near 
ultraviolet absorption spectra are to be observed 
preferably of unstable molecules, for example, 
the free radicals OH, NH, CN, whereas most 
stable molecules, such as He, No, Os, COs, HO, 
have absorption spectra in the extreme ultra- 
violet. The reason is that most stable molecules, 
but not the unstable ones, have electronic con- 
figurations similar to those of the rare gases; 
hence they have high values of excitation energy. 
This preferred observation of free radicals by 
absorption spectra contrasts with the methods 
of chemical analysis. Although they have been 
well developed for stable molecules, there is no 
direct chemical test for the simplest free radicals 
like OH or NH. These radicals, however, play a 
major hypothetical part in theories of chemical 
reactions. Here evidently is a chance for the 
spectroscopist to supplement the chemical analy- 
sis of gases, that is, to observe radicals during 
the progress of a reaction. Although the absorp- 
tion spectra of most simple radicals have not 
been observed as yet, their wavelengths can be 
predicted since their emission spectra are known 
as well as the spectra of stable molecules. 
The OH group and its reactions present a 
problem of major interest because of the part 


1 Lecture at the Harvard Observatory, 1931; Astrophys. 
J. 79, 317 (1934). 
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they play in the hydrogen oxygen reaction.’ 
Furthermore, there is a particular interest in 
recent results of Rodebush and Wahl’ who sug- 
gest a bimolecular association OH+OH—-H.0O, 
as the principal process by which the OH 
radicals are consumed. 

In the present paper the life of OH radicals, 
produced by an electric discharge through H.O 
vapor as they disappear by a chemical reaction 
after the interruption of the discharge, is being 
investigated. The concentration of OH is esti- 
mated by the intensity of the absorption 
spectrum. 

Quantitative results regarding the concentra- 
tion can be derived from absorption, not from 
emission spectra. The reason is that within 
certain limits, in particular at low pressures, the 
absorption coefficient—not the intensity of the 
absorption line measured against the background 
—is proportional to the concentration of the 
absorbing molecules (Beer’s law)* whereas the 
intensity of the emission spectrum is very largely 
determined by other uncontrollable factors.’ 


2 Cf. K. F. Bonhoeffer and P. Harteck, Grundlagen der 
Photochemie, p. 257. 

3 W.H. Rodebush and M. H. Wahl, J. Chem. Phys. 1, 
696 (1933); cf. G. I. Lavin and F. B. Stewart, Proc. Nat. 
Acad. Sci. 15, 829 (1929). 

4R. W. Wood, Physical Optics, 3d ed. p. 101; Macmillan 
Company (1934). 

5 The method of Bay and Steiner (Zeits. f. physik. 
Chemie 3, 149 (1929)), based on the intensity of emission 
spectra, does not apply to the OH radical, since in the 
electric discharges the OH’ is largely generated directly from 
H.O so that the intensity of the OH band fails to indicate 
the concentration of OH. 
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Hence by photometry the absorption spectrum 
gives evidence of the concentration of the ab- 
sorbing gas, but—except after special calibra- 
tion—only of relative values, since an unknown 
probability of transition goes in determining the 
absolute values. It is intended to calibrate the 
absorption spectra thus obtained with a known 
partial pressure of free OH radicals and so to 
proceed to the measurement of the absolute 
value of their concentration. This will be 
attempted by thermal dissociation of HeO vapor 
according to Bonhoeffer and Reichardt.*® 

For rapid reactions the absorption spectra 
have to be taken in snapshots. In order to secure 
sufficient exposure, the chemical process and the 
snapshot have to be repeated periodically. 

In recent years there have been reports of 
numerous attempts to observe the absorption 
spectrum of OH in chemical reactions or electric 
discharges." These attempts, however, have 
failed,* although in the electric discharge through 
water vapor the emission spectrum, connected 
with the normal state of the OH molecule, 
appears very strongly; sO strongly, as a matter 
of fact, that the emission of OH forms one of 
the most annoying spectral impurities. The 
absorption spectrum of OH has been observed, 
however, in thermally dissociated water vapor at 
temperatures up to 1600°C by Bonhoeffer and 
Reichardt.* The complete failure to observe the 
absorption spectrum of OH in chemical reactions 
or electric discharges led to the conclusion that, 
in any case, the concentration of OH is extremely 
small since the high reactivity limits the lives 
of these radicals to a very short period. Bon- 
hoeffer and Reichardt were able to establish an 
upper limit for the concentration of OH present 
in the electric discharge through water vapor. 

It seems probable that all failures to observe 


°K. F. Bonhoeffer and H. Reichardt, Zeits. f. phvsik. 
Chemie 139, 75 (1928). 

"H. C. Urey and G. I. Lavin, J. Am. Chem. Soc. 51, 
3293 (1929); K. F. Bonhoeffer and T. G. Pearson, Zeits. f. 
physik. Chemie 14, 1 (1931); W. Frankenburger and H. 
Klinkhardt, Zeits. f. physik. Chemie 15, 440 (1932); 
P. Harteck, Trans. Faraday Soc. 30, 139 (1934); K. H. 
a P. Harteck, Zeits. f. physik. Chemie A170, 1 

3+) 

_ > The absorption spectrum of OH, occasionally observed 
in the under water spark by H. Stiicklen (Zeits. f. Physik 
30, 27 (1924)) is an exception. 

_,. K. F. Bonhoeffer and H. Reichardt, Zeits. f. physik. 
Chemie 139, 75 (1928): cf. E. Hulthen and R. V. Zumstein, 
Phys. Rev. 28, 14 (1926). 
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the absorption spectrum of OH in reactions or 
discharges have in common that the experiments 
have been undertaken with medium sized spec- 
troscopic instruments. To observe faint and 
sharp absorption lines, however, is much more 
difficult than is the corresponding problem for 
emission. A faint, sharp emission line, when 
photographed with an instrument like a medium 
size quartz spectrograph, appears broader than 
it actually is emitted from the source. A faint, 
sharp absorption line, however, when photo- 
graphed with the same instrument, disappears 
entirely, since its observation is about equivalent 
to the resolution of two closely adjacent emission 
lines.'!° For the most sensitive investigation of 
absorption spectra a spectrograph is required 
with a resolving power comparable with the 
width of the absorption lines to be observed. 
The great increase of sensitivity obtained by 
observing absorption lines with high resolving 
power becomes evident by the following com- 
parison: While the quartz spectrograph fails to 
show any OH absorption spectrum in the electric 
discharge through HO vapor with the carbon 
arc furnishing the continuous background, the 
21-foot grating reveals an undesirable OH ab- 
sorption spectrum even without the discharge 
through H.O in any carbon arc in air although 
operated in carefully dried air with purest 
carbons just baked out. 

There is an essential difference between chem- 
ical processes in the discharge and other chemical 
processes in that the discharge produces excited 
and ionized states of all possible energy values, 
some of them much higher than the energies 
available in most chemical reactions. It can be 
assumed, however, that the unstable excited 
states of atoms and molecules radiate within so 
short a time interval that they are not to be 
observed after interrupting the discharge with 
the rotating disk method. Jons and electrons 
recombine very probably at such a high rate 
because of the Coulomb force of attraction that 
they may be just barely observable with the 
same method.!! The recombination of atoms and 
radicals is probably by far the slowest process 


1 Cf. W.H. J. Childs and R. Mecke, Zeits. f. Physik 68, 
353 (1931). J 
"C, Kenty, Phys. Rev. 32, 624 (1928); L. 


J. Hayner, 
Zeits. f. Physik 35, 365 (1925). 
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after interrupting the discharge. This chemical 
process is presumably the only one actually 
observed in the present experiment. There might 
be exceptional cases, however, in which this 
chemical process, following an electric discharge, 
would differ from the ordinary chemical processes 
if metastable states of a considerable lifetime are 
produced. Such states do exist in very pure 
helium and very pure mercury vapor. Since 
small traces of impurities, in particular hydrogen, 
are very efficient in destroying metastable levels, 
it is expected that the metastables do not play 
a part for any appreciable period after the 
interruption of the discharge through H.O. 
Hence the results of the present method are 
directly comparable with the results of chemical 
methods for the investigation of the rate of 
chemical reactions in gases. 

Since the rate of reaction is strongly affected 
by the temperature, it is of interest that the 
temperature—also a rapid change of tempera- 
ture—can be observed, as well as the concentra- 
tion from absorption spectra, by measuring the 
intensity distribution within one branch of a 
band. 

The advantage of the method of absorption 
spectra over those of chemical analyses is its 
power to recognize free radicals and trace the 
progress of a fast reaction by snapshots. 


II. PERSISTENCE OF OH ABSORPTION SPECTRUM 


a. Apparatus 


The apparatus consisted of an electric dis- 
charge through water vapor, the absorption 
spectrum of which is to be observed; a source 
for the continuous spectrum; a synchronous 
motor, operating the electric switch and the 
optical shutter; and the spectrograph (Fig. 1). 

The discharge tube with large nickel electrodes 
and quartz windows, containing water vapor, 
afforded a column of 1.27 m length and 5 cm 
diameter for the absorption experiment. The 
electrodes, fastened in side tubes, had a distance 
of 10 cm from the main body of the tube which 
served as the absorption tube. Before the final 
experiments the glass walls were carefully freed 
of sputtered metal. ‘‘Conductivity water,’’ dis- 
tilled,several times in high vacuum, was used. 
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Fic. 1. Absorption experiment. A, carbon arc; L, quartz 
lens; D, discharge tubes; M, synchronous motor; X&, re- 
volving switch; D, sector disk; S, slit of spectrograph. 


The source for the continuous spectrum offered 
an unexpectedly difficult problem. At 3100A the 
hydrogen discharge tube, emitting a perfectly 
uniform spectrum, gave an intensity, even with 
the highest current density obtainable, far in- 
ferior to the carbon arc. The carbon arc in air, 
however, was useless since it contributed a 
pronounced absorption spectrum of OH. By 
various experiments it was found that this 
absorption spectrum was due to the combustion 
in air of hydrogen compounds serving as binding 
material for the carbons. The carbon arc in 
flowing nitrogen, carefully dried, was applied. 
Its spectrum, however, was not quite continuous 
because of superimposed nitrogen bands, nor 
was its intensity quite uniform with time because 
of the irregular shaping of the carbon electrodes. 

A synchronous motor carried on its axis a 
brass and Bakelite sector disk as a switch. It 
connected the discharge tube with the secondary 
of a transformer in most experiments for a half- 
period of 60 cycles a.c. and disconnected the 
tube for the next half-period or three half- 
periods. During this time interval the chemical 
recombination was investigated by the absorp- 
tion spectrum. In order to determine the maxi- 
mum time interval through which the absorption 
spectrum of OH could be traced, the speed of 
the revolving switch was reduced by gears to 
one-fifth of its former rate, the switch trans- 
mitting one full period of current and cutting 
out nine full periods. 

An optical sector disk of Bakelite rotating on 
the same axis with the switch, immediately in 
front of the slit of the spectrograph, screened off 
the inevitable radiation emitted from the water 
vapor during the electric discharge running and 
let through only emission flashes of the con- 
tinuous spectrum passing the water vapor within 
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short time intervals after interrupting the dis- 
charge. By varying the phase difference between 
electrical and optical sectors and observing the 
intensity of the absorption spectrum, decaying 
after the discharge, the lifetime of the OH 
molecules was investigated. With an optical 
sector of 12°, each individual flash lasted 1/900 
sec. 

A 21-foot Rowland concave grating servéd as a 
spectrograph. Its mounting had recently been 
rebuilt in a constant temperature room in the 
basement of the new Research Laboratory of 
Physics by Professor F. H. Crawford and Mr. 
I). W. Mann. Its second order at 3100A has a 
resolving power of 0.014A. On the other hand, 
the width of OH lines, when computed on the 
basis of the Doppler broadening at 400°C, is only 
(.005A, considerably smaller than the resolving 
power. Hence an apparatus with higher resolving 
power would afford a still more sensitive spec- 
troscopic test for OH radicals. In order to secure 
the highest resolving power a narrow slit and 
the fine grain Eastman IV —O plates were used 
in spite of their very small speed. 

Since an intermittent continuous spectrum 
was to be photographed with very high dis- 
persion on slow photographic plates, a serious 
difficulty was expected from the lack of light 
intensity. Hence the experiment started from a 
systematic comparison of the various methods by 
which the intensity of illumination could be 
increased, in particular by compensating for the 
astigmatism. The results were published in a 
previous paper.” A cylindrical lens in front of 
the grating proved useful.'* The times of ex- 
posure, including intermissions, were of the order 
of magnitude of 30 minutes. . 


b. Results 


By using the apparatus described, the absorp- 
tion spectrum of OH radicals came out clearly. 
(Fis. 2.) It was strongest when the water vapor 
pressure was high (2 mm) and the current strong. 
Beciuse of the inevitable intense emission spec- 
trun of the discharge, it was best observed 
imn ediately after the interruption of the current. 


. Oldenberg, J. Opt. Soc. Am. 22, 441 (1932); cf. G. 
J. Opt. Soc. Am. 23, 274 (1933) and J. B. Green and 
A. Loring, J. Opt. Soc. Am. 24, 348 (1934). 

\\. J. Humphreys, Astrophys. J. 18, 324 (1993); E. 


Diel 


Gehi ke, Zeits. f. Instr. Kde. 31, 217 (1911). 
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It could be traced down to about 0.3 mm. At a 
pressure of 2 mm, by means of the slow inter- 
rupter, its decaying intensity could be followed 
after the interruption over the full period of 
1/8 sec. Within this interval a considerable decay 
was observed (Fig. 3). The figure, however, does 
not express an average lifetime of OH radicals 
since it is largely dependent upon the sensitivity 
of the test. The time interval of 1/8 sec. holds 
for favorable conditions. 

In order to compare the absorption spectra 
during the discharge and after the interruption, 
care must be taken that during the discharge 
the emission spectrum is faint as compared with 
the continuous background or else it would 
mask the absorption spectrum. Hence a small 
current (30 ma) through the water vapor 
and, in order to secure sufficient absorption, 
rather high pressure (2 mm) was applied. The 
intensity of absorption was not much stronger 
during the discharge running than after the 
interruption. Hardly any change was noticed 
within one half-period—in agreement with the 
total observed as 
amounting to many half-periods of 60-cycle 
current. 

The intensity distribution within the absorp- 
tion band was very different from the one within 
the well-known OH emission band. The strong 
heads, characteristic in emission, disappeared in 
absorption. The reason is that in absorption the 
molecular rotation is low, mainly determined by 
the temperature, whereas in emission much 
higher rotational quanta come out with great 
intensity. The process exciting the abnormal 
rotation in emission has been discussed in a 
previous paper." 

Up to the present no quantitative photometric 
work, comparing concentrations of OH, could be 
carried through since the source of the continuous 
spectrum, the carbon arc in Ne, contributed 
traces of OH absorption spectra in a rather 
erratic way. It is easy to recognize the OH 
absorption due to the discharge which is unlike 
the one caused by the arc because of the widely 
different widths of the lines and intensity dis- 
tributions of rotation. For photometric work, 
however, a perfectly uniform background is 


persistence of absorption 


1440, Oldenberg, Phys. Rev. 46, 210 (1934). 
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Fic. 2. Spectra of OH radicals (enlargement 9-fold); part of band 3064; (1) emission; (2) absorption. 
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Fic. 3. Decay of absorption spectrum after interruption of 
discharge; (1) 1/300 sec.; (2) 1/8 sec. after interruption. 


indispensible. Therefore, only the estimate given 
above of the persistence of the OH absorption 
spectrum can be made as yet. 


III. NATURE OF THE LONG LIFE PARTICLES 


As described in the preceding section, the 
absorption spectrum of OH radicals can be 
traced over a period as long as 1/8 sec. after the 
discharge is interrupted. This, however, does not 
necessarily indicate a correspondingly long life 
for the free OH radicals. Bonhoeffer and Pearson, 
discussing the afterglow, suggested an alter- 
native; namely, that the long lives of the O and 
II atoms by slow recombination produce a fresh 
supply of OH radicals over a considerable period 
of time; these radicals are assumed to have a 
much shorter life. This is analogous to the per- 
sistence of the molecular emission spectrum in 
the well-known afterglow of nitrogen which is 
ascribed not to the long life of the molecules but 
to slowly recombining atoms." 


a. Lifetime not attributed to the recombination of 
O and H atoms 

First, we shall compare the various possibilities 
of slow recombination of O and H atoms with 
the spectroscopic observations. 

(1) Triple collisions forming OII'. Bonhoeffer 
and Pearson investigated various processes by 
which the band spectrum of OH might be 
excited. They proved by experiments that 7 
the discharge the emission of OH is largely due to 


15 G, Cario and J. Kaplan, Zeits. f. Physik 58, 769 (1929). 


dissociation of H2O and excitation of OH in one 
process (impact by fast electrons of 9 volts or 
more). This process has recently been confirmed 
by an independent argument." The afterglow of 
the OH bands, however, discovered by Lavin 
and Stewart!’ must originate from another ex- 
citation process since outside of the discharge no 
fast electrons are present. Bonhoeffer and 
Pearson proved that H atoms in Oz or O atoms 
in He fail to excite the OH band but that H and 
O atoms when mixed produce the OH band with 
high intensity. They concluded that the after- 
glow is caused by O+H-—-OH’, the energy of 
combination changing largely to excitation of 
OH; for known reasons triple collisions have to 
be assumed. Furthermore, they assumed a rela- 
tively short life of the OH radicals, estimated to 
10°? sec. 

This theory, which well explains previous 
results, is to be correlated with the absorption 
spectra described. The underlying hypothesis 
(short life of OH, long lives of O and H forming 
a fresh supply of OH) leads us to predict that any 
change in the supply formed per sec. will be 
followed by a change in the concentration of OH 
within a time limited approximately by the life 
of OH(~1/1000 sec.). Therefore a somewhat 
slower fluctuation of the supply ought to be 
accompanied by a corresponding fluctuation of 
the concentration. The cutting out of alternate 
half-periods of 60 cycle a.c. will cause a fluctua- 
tion of the discharge by which this theoretical 
prediction will be tested. By spectra taken in 
snapshots we can trace both the supply formed 
per sec. and the concentration. The supply is to 
be traced since all formation processes mentioned 
in the theory lead to excited OH radicals (during 
the discharge running by dissociation of |1,0 
into OH’+H as the process of major probal) lity 


' O. Oldenberg, Phys. Rev. 46, 210 (1934). - 
17 G.I. Lavin and F. B. Stewart, Proc. Nat. Acad. Sc . 15 
829 (1929). 
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and as afterglow between discharges by O--H 
+M—OH’+M); hence each formation process 
of an OH radical is accompanied by an emission 
process of the OH band. Therefore, by measuring 
the intensity of the emission spectrum we 
actually register the supply formed per sec. of 
OH radicals. On the other hand, independently 
of any theory, by the absorption spectrum we 
can also observe the concentration of OH present. 
The theory under discussion predicts that the 
concentration of OH will fluctuate at the same 
rate as supply formed per sec. This becomes 
equivalent to the absorption spectrum of OH 
fluctuating at the same rate as the emission 
spectrum. This prediction is to be tested by 
snapshots of emission and absorption spectra 
during and between discharges. 

With the same discharge through water vapor 
and the same sector disk and interrupter, letting 
through and cutting out alternate half-periods, 
and with a quartz spectrograph, the emission of 
OH was compared during and between dis- 
charges. It faded to an inferior order of magni- 
tude after the interruption. For example, ap- 
proximately the same photographic effect was 
produced during the discharge in 1/2 sec. and 
between discharges in 5 hours, the ratio being 
about 1 : 40,000. Therefore according to the 
theory under discussion—concentration rapidly 
changing with supply formed per second—we 
expect a striking decay of the concentration 
after the interruption. The absorption spectrum 
of OH, however, indicating the concentration, 
persisted during the intermission of one half- 
period without any striking decay as discussed 
above (II 2). This contradicts the theoretical 
prediction based on a very short life of OH, but 
long lives of O and H. 

’) Triple collisions forming OI. Secondly, we 
must consider O+H combining in triple collisions 
leading to OH instead of OH’. (In this case we 
cannot estimate the supply formed per sec. of OH 
from the intensity of emission.) There are two 
arguments against this process. Bonhoeffer and 
Person actually proved that two streams of 
gas’s, one containing H atoms, the other O 
atons when flowing together, emit OH bands 
wit: high intensity. Hence the very small in- 
ten-ity of these bands, observed during the inter- 
mis-ion, indicated that no large concentration of 
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both atoms is present. Hence we fail to explain 
the absorption between discharges by the forma- 
tion of OH from O and H. The other argument 
leading to the same conclusion is as follows: 
Granted that some OH radicals are supplied 
during and between discharges by recombination 
of O and H atoms, we still positively know that 
during the discharge another powerful source of 
OH radicals exists, the simultaneous dissociation 
of H2O and excitation of OH by electron impact. 
Hence we should expect the concentration of OH 
to be considerably stronger during than after 
discharges—contrary to the observation. There- 
fore, the long lifetime of 1/85 sec. cannot be 
attributed to O and H atoms forming short life 
OH radicals. This argument applies against the 
catalytic effect of the wall as well as against 
triple collisions in the gas. 


b. Lifetime attributed to OH radicals 


Since the formation of OH from the atoms is 
excluded by the arguments just given, the 
plausible conclusion is that the persistence of the 
absorption spectrum of OH observed is due to a 
correspondingl:, long life of OH radicals. 

Another argument against a very short life- 
time of OH radicals is given by the fluorescence 
radiation of H2O vapor, sensitized by Hg vapor. 
When excited by the Hg resonance line 2537 the 
OH band 3064 is strongly emitted.'* The energy 
of 4.9 volt electrons, corresponding to the line 
2537, is unable to dissociate the H»O molecule 
and excite the OH radical in one process; hence 
we must assume two consecutive processes. It is 
most plausible to assume one impact of the 
second kind leading to dissociation of HO into 
OH-+H and another to excitation of OH. The 
second process is probable only if an appreciable 
concentration of OH has been built up by the 
effect of the first process; or in other words, 
if the OH radicals have an appreciable lifetime. 


c. Contribution of decomposing H.O, 


In addition to the association of OH+OH, 
suggested by Rodebush and Wahl as the main 
reaction by which OH radicals disappear, one 


18 FE, Gaviola and R. W. Wood, Phil. Mag. 6, 1202 (1928). 
Under somewhat different conditions (a few mm H,O 
vapor, sensitized fluorescence), Gaviola and Wood est!- 
mated the lifetime of free OH to the order of magnitude of 
1/100 sec. 
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must consider the effect of another process 
regenerating these radicals. According to Geib 
H.O2 is decomposed by H atoms probably in 
the following process H+H:0.—-H.2O+0OH.’® Is 
there a chance that the persistence of the absorp- 
tion spectrum is caused by the slow liberation of 
OH radicals in this process from H2Os, while 
the radicals themselves have a much shorter 
life? Since so large a concentration of HsO» was 
observed by Rodebush and Wahl among the 
gases pumped out of the discharge, it is hard to 
assume that a still much larger amount of H2O. 
was originally produced which would give rise 
to the full amount of OH observed in the absorp- 
tion spectrum. Moreover it is hard to assume that 
OH is stored up and bound in the shape of H2O» 
from which it is gradually liberated only for very 
brief moments by H atoms, and yet the perma- 
nent concentration of free OH is so high, over a 
considerable period of time, that its absorption 
spectrum is well observable. It is probable, 
therefore, that the lifetime of OH observed in the 
absorption spectrum is largely determined by the 
real life of OH slowly recombining. 


IV. ESTIMATE OF CONCENTRATION OF OH FROM 
THE EMISSION SPECTRUM 


The absorption spectra described indicate only 
relative values of the concentration of OH under 
different conditions. They fail, however, to give 
evidence of the OH concentration as compared 
with the concentration of H:O molecules. For 
the stationary concentration caused by a heavy 
current, the emission spectrum indicates an 
appreciable concentration of OH in the following 
experiment. 

In a previous paper it was pointed out that 
the excitation process (simultaneous dissociation 
of H.O and excitation of OH) leads to a charac- 
teristic abnormal rotation of the radiating OH, 
observed by the intensity distribution in the 
band spectrum. Conversely, we can recognize 


19 K. H. Geib, Zeits. f. physik. Chemie A169, 166 (1934). 
Mr. N. D. Smith called my attention to the possible effect 
of this reaction on the lifetime of OH. The next reaction, 
suggested by Geib, OH+OH—H.0+0 does not agree 
with the results of Rodebush and Wahl who proved the 
absence of atomic O from the discharge through H.O 
vapor (cf. Section Vla). The alternative OH+OH+M 
—H.02+M seems preferable. Since, however, this reaction 
regenerates half the H,O2, just decomposed, the yield of the 
first reaction would be higher than computed Dy Geib. 


Fic. 4. Part of OH band (negative, quartz spectrograph, 
enlargement 55-fold). (1) weak current; (2) strong current. 
Abnormal rotation is pronounced at weak current; the 
following lines are enhanced; a=R,(173) and R,.(153); 
b= R,(143); c= R2(134); rotational quantum number 3 or 4 
would be most probable according to the temperature. 


this excitation process from the strongly abnormal 
intensity distribution observed. On the other 
hand, free OH radicals in the discharge tube are 
certain to contribute a band of normal intensity 
distribution in emission since the absorption 
spectrum reveals normal rotation of the un- 
excited OH radicals present. (Section II b.) 
Therefore we may expect for a very small current 
(hardly any free OH present, hence all excited 
OH originating directly from H2O) the greatest 
abnormal rotation and, on the other hand, for a 
very strong current a tendency towards normal 
rotation (smaller intensity of the very high 
rotational quanta) provided that an appreciable 
part of the intensity is due to free OH radicals 
present under the action of the current. 

For a test the intensity distributions are com- 
pared for currents of 3 ma and 250 ma by spectra 
taken on the same plate, and with the same 
times of exposure in order to make them com- 
parable, the intensity with the heavy current 
being reduced by dense gauze screens. Actually, 
as expected, heavy currents tend to reduce the 
high rotational lines. This effect can best 
observed on lines of an R branch which hev 
neighboring wavelengths although they bel: 
to widely different rotational quanta (Fig. 
This observation proves that for strong curr 
the partial pressure of OH radicals is appreci« 
Since in the discharge the OH radical is prefe: 
by its small energy of excitation (4.0 volt 
small concentration of free OH would be + 
cient for the observed effect. In this argumen 
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unknown contribution of H2O: to the spectrum is 
neglected. 


V. THE RATE OF PRODUCTION OF OH BY THE 
Process H.O+el—-OH’+H 


It might seem unsatisfactory to attribute the 
persistence of the OH absorption, observed as 
1/8 sec., to the long life of OH radicals (Section 
IIIb) since all previous papers agree to a very 
short life of OH, estimated to 1/1000 sec. The 
reason for this discrepancy must be discussed.?° 

The main argument for a short life has been 
that the striking intensity of the emission band 
3064 indicates a supply of OH formed per sec. so 
rich that its absorption spectrum could be 
expected to be observable with the quartz spec- 
trograph (Section I).”! It is difficult, however, to 
base this estimate of intensity upon a correct 
comparison. The observation shows only that the 
intensity of 3064 is strong when compared to 
other spectra. 

The following comparison seems preferable: 
Assuming that each emission process of the band 
3064 indicates the production of one OH radical 
and furthermore, assuming a lifetime as long as 
1/8 sec. (Section IIb), the intensity of the 
emission band, measured in absolute units, leads 
to the concentration of OH (=number produced 
per sec. Xlifetime) as far as it is produced by 
this observed process. This concentration is to be 
compared with a reasonable minimum concentra- 
tion required for an absorption spectrum. The 
previous estimates (leading to a short life of OH) 
make us believe that this measurement of supply 
combined with a Jong lifetime of OH would lead 
to a concentration so high as to render the 
previous failures to observe-the absorption spec- 
trum inexplicable. 

For this estimate the total number per sec. of 
emission processes of the band in the whole dis- 


_.*° The estimate of Bonhoeffer and Pearson, based on the 
idea that a considerable fraction of the total electric energy 
serves for producing OH’, need not be discussed here since 
it leads only to a lower limit of the lifetime of OH. 

*1 From the failure to observe the absorption spectrum 
Bonhoeffer and Pearson derived the concentration of OH 
as being less than 0.03 mm. It must be taken into account 
that the calibration of an absorption spectrum at four 
times the temperature and about five hundred times the 
pressure of the discharge cannot be quite accurate unless 
the spectrograph has very large resolving power, since in a 
medium size spectrograph a narrow line has a much larger 
tendency to escape from observation than a broad one. 
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charge must be measured. The densitometer 
trace of the band—with proper calibration of the 
photographic plate—yields a curve representing 
intensity in arbitrary units as a function of wave- 
lengths. The emissive power of the surface of 
the discharge is obtained in absolute units from 
a quantitative comparison with the black body 
radiation at the same wavelength.” (In order to 
take the spectra of the discharge and the black 
body at 1415°K on the same plate with the same 
time of exposure, the intensity emitted from the 
discharge was reduced by the factor 950.) By 
integration over the whole surface of the dis- 
charge the total energy or number of quanta 
emitted per sec. from the whole discharge can be 
derived; that is, the total number of OH radicals 
generated per sec. by this observed process, the 
yield of which is to be estimated. 

With a lifetime of 1/8 sec. a concentration of 
410" OH radicals per cm® results, or of 1 OH 
radical, as produced by the process considered, 
among 8000 H.O molecules. (Apart from the 
uncertainty of the estimate caused in particular 
by the rather uncertain lifetime, it is easy to 
explain why the total concentration of OH must 
be considerably higher than the contribution of 
one process estimated here.) The result is that 
the intensity of the OH band observed in absolute 
units combined with a lifetime as long as 1/8 sec. 
does not lead to absurdly high values of the con- 
centration of OH. 


VI. PROCESSES IN THE DISCHARGE THROUGH 
H.O Vapor 


a. Reactions consuming OH radicals 


The reactions by which OH radicals disappear 
have been much discussed. It seems that the 
chemical processes are simpler in the electric 
discharge through water vapor than in the photo- 
chemical or thermal combination of He and Os. 
The reason is that in the discharge at moderate 
currents no free oxygen is formed. This is proved 
by the following three arguments: (1) Rodebush 
and Wahl proved by the Stern-Gerlach experi- 
ment that free O atoms are absent from the gases 
flowing out of the discharge. (2) They found no 
molecular oxygen in the discharge. (3) As re- 


22 | am obliged to Dr. F. F. Rieke for the construction of 
the black body. 
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ported in Section III, the afterglow was exceed- 
ingly faint (in agreement with the results of 
Rodebush and Wahl and of Stoddart).” Its 
intensity, however, is known to be high if both 
O and H are present. Since there is no doubt 
about the presence of H, it follows that O is 
absent." The conclusion agrees with the result 
of Stoddart who enhanced the afterglow by 
introducing Oz to the H2O vapor. Hence in the 
discharge through H2O vapor with moderate 
current we are allowed to discard the considera- 
tion of O or Ox. 

Among the gases pumped out of the discharge 
Rodebush and Wahl found a large concentration, 
about 50 percent, of H2O2. (Considering that 
H.O2 decomposes by the effect of the electric 
current as well as in the presence of H atoms, 
the concentration of HO. originally formed 
must have been still higher.) They explained 
previous failures to observe H2Oz by the catalytic 
effect of metal sputtered in the discharge tube, 
decomposing the H2Oe. After a detailed dis- 
cussion of all possible reactions leading to H2O» 
they concluded that it is formed from OH+OH 
and that this is the main reaction by which OH 
disappears. 

Harteck proved that the reaction by which OH 
disappears is homogeneous.” The rate of this 
reaction, however, seemed to bring in a difficulty. 


23 E. M. Stoddart, Phil. Mag. [7] 18, 410 (1934). 

24 The absence of free oxygen might be interpreted by 
the assumption that in the electric discharge at moderate 
energies of electrons, OH radicals (as well as N»2 molecules) 
have a tendency to become excited or ionized rather than 
dissociated. This is to be expected from the high relative 
intensity of the 0-0 band emitted by the neutral molecule 
and the small change of internuclear distance that takes 
place by ionization. The values of these distances were 
kindly communicated to me by Professor F. W. Loomis; 
they are for the neutral molecule 7’ = 1.022, 7’’=0.979 and 
for the ionized molecule r’=1.139 and r’’=1.027. The 
absence of free oxygen is not easy to reconcile with the ob- 
servation of G. von Elbe (J. Am. Chem. Soc. 55, 62 (1933)), 
who found free oxygen in HxO2+He2 when illuminated with 
the light of the zinc spark. It has been proved by Urey, 
Dawsey and Rice (J. Am. Chem. Soc. 51, 1371 (1929)) that 
the intense zinc lines around 2100A excite the OH band in 
the fluorescence of HO». Hence OH is certainly produced. 
Nevertheless, the interpretation of the experiment is not 
unique since the continuous absorption spectrum of H.Osz, 
observed between 2000 and 3100A is caused by the super- 
position of different processes; for instance light of 3100A, 
just coinciding with the band, clearly cannot produce dis- 
sociation of H,O2 and simultaneous excitation of OH which 
are produced by 2100A. The light of the zinc spark, there- 
fore, is certain to produce other additional particles, making 
the interpretation rather uncertain. 

°> Cf. K. F. Bonhoeffer and T. G. Pearson, Zeits. f. physik. 
Chemie 14, 5 (1931). 
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Rodebush and Wahl found the concentration of 
H.O2, pumped out of the discharge and con- 
densed in a liquid air trap, to be independent of 
the time interval of the flow between discharge 
and trap, even for intervals as short as 1/300 sec. 
They inferred that the process OH+OH—-H.0, 
is finished within a still shorter time. Considering 
this high rate of reaction, they discarded the idea 
of a triple collision and rather assumed a bi- 
molecular association.2*> The argument for this 
rather unusual process becomes questionable, 
however, when considering the longer lifetime of 
OH reported in Section IIb. Instead, the follow- 
ing alternative, which was mentioned by Rode- 
bush and Wahl themselves, seems preferable. 
It is assumed that the combination of OH+OH 
—H.Oz goes on in the gas as well as within the 
liquid air trap so that for the measurement of the 
H.O2 concentration it does not matter whether 
a molecule is formed before or during condensa- 
tion. Because of the reported much _ longer 
life of OH there is no argument left in favor of a 
bimolecular association. Consequently, the usual 
assumption of a triple collision for the process 
within the gas seems preferable. The experiments 
reported above are not sufficiently quantitative 
as yet to give direct evidence of the order of the 
reaction by which OH disappears. 

More questionable than the homogeneous re- 
actions discussed are the reactions of OH on the 
wall. Urey and Lavin*’ discovered the high re- 
activity of H2O vapor when partly dissociated by 
an electric discharge. The present results agree 
with their assumption that OH radicals can be 
pumped out over considerable distances from 
the discharge. Yet it seems difficult to attribute 
the reactivity of the gas mixture to one or the 
other particle, for instance OH, formed in the dis- 
charge. According to Rodebush and Wahl it is 
true that in pure glass vessels only H, OH and 
H.Oz are produced by the discharge. But since 
H2O2 decomposes readily in the presence of a 


6 A triple collision is fundamentally required for the asso- 
ciation of two atoms (like H-++H); not, however, for the 
association of two radicals forming a fouratomic molecule 
(like OH+OH—H,0.). A process of this type without a 
triple collision was reported by Farkas and Harteck (K. F. 
Bonhoeffer and P. Harteck, Grundlagen der Photochemie, 
p. 257). The probability was estimated on the basis of wave 
mechanics by N. Rosen (J. Chem. Phys. 1, 319 (1933)). 

27H. C. Urey and G. |. Lavin, J. Am. Chem. Soc. 51, 
3290 (1929). 
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catalytic surface, oxygen also must be taken into 
consideration. There seems to be no chance to 
investigate individual effects of one or the other 
particle since Stoddart reports a failure to 
separate the hydroxyl from atomic hydrogen and 
oxygen. 


b. Survey of processes in the discharge through 
H.O vapor 

There is no doubt about the principal primary 
process taking place in the discharge, the water 
molecule being broken up by electron impact: 
H.O+electron—OH’+H. (Instead of OH, the 
H atom might become excited; or both particles 
might be formed in their normal states.?*) 
Since at moderate currents free atomic or 
molecular oxygen does not play a part, the 
processes to follow are the combinations of 
the particles mentioned, namely, H+H+M-—H, 
+M and OH+OH+M-—H,.0.+M. The follow- 
ing competing associations must be considered: 
H+OH+M-—H.0+M obviously takes place 
but at a lower rate, or else the concentration of 
HO» could not be so high as actually observed 
in glass tubes with no metal present. For the 
same reason the reaction OH+OH—-H.,0+0 is 
supposed to have a lower rate; moreover, it 
would lead to atomic oxygen, the absence of 
which has been proved. The H2Oe formed in the 
discharge is partly broken up again by H.Oz 
+H—-H,O+0OH. The reaction OH+H:—-H:,0 
+H is known to have an energy of activation 
=14 kcal. so that it does not take place at room 
temperature.”? 
A fundamental uncertainty is unavoidable in 


*S J. Franck, H. Kuhn and G. Rollefson (Zeits. f. Physik 
43, 155 (1927)) describe analogous processes caused by 
absorption of light. 

*Y W. Frankenburger and H. Klinkhardt, Zeits. f. physik. 
Chemie 15, 436 (1932). : 
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discussing processes in discharges. Excited atoms 
and molecules are being produced which, due to 
the high energy content of their electronic 
structure, might have a special reactivity, so high 
that even during the life of the excited state 
reactions are probable. Obviously the processes 
after the interruption of the discharge are easier 
to correlate to known chemical processes. 

In the combination of He and Oz, started by 
photochemically produced H atoms, H2Qz¢ is 
formed, too. If the most recent theory is correct, 
in this photochemical experiment the H»2Qz is 
formed over the intermediate compound HO: by 
association of H atoms with Os molecules.*° 
It is of interest that in the electric discharge 
through H.O vapor the formation of H2O, takes 
place by an entirely different process. 

The main results of the present paper are as 
follows: As a test for free OH radicals absorption 
spectra taken in snapshots were applied. They 
indicate in particular the change of concentra- 
tion with time. The lifetime of free OH radicals 
is much longer than previously assumed, so 
long that no bimolecular association must be 
assumed for the process by which they dis- 
appear. Attempts will be made to obtain quanti- 
tative results by photographic photometry and 
to search for free OH radicals in other reactions. 

The author wishes to express his sincere 
appreciation of the work done at the machine 
shop of this laboratory under the supervision of 
Mr. D. W. Mann, in particular by Mr. L. Cail, 
and of the glassblowing work of Mr. H. W. 
Leighton. He is greatly indebted to Professor G. 
B. Kistiakowsky and Dr. A. A. Frost for helpful 
discussions. 


3° K. H. Geib, Zeits. f. physik. Chemie A169, 161 (1934); 
K. F. Bonhoeffer and P. Harteck, Grundlagen der Photo- 
chemie, p. 262. 
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The Statistical Weights of the Rotational Levels of Polyatomic Molecules, Including 
Methane, Ammonia, Benzene, Cyclopropane and Ethylene 


E. Bricut WILSON, JR.,* Mallinckrodt Chemical Laboratory, Harvard University 
(Received March 5, 1935) 


A general method is described for calculating the statis- 
tical weights (degeneracies) of the energy levels of poly- 
atomic molecules. Wave functions for a molecule are as- 
sumed to be expressible as linear combinations of products 
of the electronic, vibrational, rotational and nuclear spin 
functions. By using standard methods of group theory, 
the number of linear combinations of these products are 
found having the correct symmetry with respect to those 
permutations of identical atoms which are equivalent to 
rotations of the molecule. It is not necessary to find the 


combinations themselves. The molecules CHy, CDs, CH;D, 
CHD;, CH2D,, CH;X, CD;X, NHs, NDs, CeHe, CsHe and 
C,H, are treated. In addition noncombining species in 
polyatomic molecules and the splitting of energy levels due 
to the multiplicity of equilibrium configurations are dis- 
cussed. Tables of statistical weights are given for the above 
molecules which could be used to interpret alternating 
intensities in rotation-vibration spectra or for more exact 
calculations of thermodynamic quantities than are usually 
made. 





EVERAL investigators! have calculated the 

quantum weights of the rotational states of a 
number of polyatomic molecules and the prin- 
ciples involved are well known. This paper 
indicates how such calculations can be con- 
siderably simplified (especially for complicated 
molecules) by the application of group theory. 
Methane is used as an illustration. 


CONSTRUCTION OF APPROXIMATE WAVE 
FUNCTIONS 


The wave function for a molecule is usually 
approximately expressed as a product of func- 
tions of the various coordinates; i.e., 


Y=Velvvbeyrys, (1) 


where the subscripts E, V, R, T and S indicate 
that the corresponding factor is a function of the 
electronic (including electron spin), vibrational, 
rotational, translational and nuclear spin coordi- 
nates. The electronic, vibrational and rotational 
wave functions depend on the choice of equi- 
librium configuration of the nuclei, regarded as 
distinguishable. Thus in defining the Eulerian 
angles 0, y, and x for Wr(8, ¢, x) it is necessary 
to define the moving Z’ axis in terms of certain 
specified nuclei. 


* Junior Fellow of the Society of Fellows. 

1F. Hund, Zeits. f. Physik 43, 788 (1927) ; 43, 805 (1927). 
These papers describe Hund’s general method and treat 
NH;; W. Elert, Zeits. f. Physik 51, 6 (1928). This paper 
treats CH, by Hund’s method which in this case is much 
more difficult than the group theory method; D. M. 
Dennison, Rev. Mod. Phys. 3, 280 (1931). 


For this reason the wave function of Eq. (1) 
is not unique; there will exist a number of others, 
obtained from it by permuting identical nuclei. 
Some of these permutations may not, however, 
yield independent wave functions, since they 
may correspond to rotation of the molecule and 
therefore can be expressed in terms of other 
wave functions of the set in a manner which 
will be discussed later. Fig. 1 illustrates the two 
independent configurations for methane. 

In addition to the degeneracy due to this cause, 
there may exist additional degeneracy due to the 
fact that yr, Ys and sometimes Wy and Wz may 
themselves be degenerate. Thus for CH, ws is 2! 
or 16-fold degenerate while ys for CD, is 34 or 81- 
fold degenerate. Moreover, Wz for a spherical top 
(such as methane) has a degeneracy (2/+1)? for 
the energy level with quantum number J. 

Although no one of the functions y of the 
type of Eq. (1) possesses such symmetry, the cor- 
rect wave function for the molecule must be 
either completely symmetrical to interchange of 
members of a given set of identical particles or 
it must be completely antisymmetrical to such an 
interchange. Since protons obey Pauli’s principle, 


@ @ 
aL.» Ss @ © ~e 
(3) (3) 
a b 


Fic. 1. The two different configurations or “frameworks” 
for methane. 
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CH, has antisymmetrical wave functions whereas 
CD, has symmetrical wave functions. In order 
to obtain a function with the proper symmetry 
character it is necessary to form linear combina- 
tions of the degenerate functions discussed above. 
The number of such combinations with the 
correct symmetry properties which can be formed 
is the quantum weight for that energy level. 

In forming these combinations we first com- 
bine those functions for a given energy which 
correspond to the same configuration of the 
nuclei; i.e., the same “framework” of atoms, 
say a of Fig. 1. We form the combinations so 
that each one has the correct symmetry with 
respect to those permutations of identical par- 
ticles which correspond to rotations of the 
molecule without change of framework. When we 
have done this for each framework, combinations 
involving functions based on different frame- 
works are made in such a manner that the 
resulting wave functions have the correct sym- 
metry with respect to all possible permutations 
of identical nuclei. 

In finding the number of combinations of 
functions based on one framework and possessing 
the correct symmetry properties, we use the 
methods of group theory.? The permutations of 
identical particles which are equivalent to rota- 
tions of the molecule form a group. The set of 
functions PeWvWrvs based on one framework and 
belonging to one energy level constitute the 
basis of a reducible representation of this group 
which, by the construction of linear combina- 
tions, can be reduced so that the combinations 
divide into sets each of which forms an irre- 
ducible representation of the group. Only one 
certain irreducible representation will possess the 
required symmetry properties; we therefore seek 
the number of combinations of the functions 
YvevvWevs which correspond to this particular 
irreducible representation. (In methane this is 
the representation A of Table I, since all the 
rotations are equivalent to even permutations.) 
‘It is this number which is given by the formal 
methods of group theory without the need of 


*See, for example, E. Wigner, Gruppentheorie, Vieweg 
and Sohn, 1931; B. L. van der Waerden, Die gruppentheore- 
tusche Methode in der Quantenmechanik, Springer, Berlin, 
1932; A. Speiser, Die Theorie der Gruppen von endlichen 
Ordnung, Springer, Berlin, 1927; Also R. S. Mulliken, 
Phys. Rev. 43, 279 (1933). 
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TABLE I. Irreducible representations of the group T. N=12. 








E 3C, 4C; 4C;’ 
1 





—1 


w =e2til3 








actually constructing the combinations them- 
selves. 

The work is further simplified by the fact that 
the symmetry of the set Pe-vrws can be ob- 
tained by finding the symmetry of yz, yv, ve 
and ws separately and then taking the ‘direct 
product” of their irreducible components. We 
shall therefore study each of these functions 
separately. 

The electronic wave function yz is ordinarily 
completely symmetrical for the normal state and 
the modifications in the method which are 
necessary in case this is not true will be obvious. 
The translational function yr does not involve 
internal coordinates and is therefore symmetrical. 


SYMMETRY OF THE ROTATIONAL FUNCTIONS 
FOR METHANE 


The functions ~re= Wsxmu(6, ¢, x) for methane 
are the solutions of the spherical top wave 
equation. The angles 0, g and x are the Eulerian 
angles, 6 and g being the ordinary polar co- 
ordinates of an arbitrarily chosen axis in the 
molecule. The group of permutations which 
correspond to rotations is in methane iso- 
morphous with the symmetry point group T, 
which consists of those pure rotations which 
carry methane into itself. The irreducible repre- 
sentations of T have the characters’ shown in 
Table I. 

The symbols E, C2, C2 refer to the operations 
identity, twofold rotation, and threefold rota- 
tion, respectively. C3’ differs from C3; in the 
sense of rotation. The symbols A, E,, E, and T 
represent the irreducible representations. Applied 
to Fig. 1a, one C2 is equivalent to the permuta- 
tion (12)(34) while one C; is equivalent to 
(234), both of which are even permutations. 

3 This table is taken from R. S. Mulliken, Phys. Rev. 43, 
279 (1933) which contains similar tables for all the 
crystallographic point groups; See also E. Wigner, Gott. 


Nachr., p. 133 (1930) and L. Tisza, Zeits. f. Physik 82, 48 
(1933). 
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TABLE II. Characters for the representation of T formed by Wr. 


(p, k=0, 1, 2, 


++, Multiply all characters by 2J+1.) 





, , 
XE Xo, 


, , 
x C3 x C3" 








J=6p+k | 12p+2k+1 (—1)é 


(2/¥ 3) sin 3r(1 —k) 2/3) sin 3x(1—R)| 





In order to reduce* the representation formed 
by the functions Yyxu(6, ¢, x) for a given value 
of J it is necessary to know how these functions 
transform under the operations C2 and C3. If we 
choose the Z’ axis in the molecule so that it 
coincides with one of the threefold axes, the 
operation C; will affect only the angle x, changing 
it to x +3. Since x enters Ysxu(O, ¢, x) only 
by way of the factor* e‘**, the transformation 
C; has the effect 


C3 
x)vsKxul9, 9, x+57) 
= et (2ri/KY 7 (0, ¢, x). (2) 


The character’ x’c, (which is the sum of the 
diagonal elements in the transformation matrix 
for C3) is therefore® 


vik (8, ¢, 


(2J+1) > ef ri/3)K 


Ka—J 


which has the values given in Table II. 

With the choice of axis used above it is awk- 
ward to obtain x’c,. However, we may make use 
of the fact that the choice of axis in a spherical 
top molecule is entirely arbitrary. Changing the 
axis is equivalent to a linear transformation of 
the functions Yyxu(6, ¢, x), with J and M fixed. 
Furthermore the diagonal sum of a matrix is 
invariant under such transformations. It is 
therefore convenient to calculate x’c, with the 
axis coinciding with C2, so that 


C2 
Yikul9, 9; x) WsKu(O, ¢, x +7) 


=e*riv xu (8, ¢, x), (3) 
which yields 


z 
c,=(2J+1) 2 (-1)*, 
K=a—J 

3a The results of this paragraph may also be obtained 
from C. Eckart, Rev. Mod. Phys. 2, 344 (1930) Eq. (111). 

4 Reiche and Rademaker, Zeits. f. Physik 39, 444 (1926). 

5 See reference 2. 

6 The factor (2J+1) comes from the 2/+1 values of M 
associated with a given J. 


whose values are also given in Table II. 

These characters and those in Table I are 
used in the standard group theory formula’ to 
obtain the symmetries of the rotational functions 
Wr. These are given in Table III. (Since E and 7 
TABLE III. 


Symmetry of the rotational functions for methane. 





J=6p 

J=6p+1 
J =6p+2 
J =6p+3 
J=6p+4 
J =6p+5 


(27-41) ‘(4 

(27 +1) 

(27 +1) p 

(27 +1) y, 

(2J+1) (p 
rb 


+1)E+(3p+1) 
p E+(3p4+2) 
+1)E+(3p42) 
isnt 








are degenerate representations their coefficients 
in Table III need to be multiplied by 2 and 3, 
respectively, to give the number of linear com- 
binations with these symmetries. ) 


SYMMETRY OF VIBRATIONAL FUNCTIONS 


A permutation of identical particles equivalent 
to a rotation of the molecule may affect Wy also. 
A simple illustration is shown in Fig. 2. Before 
interchange the distance (2-3) is greater than 
(1-3); after interchange these distances have 
exchanged magnitudes. But the normal coordi- 
nate for this motion is 


Q= (2-3) — (1-3), (4) 


so that interchanging the magnitudes of these 
distances changes the sign of Q. 

Since the wave function yy is a function of 
the normal vibrations Q;, the functions yy may 
be transformed when the Q,’s are. In particular, 
the wave functions for a state with one quantum 
of one normal vibration excited transform in the 
same manner as the corresponding normal Co- - 
ordinates Q for that vibration.”* It is known’ 

7 See reference 2. 

7a The wave functions corresponding to overtone and 
combination levels usually have transformation properties 
different from those possessed by the normal coordinates, 
but obtainable from them. See Tisza, reference 3. 


*For the application of group theory to molecular 
vibrations see E. Wigner, Gott. Nachr., p. 133 (1930); L. 
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Fic. 2. The effect of the permutation (12) on one of the 
normal coordinates of a linear triatomic molecule. 


LEVELS OF 


that the normal coordinates Q; and the vibra- 
tional wave functions yy have symmetry prop- 
erties corresponding to the various irreducible 
representations of the symmetry point group for 
the molecule. This point group may contain 
reflections, inversions, etc., besides pure rota- 
tions. We are only interested in the pure ro- 
tations, but the effect of these on yy may be 
easily obtained from a knowledge of the irre- 
ducible representation of the complete symmetry 
group to which yy belongs. 

Methane belongs to the point group® Ty, (of 
which T is the rotational subgroup). Using the 
notation of Mulliken,’® we can find! that there 
is one non-degenerate normal vibration (irre- 
ducible representation A, of Tz), one doubly- 
degenerate vibration (£), and two triply-de- 
generate vibrations (72). The fundamental levels 
of yy therefore have these symmetries. In this 
case it is very simple to find to what representa- 
tions of T these various representations of Ty, 
belong. One finds'® that A; and Az of Ty belong 
to A of T, E belongs to E= E,+ Ee, and 7; and 
T, belong to T. Therefore, knowing*® the sym- 
metry of yy with respect to Tu, we can find its 
symmetry with respect to T. 


SYMMETRY OF THE NUCLEAR SPIN FUNCTIONS 


The nuclear spin functions ~s form a basis of 
a representation of the group of permutations of 
identical atoms which are equivalent to rotations. 
For CH, there are’ 16 fynctions ws, typical 
members of which are listed in Table IV, which 
gives the quantum number m, for each atom. 
To obtain the characters for the representation 
formed by these functions, use is made of the 
fact that if a given permutation transforms one 
spin function into a different function that row 
of the transformation matrix has a zero diagonal 


Tisza, Zeits. f. Physik 82, 48 (1933); E. B. Wilson, Jr., J. 
Chem. Phys. 2, 432 (1934). 
*T, contains planes of symmetry, etc., in addition to the 
pure rotations of T. 
R.S. Mulliken, reference 2. 
~" ry for example, E. B. Wilson, Jr., reference 8, Table 
, 4.4a, 
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TABLE IV. Nuclear spin functions for CH. 
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element; if it transforms a function into itself, 
the diagonal element is unity. The character x’ 
for a given permutation is therefore equal to 
the number of spin functions unchanged by 
the given permutation. Thus the permutation 
(14)(23) (which is equivalent to a twofold rota- 
tion) leaves I and XVI invariant and also two 
functions from VI-XI; the character x’c, is 
therefore 4. Similarly x’c, is found to be 4 
and evidently x’ is 16. Using these data and 
Table I we find that for CH, 


Ys=5A+E+3T. (5) 


For CD,, m,=1, 0, or —1 so that there are 
81 functions ~s, of which typical members are 
given in Table V. From these the characters are 


TABLE V. Nuclear spin functions for CDs. 





—1 
—1 
—1 


eR OCC | 
eee eee 
RAPE RAL 











found to be x’z=81, x’c,=9, x'c,=9 so that 
for CD, 
Ws=15A+6E+187. (6) 


SYMMETRY OF Wrvevvvevs 


The symmetry of the product of two functions 
is obtained by the use of the rule that the 
character x’ for the product of two independent 
functions is equal to the product of the char- 
acters for the two functions. In this way a 
multiplication table for the irreducible repre- 
sentations A, FE, T can be built up. For T it is 
given in Table VI. 

If a function x has the symmetry F,, for 
example, while y has the symmetry 7, the 
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TABLE VI. Multiplication table for irreducible representations 
of T 








Ey E, : 
r 


T 
A+E,+E:+2T 





A 
Ei 
i 


E2 
A 
ig 








product xy has the symmetry given under EF; and 
opposite 7 in Table VI; namely, 7. 

By using the results of Table III, and Eqs. (5) 
or (6) in conjunction with Table VI, the sym- 
metry of the functions YrvevvWevs can be 
obtained. It will be later shown that states in- 
volving nuclear spin functions of different per- 
mutation symmetries do not readily combine. 
It is therefore convenient to tabulate separately 
the states with different spin symmetry. The 
products yrve-vvevs must have the symmetry 
A in order to occur in nature. The number of 
linear combinations of these products (based on 
one framework) with the symmetry A and with 
the spin symmetry and vibration symmetry 
shown are listed in Table VII for CH, and in 
Table VIII for CD,. It will be shown that these 
numbers are the statistical or quantum weights 
of the corresponding energy levels of the mole- 
cule. 


CONSTRUCTION OF THE COMPLETE WAVE 
FUNCTIONS 


The linear combinations enumerated in Tables 
VII and VIII are functions of all the coordinates, 
including spin, and have the correct symmetry 
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properties with respect to those permutations of 
identical nuclei which correspond to rotations of 
the framework of nuclei upon which these func- 
tions are based. They do not have the correct 
symmetry with respect to other permutations of 
the identical atoms, since such permutations 
alter the framework. However, by adding (with 
the proper sign) the correct linear combination 
(call it Ys) for the framework A and the corre- 
sponding linear combination yy, for the frame- 
work B, a function can be constructed which 
has the correct symmetry properties with respect 
to any permutation of identical atoms. Thus the 
function 
V=Ya-—ve (7) 
will be completely antisymmetric for any inter- 
change of protons; for example 
(12)¥W=yYs-—Ypa=—-YV, (8) 


while 


(123) W = (123) Pa — (123) Ye 
=Ya—vs=V,etc. (9) 
Furthermore no new degeneracy is introduced 
by this step, because even when there are many 
independent frameworks for a molecule there 
will be only one combination of the functions 
based on the different frameworks with the 
correct symmetry. 


To prove this, let us consider the representation of the 
entire permutation group (not the subgroup we have pre- 
viously used) by the set of functions Ya, Wz, etc., one for 
each framework, the rotational properties of each being the 


TABLE VII. Weights of rotation-vibration levels for CH. 
(p=0, 1, 2, ---. Multiply each weight by 2/+1.) 








YE¥y Ws J=6p 6p+1 
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3(3p+2) 
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10p 10p 

4p+2 4p 
18p 18p+6 
32p+2 32p+6 


10(p+1) 
4p+2 
18p+6 

32p+18 


10(p+1) 
4p+2 

18p+18 

32p+30 


10p 10(p+1) 
4p+2 
18p+12 


32p+14 32p+26 





15p 15p+5 
6p 6p+2 
27p+3 


27p+6 
48p+3 48p+13 


15p+5 
6p+2 

27p+12 

48p+19 


15p+15 
6p+6 
27p+ 24 
48p+45 


15p+10 
6p+4 

27p+21 

48p435 


15p+10 
6p+4 

27p+15 

48p+29 
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TABLE VIII. Weights of rotation-vibration levels for CDs. 
(p=0, 1, 2, ---. Multiply each weight by 2/+1.) 








vey vs J=6p 6p+1 


6p+2 


6p+3 6p+4 6p+5 





15A 15p+15 

‘ 6E 12p 
“) 18T 54p 54p+18 
| Total 81p+15 81p+18 


15p 
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12p+12 
54p+18 
81p+30 


15p+15 
12p+12 
54p+36 
81p+63 


15p 
12p+12 
54p+54 
81p+66 


15p+15 
12 


54p+36 
81p+Si 





6E 24p+12 24p 
18T 108p 108p+36 
| Total 162p+12 162p+36 


| 15A 30p 30p 
E 


30p+30 30p 
24p+12 
108p+36 
162p+78 


30p+30 

24p+24 
108p+72 
162p+126 


30p +30 

24p+12 
108p +108 
162p+150 


24p+12 
108p+72 
162p+84 





15A 45p 45p+15 
r 6E 36p 36p+12 
18T 162p+18 162p+36 
Total 243p+18 243p+63 


45p+15 

36p+12 
162p+72 
243p+99 


45p+30 

36p+24 
162p+90 
243p+144 


45p+30 

36p+24 
162p +126 
243p+180 


45p+45 

36p+36 
162p+144 
243p+225 








same. There will be N/o of these functions if N is the 
number of permutations in the group and o the number of 
permutations equivalent to rotations (the symmetry 
number). Assume first that we seek symmetric W’s. For 
each Wa, etc., there will be « permutations transforming Wa 
into itself. Each function will therefore contribute the 
quantity o to 2(j)x’; and since there are N/o functions, 
2(j)x';=N. The number of linear combinations in the 
symmetrical irreducible representation is m= (1/N)=(j)x’; 
since every permutation has the character unity in this 
irreducible representation. We thus have m,=1 as stated. 
For the antisymmetrical case, the same result is true, 
because whereas odd permutations have the character 
minus one in the antisymmetric irreducible representation, 
these same permutations will either change one function 
into another, contributing zero to the character, or they 
will transform ya, etc., into —wWa contributing —1 for 
each function so transformed. The minus signs cancel and 
we have the result as before that there is only one com- 
pletely antisymmetric combination. 


The numbers in Tables VII and VIII, which are 
the degeneracies of 4 or Wz are therefore the 
correct degeneracies or weights of WV. 

As has been pointed out by other writers,” 
statistical mechanical calculations for high tem- 
peratures, at which many rotational levels are 
occupied, can be made with sufficient accuracy 
by disregarding the alternation in the weights of 
the various rotational states, assigning to each 
state the average weight obtained by dividing 
the total nuclear spin multiplicity (16 for CH,, 
81 for CD,4) by the classical symmetry number 
(here 12). That this is the correct average 
weight can be checked by the use of the numbers 
in Tables VII and VIII. 

?Ludloff, Zeits. f. Physik 57, 227 (1929); Mayer, 


Brunauer and Mayer, J. Am. Chem. Soc. 55, 37 (1933); 
Kassel, J. Am. Chem. Soc. 55, 1351 (1933). 


NONCOMBINING SPECIES 


Most types of perturbation which cause transi- 
tions do not involve the nuclear spin coordinates. 
Consequently, two states of a molecule with 
mutually orthogonal spin functions will be non- 
combining if only such perturbations occur. A 
well-known example is the case of the hydrogen 
molecule which may have either singlet or 
triplet spin functions and therefore exists in two 
noncombining forms (in the absence of catalysts). 
The deuterium molecule De, on the other hand, 
can have singlet, triplet, and quintet nuclear 
spin states, which presumably have long life- 
times. However, the singlet and quintet states 
have the same permutation symmetry and there- 
fore are associated with the same set of rotational 
levels. Their equilibrium proportions will thus 
be the same at all attainable temperatures and 
there will be no way of changing their relative 
proportions. 

A similar situation exists for polyatomic mole- 
cules. The spin functions can be divided accord- 
ing to their permutation symmetry and according 
to their multiplicity. States with spin functions 
differing either in symmetry or multiplicity will 
rarely combine unless a catalyst is present. 
However, as in De, the equilibrium proportions 
of certain sets of noncombining states will be 
unaffected by temperature so that for all prac- 
tical purposes these sets of states constitute the 
noncombining states of interest. If two spin 
functions belong to the same irreducible repre- 
sentation of the group of permutations equivalent 
to rotations of the molecule (a subgroup of the 
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complete permutation group) the corresponding 
species will always exist in the same proportions 
and may be considered for practical purposes as 
one species. If they belong to different irreducible 
representations, they may have equilibrium pro- 
portions varying with the temperature. This is 
the justification for the statement previously 
made that there are three noncombining species 
in CH, and CD,, having nuclear spin functions 
with symmetries A, E and 7, respectively. 


SYMMETRY OF ROTATIONAL FUNCTIONS FOR SYM- 
METRICAL Top MOLECULES 


In discussing other molecules it is necessary to 
know the transformation properties of sym- 
metrical top functions. There are two types of 
rotations of interest. One is a rotation by an 
angle 27j/n about the symmetry axis, where j 
and » are integers. Such a rotation changes the 
Eulerian angle x into x+27j/n and therefore 
has the effect 


Cri 
Yre=WsxulO, 9, x) vsKmM(O, ¢, x+27j/n) 


=e? Kind yeu (8, 9, x); 


(10) 


since the angle x enters Ye only through the 
factor* e'**, Since there are only two functions 
for a given value of J and M and for a given 
energy; namely, ~yxw and ~y_xm, the character 
for such a rotation C is 


x’ co =erriKin4 e-2riiKin=) cos (2rjK/n), (11) 
if K¥0, and unity if K=0. The other type of 
rotation is a twofold rotation about an axis 
perpendicular to the symmetry axis and making 
an angle a with the axis X’ which is the origin 
of the angle x measuring the rotation of the 
molecule about its symmetry axis. This rotation 
corresponds to changing @ into r—0, gy into r+¢ 
and as can be seen from Fig. 3, x into r— x+2a. 
Using the properties‘ of the symmetrical top 
functions we can find that 


Voxu(r—8, +e, m™—x+2a) 


= (—1) Ith MeiK(rt20)) 5 _ ey (0, 9, x), 


(12) 


where ¢ is the larger of |K| and | M|. Therefore 
if K+0, the character for this rotation is zero; 
if K=0, the character is (—1)/. 


x A 


» N fy N 


Se 
—— 


A A 








Fic. 3. The X’ Y’ plane of a symmetrical top, showing the 
effect on the angle x of a twofold rotation about an axis AA 
in the X’Y’ plane. NN is the nodal line. x’, the angle after 
the rotation, equals —x+2a. 


REsULTs FoR CH;X, CD;X, NH; Ano ND, 


These molecules belong to the rotational point 
group C;, which has the elements E, C3 and C;?. 
Symmetrical top functions have the characters 
given in Table IX for these symmetry elements, 
obtained by using Eq. (11). Using these and the 
characters x,‘ for C3 we find that the rotational 
functions have the symmetries shown in Table X. 

The spin functions for CH;X (if X has no 
spin) are given in Table XI. From this it is 
seen that x’n2=8, x’c,=2, x’c,=2, so that ys 
has the symmetry 44 +2E. Combination of this 
result with the data of Table X yields Table XII 


TABLE IX. Characters for the representation of C; formed ° 
by WR. 
(p=0, 1, 2, ---. Multiply all characters by 2J+1.) 











K=0 
|K| =3p 
|K| =3p+1 








TABLE X. Symmetry of the rotational functions for a molecule 
with rotational symmetry C3. 


(p=0, 1, 2, ---.) 





K=0 
|K| =3p 
|K| =3p+1 


(2J+1)A 
2(2J+1)A 
(2IJ+1)E 


TABLE XI. Nuclear spin functions for CH;X. 
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13 See Mulliken, reference 2, for values of x; for C.. 
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TABLE XII. Weights of the vibration-rotation levels of CH;X. 
(p=0, 1, 2, ---. Multiply all weights by 2J+1. 
No spin on X.) 











A ——_—— -k——— 
4A 2E Total 4A 2E Total 





4 0 4 4 

8 0 8 

12 

when the multiplication properties of C,; are 

used, which are: AXA=A, AXE=E, EXE 

=2A+E. The complete wave function must 
have the symmetry A. 

If X has a spin s, these weights must be 
multiplied by 2s+1. Thus if X is D, a deuterium 
nucleus, 2s+1 has the value 3. Table XII also 
can be used for NHsg, if all weights are multiplied 
by three, since nitrogen has a unit nuclear spin. 

A similar treatment for CD3;X shows that if X 
has no spin Ys has the symmetry 114+8E, 
from which the weights in Table XIII may be 
deduced. If X has a spin s, these weights must 
be multiplied by 2s+1, so that the factor for 
CD;H is 2. By multiplying all weights by 3, 
Table XIII can also be used for ND3. 


TABLE XIII. Weights for the vibration-rotation levels of CD;X. 
, 2, +++. Multiply all weights by 2J+1. 
X has no spin.) 
en Sees 
114 8E Total 


Total 





K =3p+1 








THE SPLITTING OF ENERGY LEVELS 


In certain cases the energy levels which we 
have previously considered as degenerate are 
really split into two or more experimentally 
separable components. A well-known example is 
NH; in which the energy levels are nearly all 
doublets.* This splitting is due to the fact that 
the nitrogen atom can pass through the plane of 
the three hydrogen atoms without great difficulty 
so that it has two equivalent positions of 
minimum potential energy. Such an inversion 
corresponds to an interchange of identical hydro- 
gen atoms (which changes the sign of W), 


_—_ 


's Dennison and Hardy, Phys. Rev. 39, 938 (1932). 
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followed by a rotation of the molecule, which 
may or may not again change the sign of W, 
depending on the nature of the functions yr. The 
rotation required is equivalent to the trans- 
formation 6-7-8, g-r+¢, x-2r-—x, which 
has the following effect on the symmetric top 
function YyKwv: 


Vakul, ¢, x) WsKku(r—8, r+ ¢, r—x) 
=(—1)74 AMY 7 (0, 9, x), 


where ¢ is the larger of K)| and J/'. From this 
one finds that the levels with K=0 are not 
split but have total wave functions W which are 
antisymmetric with respect to inversion when J 
is even and symmetric when J is odd. When 
K+#0, the energy level is split into two com- 
ponents to each of which corresponds an equal 
number of wave functions, one component having 
functions which are symmetric towards inversion, 
the other component having antisymmetric func- 
tions. 

The criterion as to whether or not the splitting 
is appreciable is the magnitude of the potential 
hill separating the two forms. Methane also 
possesses two forms but to pass from one to the 
other would require passing through too great a 
potential barrier so that no doubling is ob- 
served. 

Since the electric moment changes sign on 
inversion, transitions in absorption occur only 
between states of opposite symmetry character 
with respect to inversion, whereas Raman transi- 
tions occur between states of like symmetry. 


WEIGHTS FOR BENZENE 


If benzene is a regular plane hexagon it has 
the point group symmetry Dy,,, the rotational 
subgroup for which is Dy, with the elements E£, 
Co, 2C3, 2Cs, 3C2’, 3C2’’. Benzene is a sym- 
metrical top, and so by using Eqs. (11) and 
(12) we can obtain the characters x’ shown in 
Table XIV for the representation of D, formed 
by Wr. Using the character table for Ds given by 
Mulliken? and the values in Table XIV, we 
obtain the symmetries of the linear combinations 
of wr given in Table XV. 

The 2° nuclear spin functions for CH, are 
easily shown by the methods employed for the 
other molecules to combine to functions with the 
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TABLE XIV. Characters for the representation of Ds formed by vr. 


(p, k=0, 1, 2, 3, 


++, Multiply all characters by 2J/+1.) 





, 


X ¢ 


if 


X 4 








K=0 . ; 
K =6p+k 


2 cos 2rk/3 


1 
(—1)*2 cos 27k /3 





symmetries 13A,+A.+7B\+3B.+11F*+9EFx*. 
The total wave function for the molecule must 
fall in class B,; since this represents complete 
antisymmetry. 


TABLE XV. Symmetry of the rotational functions for a mole- 
cule with rotational symmetry Dg. 


(p=0, 1, 2, 3, ---.) 
K=0, J even: 
K=0, J odd: 
K=6p: 
K=6p+1: 
K=6p+2: 
K=6p+3: 





(2J+1)A, 
(2J+1)A2 
(2J+1)(A,+A2) 
(2J+1)E,* 
(2J+1)E* 
(2J+1)(Bi+Bz2) 


The multiplication properties of the repre- 
sentations of D, are 


AzXA:=A;, AcXB;=B;, A2XE*=E*, AeXE,*=E-*, 

BiXBs=A,, BsXBe=As, BsXE*=E,*, Bs XE,*=E"*, 

E*X E* = Ex* X Ey* =A, +A2.+E*, E* XE,* =Bi+Be 
+E,*; where i=1, 2, 7 =2, 1. 


With the use of these we can find the symmetry 
of Ye’vWrwvs and from that the number of com- 
binations with the total symmetry B;; i.e., the 
quantum weights. These are given in Table XVI. 
Only the total weights are given. 


TABLE XVI. Weights of the vibration-rotation levels of CeH,. 
(p=0, 1, 2,3, --+. Multiply all weights by 2/+1.) 








wevr: y Ao 


K=0, J even 
K=0, J odd 
K=6p 

K =6p+1 
K=6p+2 
K=6p+3 











CYCLOPROPANE 


The most reasonable model for C3;H¢, cyclo- 
propane, has the point group symmetry Dz, 
for which the rotational subgroup is D3. Since 
the work of finding the weights for this molecule 
follows exactly the method previously used, only 


the results will be listed. These are given in 

Table XVII. 

TABLE XVII. Weights of the vibration-rotation levels of C;H,. 
(p=0, 1, 2, 3, ---. Multiply all weights by 2J/+1.) 








vevyv: A, 





& 
odd 16 
24 





ETHYLENE—AN ASYMMETRIC Top 


Ethylene!’ will serve as an illustration of how 
this method is to be used for asymmetric top 
molecules. The rotational functions. for the 
asymmetrical top are of four kinds, differing in 
their symmetry with respect to the twofold 
rotations about the three principal axes of the 
molecule. Table XVIII gives the symmetries of 
TaBLeE XVIII. Characters for the representation of V formed 


by the four kinds of asymmetric top functions. 








| Irred. 
| Rep. 











these four kinds, which are labeled by the nota- 
tion ++, + —, etc., introduced by Dennison.” 

The spin functions for CsH, are 16 in number 
and have the symmetry 7A,;+3B,+3B2+3B;, 
as can be verified by using the methods pre- 
viously described. The multiplication table fer 
V is: AiXA\=A,, A,XB;=B,;, B;X B;=A\, B, 
< B.=B;3, BiXB3;=B2, B,_XB3=B;. The sym- 
metry of veWvrws can be found by the use of 
this information and the number of combinations 
with the symmetry A, (which corresponds to 


antisymmetric or symmetric wave functions) 


1s These weights have been previously calculated by 
Professor R. M. Badger but have not been published. 
16 PD), M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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TABLE XIX. Weights for the vibration-rotation levels of CIT. 
(Multiply all weights by 2/+1.) 





CONSTANT 


STUDIES. 


TABLE XX. Weights for the vibration-rotation levels of CH2D2. 
(Multiply all weights by 2/+1.) 

















obtained. These are the weights given in Table 
XIX. 

The molecule CHD». (methane) is also an 
asymmetric top with the rotational symmetry 
C,. There will be only two symmetry classes of 
rotational levels (with respect to Ce), those 


sym. (K even) 
antisym. (K odd) 








symmetrical to the twofold axis and those anti- 
symmetrical. The vibrational and spin functions 
can likewise be divided into these two classes. 
The weights obtained are given in Table XX. 

In conclusion I should like to express my 
appreciation to Professor J. H. Van Vleck for 
suggesting this problem and for valuable criticism 
in connection with it. 
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Dielectric Constant Studies 


V. Anomalous Dispersion of Lecithin in Viscous Mineral Oils 


A. L. FerGuson, L. O. Case And G. HARLOWE Evans,* Chemical Laboratory, University of Michigan 


(Received January 14, 1935) 


A unique receiver and accompanying cell for determining 
dielectric constants by the resonance method were con- 
structed. Their design is such as to reduce to a minimum 
several types of error inherent in the resonance method. 
Dielectric constant measurements were made on five solu- 
tions of lecithin in very viscous mineral oils at three tem- 
peratures and at six wavelengths. Three distinct methods 
were employed: (1) refractometer method—frequency of 
D line of sodium, (2) resonance method—32.8 m, 94 m, 
320 m, 1308 m, (3) bridge method—precision Wheatstone 


VIDENCE is becoming more convincing 

that the viscosity affecting the rotation of 
dipole molecules in an alternating field, referred 
to as microscopic viscosity, is not the same as 
that measured by the ordinary method, known as 
macroscopic viscosity. 

The determination of microscopic viscosity by 
means of dielectric constant measurements is 
based upon the theoretical work of Debye on 
anomalous dispersion, the decrease in dielectric 
constant with increasing frequency. A large ab- 
aa 


* Rewritten from a thesis by G. Harlowe Evans in 
partial fulfilment of the requirements for the degree of 
Doctor of Philosophy at the University of Michigan. 


bridge, 1036 cycles. Analyses of the experimental curves 
and comparisons with theoretical Debye curves definitely 
indicate dispersion though not of a simple type. Another 
series of measurements made on one solution at 16 different 
wavelengths covering the range from 40 m to 5,000,000 m 
still failed to show any individual simple dispersion regions. 
A qualitative explanation based on varying viscosities of 
the components of the solvent, the presence of polar radicals 
in the solute molecules, and the extensive association of the 
solute is offered. 


sorption of energy occurs simultaneously with 


anomalous dispersion. 

In an alternating electric field four things in- 
creasingly oppose the turning of molecular di- 
poles. These are: Increasing viscosity, increasing 
particle size, increasing frequency, and decreas- 
ing temperature, the last due largely to its effect 
on the viscosity. 

The effect of these factors is to throw the 
orientation of the molecules out of phase with the 
applied field. It is this lagging of the molecules 
behind the field that causes the great absorption 
of energy in the region of anomalous disper- 
sion. 
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Fic. 1. Cell and receiver assembly. 


It thus appears that by measurements of dielec- 
tric constants in the region of anomalous disper- 
sion one should be able to measure the viscosity 
which is really involved in the turning of the 
dipole molecules. 

Measurements of this kind are reported in 
the literature for extremely large molecules of 
uncertain molecular weight in fairly fluid sol- 
vents!’ ?»* and also for small well-defined mole- 
cules in very viscous solvents.*: > However, no 
information seems to be available on systems of 
molecules of large but definite molecular weight in 
solvents of moderately high viscosities. It was 
the original purpose of this investigation to study 


the relation between microscopic and macro- 


scopic viscosity in such a system. 

Solutions of lecithin (molecular weight 804) in 
Sun Oils S. A. E. 30, 40, 50, and 60 were found to 
possess the required characteristics. These oils 


1J. Wyman, Jr., J. Biol. Chem. 90, 443 (1931). 

2N. Marinesco, J. chim. phys. 27, 454 (1930); 28, 51 
(1931); Comptes rendus 192, 625 (1931). 

3 J. Errera, J. chim. phys. 29, 577 (1932). 

4 Goldammer-Sack, Physik. Zeits. 31, 224 (1930). 

5 R. Luthi, Helv. Phys. Acta 6, 139 (1933). 
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Fic. 2. The dielectric constant cell No. 3. 


are unusually pure and cover the viscosity 
range from 3.66 to 21.5 poises at 25°C. For Sun 
Oils 30, 40 and 60, 10 percent solutions were used, 
while for Sun Oil 50 two solutions containing 3 
percent and 10 percent were employed. 

With these solutions dielectric constants were 
measured at six frequencies and at three tem- 
peratures. The results allow a study of the influ- 
ence of viscosity, temperature and concentration 
on the dielectric constant. 

The values of e) (the dielectric constant at 
optical frequencies) were obtained from refrac- 
tive indices measured by means of an Abbe 
refractometer, making use of the Maxwellian 
relation ?=e where n is the refractive index. 
Assuming that 300,000 meters (1000 cycles) is a 
long enough wavelength, values of «,, (the so- 
called “‘static’’ dielectric constant) were deter- 
mined at this wavelength by means of a precision 
Wheatstone bridge assembly developed by Mr. 
R. D. Thompson in this laboratory. 

The values of e’ (the measured dielectric con- 
stant) at wavelengths of 32.8 m, 94 m, 320 m, 
1308 m were calculated from capacitance meas- 
urements made by the substitution method em- 
ploying a simple “‘voltage tuning’”’ resonance ap- 
paratus very similar in electrical arrangement to 
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one used previously by Malone*? in this labora- 
tory. Malone found this to be the most suitable 
method of determining the dielectric constants 
of small amounts of liquids over a wide range of 
wavelengths. 

For the voltage-tuning resonance method at 
high frequencies, Lattey* has summarized the 
sources of error which should be eliminated or for 
which corrections should be made. These are: 
(1) The conductance of the dielectric or leakage, 
(2) the earth capacitance, (3) the resistance of 
the leads, (4) the capacitance of the leads and 
(5) the inductance of the leads. 

In the present work these sources of error were 
greatly reduced by the design of the apparatus, 
Fig. 1. The earth capacitance was very largely 
removed by the design of the measuring cell. 
Except for the small Bakelite insulator, the cell, 
Fig. 2, was constructed entirely of metal and 
provided with a water jacket for temperature con- 
trol. During measurements the outer plate is 
connected to the ground and only the electric 
field through the insulator can cause an earth 
capacitance. 

The errors due to resistance were practically 
eliminated by the. use of heavy brass strips for 
all leads in the resonance circuit. The errors due 
to capacitance and inductance were reduced to a 
minimum by making the leads as short as pos- 
sible, especially those connecting the cell and 
the precision condenser. Reference to Fig. 1 will 
indicate how these objectives have been at- 
tained. 

By the substitution method used, the capaci- 
tance C,* is obtained as the difference in readings 
on the precision condenser at the resonance 
points when the measuring cell is in the circuit 
and when it is out. The differences are obtained 
by a method similar to that described by 
Malone.® 

It was found that the apparent capacitance of 
the empty cell changes with waveleng'th, being 
noticeably larger at longer wavelengths. One 
might expect that the capacitance of the cell full 
of liquid would increase in approximately the 

* Malone-Ferguson-Case, J. Chem. Phys. 1, 836 (1933), 

“R. Gunn, Phil. Mag. 48, 224 (1924). 

* Lattey-Gatty, Phil. Mag. 7, 985 (1929). 

* Since the relation between capacitance and condenser 
readings is linear over the range employed, differences in 


readings were used in all calculations rather than actual 
Capacitances. 
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same ratio as that of the empty cell so that the 
relation e’=C,/C, would be maintained. How- 
ever, the dielectric constants calculated thus are 
too small indicating that the measured capaci- 
tance of the empty cell is too large. This is 
probably due to the unavoidable presence at 
some point, in any cell which might be used, of 
an insulating material having a greater dielectric 
constant than air, its relative effect being greater 
when the condenser is empty. 

A test using benzene and carbon disulphide, 
however, showed that the relation 


Cc/Cpr= €c’/en’ =a constant 


was independent of frequency. Here Cc and Cz 
are differences in condenser readings for carbon 
disulphide and benzene, non-polar liquids of 
known dielectric constants ec’ and e€,’, respec- 
tively. Since the dielectric constants measured in 
this work were intermediate between those of 
benzene and CSz, it was assumed that the relation 
would apply in all such measurements. As a 
generalized expression e,’=(Cx/Cr)en’ was used 
for the unknown material. A further confirma- 
tion of this relation was later furnished by the 
fact that the dielectric constants of Sun Oil 30 
calculated therefrom were constant throughout 
the wavelength range employed. 

The values for ¢’ at the four intermediate wave- 
lengths are plotted to scale against wavelength 
in Figs. 3-7. The values at the two extreme wave- 
lengths spoken of earlier are also shown but are 
not plotted to scale horizontally. That these 
curves represent portions of true but not simple 
dispersion regions can be best shown by an 
analysis of a set of theoretical dispersion curves, 
Fig. 8. For curve B the assumed values were 
€o= 2.46, €,=2.60 and 7=(4/r)xX10-*. Other 
necessary data appear on the figures. The effect 
of temperature, largely due to resulting changes 
in viscosity, is shown by a comparison of curves 
A, Band C. The density change with temperature 
is insignificant, but the change in viscosity 1; is 
very great if it is assumed that the temperature 
effect is roughly the same on microscopic as it is 
on macroscopic viscosity. The viscosity ratios 
nso° : N2s° : Nige= 0.1 : 1 : 2 used for the curves are 
approximately those for the macroscopic viscosi- 
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percent lecithin. 


LENG IN TER 


Fic. 4. Dispersion curves. Sun Oil 40 containing 10 
percent lecithin. 
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Fic. 5. Dispersion curves. Sun Oil 50 containing 3 
percent lecithin. 


ties of Sun Oil 30. The values of ¢)9 chosen were 
2.47 at 18° and 2.42 at 50°. In obtaining these 
values it was assumed that the effect of tem- 
perature on €9 was the same as that measured for 
n* for the solutions. 

Curves D, B and E show the effect of change 
in particle size or in viscosity at constant tem- 
perature. As seen in the figure the values of a? 
for these three curves are in the ratios of 1 : 4 : 16. 
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Fic. 6. Dispersion curves. Sun Oil 50 containing 10 
percent lecithin. 
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Fic. 7. Dispersion curves. Sun Oil 60 containing 10 
percent lecithin. 
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Fic. 8. Theoretical dispersion curves. 


It will be seen that in the region of low wave- 
lengths in all cases the values of e’ are small and 
approach ¢. In the high wavelength region the 
values of e’ are large and very rapidly become 
constant. In each of these regions the tempera- 
ture coefficient of the dielectric constant is nega- 
tive (curves A, B and C). However, in the region 
of dispersion the temperature coefficient is de- 
cidedly positive. Thus the general effect of in- 
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creasing the temperature is to shift the whole 
curve to somewhat lower values of dielectric con- 
stants and to considerably lower wavelengths. 
For instance, on curve A, by increasing tempera- 
ture, the points (a) and (b) are shifted to the 
corresponding positions indicated on curves B 
and C. The effect of increasing either the size a* 
of the particles or the viscosity 7; of the solvent 
is to shift the whole curve to higher wavelengths, 
as shown by curves B, D and E. In the light of 
these theoretical curves it can be shown that 
Figs. 3-7 represent the behavior of materials 
which are dispersing. 

On examination of Figs. 3, 4, 6 and 7 it is seen 
that the temperature coefficient is always nega- 
tive for the values of «9 =n”, while for the value 
e’ at wavelengths from 32.8 m to 1308 m it is 
positive in all cases. Thus there is a complete 
reversal of the magnitude of e’ with tempera- 
ture. This is characteristic of the lower end of 
dispersion curves as will be seen in Fig. 8. Again, 
it is observed that the values of e’ for the two less 
viscous solutions (S. A. E. 30 and 40) at 1036 
cycles (the extreme upper value) and 50° have 
dropped below the corresponding values at 25° 
and that while the corresponding value for the 
S. A. E. 50 solution has approached the value at 
25°, it has not yet dropped below it. The corre- 
sponding value for S. A. E. 60 solution has main- 
tained its relative position. 

These reversals for the less viscous solutions 
would indicate that the values at 1036 cycles are 
either near or upon the upper plateau of the dis- 
persion curves. From an examination of the 
theoretical curves, Fig. 8, it would be expected 
that the less viscous solutions at the highest 
temperatures would reach the top of the curve at 
lower wavelengths while the more viscous solu- 
tions would approach it at higher wavelengths. 
There is no indication in the range covered of the 
value at 25° on any of the curves tending to drop 
below that for 18°. 

The apparent discrepancy in the height of cor- 
responding curves for the different oils is due to 
the difference in the dielectric constants of the 
pure oils. When corrections for the latter are 
made the curves for the different solutions at 
the same temperature assume their expected 
relative positions. Furthermore, the magnitude of 
these shifts is of the right order since the change 
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in viscosity from one oil to the next is less than 
the change in viscosity of any one oil due to 
temperature in passing from 18° to 25°. 

In Fig. 5 for the 3 percent solution, dispersion 
is shown even more definitely by the fact that the 
order is completely reversed not only between 
optical frequencies and 32.8 m but also at the 
other end of the curves. In fact, the curve for 
50° begins to approach that for 25° between 200 
and 300 m, crossing it at about 1000 m and fall- 
ing much below even those for 18° at 1036 cycles 
(289,570 m). 

The fact that the value of e’ at 25° has dropped 
below that for 18° at 1036 cycles indicates, also, 
that for this concentration dispersion occurs at 
lower wavelengths than for the higher con- 
centration. 

Since the temperatures are the same and the 
viscosities greater than those of the Sun Oil 30 
and 40 solutions, the lowered dispersion range 
must be ascribed, in accordance with the Debye 
theory, at least, to the presence of smaller par- 
ticles. This clearly indicates association of the 
lecithin molecules. These phenomena all indicate 
distinct dispersion, but it must be admitted 
that considering the rapid dispersion at low 
wavelengths the dispersion covers too great a 
range to be a simple dispersion region. 

It becomes evident that it is desirable to obtain 
a larger number of points on the curve for one 
solution and to extend the measurements to 
longer wavelengths in order to determine the 
possibility of the presence of irregularities in the 
curve indicating the presence of simple overlap- 
ping dispersion curves. 
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Fic. 9. Dispersion curve. Sun Oil 3 containing 10 percent 
lecithin, at 50°C, 
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The system used for this extended work was a 
solution containing 10 percent lecithin in Sun Oil 
30 at 50°C. A set of measurements was made 
throughout the whole wavelength range on 
benzene and also on pure Sun Oil 30. The value 
at all frequencies for Sun Oil 30 was 2.293+0.003. 

The curve of Fig. 9 is a plot of the data ob- 
tained.* Curve B represents the entire range of 
wavelengths down to 6450 m while curve A shows 
the lower portion from 6450 m to zero meters 
expanded horizontally 100-fold. 

Even with these more extensive and precise 
data no evidence is obtained of separate disper- 
sion regions. As a matter of fact, except for the 
two extreme points at zero and infinite wave- 
lengths the dielectric constant is found to be, 
within experimental error, a linear function of the 
logarithm of the wavelength and is represented 
by 

e’ = 2.324+0.0619 log X. 


When the complex nature of the systems is 
considered the absence of the appearance of 
single dispersion regions is not surprising. Any 
one of the oils is doubtless a complex mixture of 
molecules of various sizes. The viscosity of each 
individual component is probably different and 
thus the relaxation time of a given polar mole- 
cule caused by each would be different. The re- 
sulting dispersion curve would therefore be a 
summation of several simple dispersion regions. 

It has been suggested by Race" and by Luthi® 
that polar radicals within the molecule as well as 
the molecule as a whole might cause dispersion. 
Lecithin is supposed to consist of several long 
chain compounds joined together in one large 
molecule. The molecule as a whole would cease 
to orient at very long wavelengths, and the 
various radicals would cease with frequency in 


* The value at 39.2 m is obviously in error, but circum- 
stances made it impossible to recheck it after plotting the 


curve. 
10H. H. Race, J. Phys. Chem. 36, 1928 (1932). 


CASE AND G. i. EVANS 
the order of their size, each having a different 
relaxation time. This likewise would result in a 
complex dispersion curve. Either or both of these 
explanations could account for the lower portion 
of the dispersion curves obtained. 

The upper portions of the complex dispersion 
curves may be caused by the presence of asso- 
ciated groups of different sizes. It was shown in 
the discussion of Fig. 5 that lecithin is highly as- 
sociated. These associated particles might readily 
approach colloidal dimensions and thus possess 
long relaxation periods.* 

Recently several other investigators have ob- 
tained dispersion curves which failed to show 
discrete dispersion regions. Among these may be 
mentioned Kitchin,'! Race and Williams” each 
of whom experienced difficulty in the interpreta- 
tion of the data obtained in terms of the Debye 
theory. In some cases the difficulty doubtless 
would have been enhanced if a wider frequency 
range had been covered. Williams states ‘In 
conclusion, it may be said that these results show 
that the Debye theory, although theoretically 
sound, is very difficult to apply quantitatively to 
the existing dielectric constant data. . . . At the 
present time it does not seem readily possible to 
obtain monodisperse systems either of large 
molecules or of particles of colloidal size, in non- 
polar solvents.”’ 

Among the unique features of the present in- 
vestigation are: (1) The design and assembly of 
the cell and receiver to minimize various errors 
present in such measurements, (2) the large 
number and extremely wide range of wavelengths 
employed, (3) the type of solute-solvent system 
employed, and (4) the relatively large change of 
dielectric constant with frequency for such 
viscous solutions. 


* An. approximate calculation shows that associated 
particles of about 1000 molecules are sufficient to account 
for the dispersion at the upper end of the curve. 

 —. W. Kitchin, J. Ind. and Eng. Chem. 24, 549 (1932). 

12 Oncley-Williams, Phys. Rev. 43, 341 (1933). 
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The Kinetics of the Decomposition of Gaseous Glyoxal 
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The decomposition of gaseous glyoxal has been investigated from 410° to 450°C. The reaction 
is homogeneous and first order. Considerable carbon and tar is deposited, however, and the 
results are erratic. The following three reactions occur to about the same extent: 


CHO —CHO =2 CO+Hb, 
CHO 


1 
CHO CHO= 
n 


CHO =C+CO+H,0, 


(CHO—CHO),. 


The results are discussed from the point of view of the occurrence of glyoxal as an inter- 


mediate in the oxidation of acetylene. 





LYOXAL has frequently been postulated as 
a product of the oxidation of gaseous or- 
ganic compounds. It has been definitely shown 
to be present in considerable quantities during 
the oxidation of acetylene.' It has been suggested 
that the primary step in the acetylene oxidation 
is the formation of glyoxal, which then either 
oxidizes further or decomposes. Although, as far 
as We are aware, no investigation of the thermal 
decomposition has ever been made, it has always 
been assumed that it occurs as follows? 


CHO—CHO=CO+HCHO, 
HCHO=CO+ Hz, 


giving the overall reaction 
CHO—CHO=2 CO+H:. 


The present investigation was therefore under- 
taken (a) because of the interest attaching to the 
decomposition of simple molecules, (b) to fill an 
experimental gap in the Bone “hydroxylation” 
scheme for the oxidation of gaseous organic 
compounds, and (c) in connection with an in- 
vestigation of the acetylene oxidation now in 
progress in this laboratory. 


EXPERIMENTAL 


Rates of reaction were followed by admitting 
the reactants to a heated Pyrex bulb, and ob- 


' Kistiakowsky and Lenher, J. Am. Chem. Soc. 52, 3785 
(1930); Spence and Kistiakowsky, ibid. 52, 4837 (1930); 
Lenher, ibid. 53, 2962 (1931). 

*See, for example, the reference 1 and also Bone, Proc. 
~4a Soc. A137, 243 (1932); J. Roy. Soc. Arts 81,550 
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serving the rate of pressure increase at constant 
volume by means of a capillary manometer. The 
apparatus was similar to one which has been 
previously described.* The reaction vessel was 
heated by means of an electric furnace, the tem- 
perature being controlled by the hand regulation 
of rheostats. The temperature could be main- 
tained constant to within 0.5°C. Temperatures 
were measured with a chromel-alumel thermo- 
couple in conjunction with a thermocouple poten- 
tiometer. The connecting tubing between the re- 
action vessel and the manometer was maintained 
at 100°C by means of an electrical heating ele- 
ment to prevent condensation. 

The preparation of pure gaseous glyoxal is 
complicated by the fact that it polymerizes 
readily to solid paraglyoxal, and the polymeriza- 
tion is not directly reversible. The gas was 
prepared in the usual way by heating paraglyoxal 
(Schuchardt) with phosphorus pentoxide.‘ Con- 
siderable dehydration and charring occurs, how- 
ever, and the gas is contaminated with a large 
amount of carbon monoxide. It was found that if 
the gas were condensed at —80° it froze out as 
the monomer, and it was possible to pump off 
foreign gases. If the solid monomer were heated 
to room temperature for a short time it evapo- 
rated with very little polymerization. During the 
course of the investigation the glyoxal was 
therefore stored as a solid at — 80°, and was only 
warmed up when a sample of gas was to be with- 
drawn. Blank experiments showed that no com- 


3 Hatcher, Steacie and Howland, Can. J. Research 7, 149 
(1932). 
4 Harries and Temme, Ber. 40, 165 (1907). 










































plication was introduced by polymerization in 
the connecting tubes. 

It should be mentioned at the outset that the 
reaction is complicated by the deposition of car- 
bon, tar, etc., and that in consequence the results 
obtained are rather erratic. It is therefore im- 
possible to investigate the reaction with very 
high accuracy, but sufficient results have been 
obtained to establish the main features of the 
process. 


THE PRODUCTS OF THE REACTION 


If the reaction proceeds according to the 
mechanism given above, we should have a pres- 
sure increase of 200 percent at completion, since 
formaldehyde would decompose at the tempera- 
tures used. Actually, however, the average pres- 
sure increase at completion was about 90 percent. 
Some typical values are given in Table I. 


TABLE I. The pressure increase accompanying the reaction. 











Temp. 
(°C) Percent pressure increases 
450 93.0, 83.8, 85.3, 103.0, 94.5, 86.6 
440 87.4, 93.2, 89.4, 92.5, 93.2, 94.0, 101.2, 94.7 
430 73.5, 71.1 
420 82.2 








The gaseous products of the reaction were 
analyzed in a Burrell gas analysis apparatus, 
hydrogen being determined by combustion. 
Almost all the condensable substances (glyoxal, 
water, etc.) condensed out in the process of 
sampling. A small amount of material (always 
less than 1 percent) was absorbed in the KOH 
burette. It is very unlikely that this could be 
carbon dioxide, and it was assumed to be a trace 
of glyoxal or other condensable material. In any 
event the error which might be involved is 
negligible. In order to refer the gaseous products 
found back to the initial amount of glyoxal from 
which they came, it is necessary to know the 
amount of condensable material present in the 
mixture being sampled. This was done in the 
usual way by adding a known amount of nitro- 
gen to the glyoxal before the reaction, and then 
determining nitrogen in the gaseous products. 
The results of analyses of this kind are given in 
Table II. 

The main conclusion to be drawn from the 
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table is that the gaseous products of the reaction 
consist entirely of hydrogen and carbon monox- 
ide, the amount of carbon monoxide being much 
greater than that of hydrogen. With decreas- 
ing temperature the H2/CO ratio decreases to 
some extent. The carbon monoxide formed is 
approximately equal to the glyoxal disappearing. 

In the case of analyses at the half-way stage 
it is difficult to draw conclusions regarding the 
amount of condensable material formed, on ac- 
count of the uncertainty as to how much glyoxal 
has decomposed. At completion, however, it is 
possible to deduce the amount of condensable 
products with fair accuracy. At 450°C the 
amount is only slightly greater than that of 
hydrogen, while at 430° it is more than twice as 
much. The significance of this will be discussed 
later. 

It is possible from the above analytical data to 
draw up a material balance for C, H and O, and 
hence to deduce the composition of the conden- 
sable material. Too much reliance cannot be 
placed in such balances, since all the errors in- 
volved in the reaction measurements and 
analyses accumulate in the result. If such a 
balance is made for all the analyses at comple- 
tion, we obtain for the average composition of 
the condensable products 


Cy.00H1.2301.00 at 450°C, 
and 
Cy.00H1.2801.00 at 430°C. 
Or, on the average, 
Cy.00H1.2601.00. 


THE MECHANISM OF THE REACTION 


It has previously been assumed? that the 
glyoxal decomposition proceeds simply as 


CHO—CHO=CO+HCHO=2 CO+H:z, (1) 


formaldehyde being completely decomposed at 
the temperatures used. It is obvious that such a 
simple scheme cannot hold, since (a) carbon and 
tar are formed, (b) the pressure increase at com- 
pletion is only 90 percent instead of 200 percent, 
(c) a large amount of condensable products are 
formed, (d) the H2/CO ratio is only 0.2 to 0.3 
instead of 0.5. 
Other obvious possibilities are: 
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TABLE II. The products of reaction. 
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CHO—CHO=C+CO+H.0, 
CHO—CHO= (1/n)(CHO—CHO),, 


(2) 
(3) 


where (CHO—CHO), represents a polymer or 
condensable substance or substances of higher 
molecular weight than glyoxal, with approxi- 
mately the same empirical composition. 

It is not possible to account for the gaseous 
products on the basis of reactions (2) and (3) 
alone. On account of the fact that both CO and He 
are formed, it appears probable that reaction (1) 
occurs to a considerable extent. A simple com- 
bination of (1) and (2) could be made to yield 
the correct H:/CO ratio, but such a combination 
would give a pressure increase far too high 
(about 150 percent). A simple combination of (1) 
and (3) would give the correct pressure increase, 
but could not yield a He/CO ratio very different 
from 0.5. 

Provided that reaction (3)- gives rise to little 
gaseous material, a combination of reactions (1), 
(2) and (3) will account satisfactorily for the 
experimental results. Thus if we assume that all 
three reactions occur to the same extent, we have 
for the overall reaction 


3 CHO—CHO=3 CO+H2+C+H,0 
+(1/n)(CHO—CHO),. 


This accounts for the main facts, thus: 


(a) The CO formed is equal to the glyoxal decomposing, 

(b) The H:/CO ratio is 0.33, 

(c) If we assume that n is about 3, the pressure increase at 
completion will be approximately 80 percent, 


(d) The composition of the condensable product is 
C+H20+C2H202 = C;.00H1.3301.00, 
which is in good agreement with experiment. 


It appears therefore that the three reactions are 
of about equal importance. 

An investigation of the photochemical decom- 
position of glyoxal has been carried out.° 
The reaction was bimolecular, about 1.5 percent 
of the initial glyoxal decomposing to CO and He, 
and the remainder yielding CO and glycerosone 
(C3H,O;). This is essentially a polymerization, 
similar to reaction 3, but splitting out CO. 
Norrish and Griffiths state that glycerosone de- 
composes at 200°C. Hence it is not likely that it is 
of importance here. It is, however, possible to 
account for the experimental results on the basis 
of a mechanism such as: 


(1’) C.H,0.—2 CO+Hz, 
(2’) 2¢«C2H20.—-C3H,03;+ CO, 
(3’) 


Considering the fact that large amounts of car- 
bon are deposited, in addition to tarry materials, 
the first mechanism seems much more plausible. 

The diminution of the H./CO ratio with de- 
creasing temperature and the corresponding in- 
crease in the proportion of the condensable 
products must be due to a change in the balance 
between reactions 1 and 2, reaction 1 predominat- 


C;H,O;—polymer (tar). 


ing at the higher temperature. 


5 Norrish and Griffiths, J. Chem. Soc. 1928, 2829-40. 





STEACIE, 


THE RATE OF REACTION 


In view of the complicated nature of the reac- 
tion as outlined above, and of the fact that carbon 
and polymerized materials are formed, it is 
hardly to be hoped that the kinetics of the reac- 
tion can be very accurately established. Actually, 
the results obtained are decidedly erratic. It is 
possible, however, to infer the main characteris- 
tics of the reaction with some certainty. 

The complete data for a typical experiment 
are given in Table III. It will be seen that excel- 
lent first order constants are obtained. This was 


TABLE III. Data for a typical experiment. 








Temperature 440°C 

Pressure 

increase 
(cm) 


Empty reaction vessel 


Percent (first order), 
reaction min,~! 


0.0 

9.2 
18.4 
26.6 
40.5 
$1.5 
65.3 
74.4 
85.4 
92.8 
95.5 
97.3 


Pressure 
(cm) 


11.70 
12.70 
13.70 
14.60 
16.10 
17.30 
18.80 
19.80 
21.00 
21.80 
22.10 
22.30 
22.50 
22.60 
22.60 


Time 
(min.) 





0.041 
0.044 
0.045 
0.045 
0.045 
0.046 
0.046 
0.046 
0.046 
0.045 
0.045 
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invariably true throughout the investigation, in 
spite of the erratic nature of the absolute rate. 

The effect of temperature and of the glyoxal 
pressure on the rate of reaction is shown by 
Table IV. As a measure of the rate we have used 
Ti.5 and T2;, i.e., the times for the pressure to 
increase by 12.5 and 25 percent, respectively. 

On account of the erratic nature of the results, 
it is not possible to determine the activation en- 
ergy with any degree of accuracy. From the mean 
values of Ti2,; and 72; at higher pressures we ob- 
tain the value E= 54,000+10,000 cals. 

It will be seen from Table IV that the times 
for a given fractional pressure increase are 
roughly independent of the initial pressure, and 
hence the reaction is apparently first order, inso- 
far as the effect of pressure on the rate is con- 
cerned. The results, however, are too erratic to 
permit of accurate conclusions. 

The last 7 runs listed above were made in a 
vessel packed with short lengths of 3/16 inch 
Pyrex tubing, the surface-volume ratio being 
about three times greater than that for the empty 
bulb. It is apparent that packing the bulb has 
very little influence on the rate, and the reaction 
is therefore homogeneous. 


DISCUSSION 


As we have seen, the results are rather erratic, 
and the reaction is far from an ideal one from the 


TABLE IV. Reaction rate data. 








Initial 
pressure, 
(cm) 


Initial 
pressure, 
(cm) 


Ti2.5, 


(min.) 


T25, 


(min.) 


T 12.5, 


(min.) 


T25, 


(min.) 


Initial 
pressure, 
(cm) 


Initial 
pressure, 
(cm) 


T25, 


(min.) 


» -_ 
7 12. 5y T 25, 


(min.) (min.) 


T 12.5, 


(min.) 





Empty vessel 410°C 
10.8 0.95 

15.8 

20.3 


20.1 


Empty vessel 
A 2.1 
a 2.4 





16.55 
13.30 
13.05 
11.70 
10.30 

9.90 


3. 7.95 
: : 6.30 
: 4.75 


- 
1. 
i. 
‘ 
E 
Ly 


Empty vessel 
22.15 y d 
21.70 
15.50 


10.40 


3.30 


Reve 


a 430°C 
13.70 
11.35 


U 


is.) 


Packe 
18.55 : 
15.10 


mio 


Packe 
12.45 d 
10.60 


mn & nw & NN dK — 
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kinetic standpoint. The main interest in the re- 
sults, however, lies in the fact that previous as- 
sumptions regarding the simple nature of the 
glyoxal decomposition have been shown to be 
incorrect. 

In the oxidation of acetylene no appreciable 
quantity of carbon or tar is formed.': ° It there- 
fore follows that while glyoxal may be formed 
and oxidized in the acetylene oxidation, it cannot 
decompose during the course of this reaction as 
has been previously assumed. 

It also follows that the scheme proposed by 


6 Steacie and McDonald, unpublished work. 
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Lenher for the oxidation of acetylene catalyzed 
by the oxides of nitrogen, viz.: 


C.H,+0,.= CHO—CHO—(CHO— CHO); 
\ CO+HCHO 

| +0, 
CO+CO,+H,0, 


is not a possible mechanism, although certain ad 
hoc modifications could be made to bring it into 
line with the present work. 

An investigation of the glyoxal oxidation is in 
progress. 
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Orientation Effects in Bimolecular Ionic Reactions 


JuLian M. Srurtevant, Sterling Chemistry Laboratory 
(Received January 15, 1935) 


It is shown on the basis of accepted views concerning second order ionic reactions in solution 
that the variation with concentration of the electrostatic orientation effects in such a reaction 
arising from the lack of symmetry of a spheroidal ion is negligible in dilute solution, Devia- 
tions from Brénsted’s law must be attributed to other causes. 


T has been observed that certain bimolecular 
ionic reactions deviate very considerably from 
the requirements of the limiting Brénsted for- 
mula, which is, for water as solvent and at 25°C, 


log k”’ =log ko+2122(u)!. (1) 


In this equation k”’ is the observed bimolecular 
specific reaction rate, ko is a quantity indepen- 
dent of concentration, 21, zz are the valences of 
the reacting ions, and y is the ionic strength. 
In many instances these deviations are probably 
due to failure of the simple Debye-Hiickel limit- 
ing law for ions of high valence, but there are a 
few cases where this explanation cannot hold. 
Outstanding among these is the one observed by 
LaMer and Kamner.! These authors found that 
the reaction between a-brompropionate and 
thiosulfate ions follows the Brénsted formula ac- 
curately, while the analogous reaction involving 
the 8-brompropionate ion leads to a slope of the 


‘LaMer and Kamner, J. Am. Chem. Soc. 53, 2832 (1931). 


log k’’ vs. (wu)? curve which does not even have the 
expected sign. After an exhaustive study had in- 
dicated that the effect is not due to side-reactions 
or other such causes, LaMer and Kamner sug- 
gested that purely electrostatic effects resulting 
from a comparatively large separation of charge 
and center of reaction (the bromine atom) in the 
unsymmetrical organic ion might lead to the 
observed results. It is the purpose of the present 
note to examine this suggestion a little more 
closely. It will appear from the study of a fairly 
general case that this explanation is not valid, 
if one accepts present views concerning second 
order ionic reactions in solution. 

Christiansen’ developed a statistical theory of 
bimolecular ionic reactions which he showed to 
be in formal agreement with Eq. (1) in very di- 
lute solutions. His treatment led to a general 
result which may be written 


k'’ =ky exp (— ezeg1/kT)-Q-(1+2 x3). (2) 


? Christiansen, Zeits. f. physik. Chemie 113, 35 (1924), 
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The term in parentheses expresses an approxima- 
tion to the effect of deactivating collisions, and 
may be considered as practically independent of 
concentration at low total concentrations. Q is a 
factor introduced to allow for orientation effects, 
and was also considered by Christiansen to be 
independent of concentration in very dilute 
solutions. The quantity ¢g; is the electrostatic 
potential due to an ion of type 1 and its atmos- 
phere, at a distance 7 from its point charge €21; 
r is to be taken as the separation of the charges of 
the ions at the moment of reaction. Christiansen 
used for the value of ¢; the value obtained by 
Debye and Hiickel* for spherically symmetrical 
ions of negligible diameter in a solution of low 
ionic strength, 


¢g1= (€2;/Dr)e—" = — ez\x/D+(€2;/Dr)(xr1). (3) 
D is the dielectric constant of the medium and x 
is defined by 

Kk? = (mre®N/125DkT)y, (4) 


N being Avogadro’s number. 
Eq. (2) will be used as the basis of an extension 
to the case that the ions of type 1 are not spheri- 
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cally symmetrical. In order to make the treatment 
reasonably simple the ions of type 1 are assumed 
to be prolate spheroidal, with semiaxes a and 
b (a>b). These dimensions define a spheroidal 
surface beyond which, on the average, the charges 
of the other ions in the solution do not penetrate. 
A coordinate system is taken with origin at the 
center of an ion of type 1, the semimajor axis of 
which lies along the z-axis. Then in the confocal 
coordinates, y, 3, ¢, defined by the relations 


Oc n<a, 
O< dcr, (5) 


x=c sinh 7 sin 3 cos ¢ 
y=c sinh 7 sin 0 sin ¢ 


z=c cosh 7 cos 3 O0<¢ g<2r, 


the “‘surface’’ of the ion corresponds to n = no such 
that 


cosh no=a/c; sinh no=b/c 


where 2c is the separation of the foci. The point 
charge €2; is located at (n’, 3’, ¢’). 

If ¥,’ is the potential inside the ion due to its 
atmosphere and y,”’ is that due to the charge e7,, 
then the total potential for n<mo is Yi=y,' 
+y,’’. Now yi’ satisfies Laplace’s equation and 
is bounded at the origin, so that it has the form‘ 


y= > > P.™(cosh )P»”(cos 3)[Cum cos me+Drun sin me | (6) 


n=0 m=0 


where the P,,” are the associated Legendre functions of the first kind and the Cum, Drm are constants. 
Expansion’ of the reciprocal of the distance between (n, 3, ¢) and (n’, 3’, ¢’) gives 


Wi" =D LV On™(cosh 7) Pn™(cos I)Enm cos me, 


n=0 m=0 


Eno= (2n+1)(ez;/Dc)P,,(cosh n’)P,(cos 8’), 


Enum = (—1)"2(2n+1)(€2:/Dce)[(n—m)!/(n+m)!]?P»”(cosh ’)P.™(cos 3’) cos mg’ 


if we take y’ =0. Here the Q,.” are the associated Legendre functions of the second kind and D is the 
dielectric constant inside the ion, assumed to be the same as that of the surrounding medium. 
The potential ¥2 for 7>n0 is a solution of Ay2= «2. Assuming that Y2= 7(n)O(8) ®(¢), we obtain 


=acosme+8 sin mg, 


(1—x?)y"” — 2xy’ + (k?2x?+u—m?/(1—x*))y =0. 


m=0, 1, 2, °°: (8) 
(9) and (10) 


In Eq. (9) y=0, x=cos 3 and —1=x~1; in Eq. (10) y=H, x=cosh n, and 1=x< @. m and » are 


’ Debye and Hiickel, Physik. Zeits. 24, 185 (1923). 


* Cf. Hobson, Spherical and Ellipsoidal Harmonics,Cambridge, 1931, p. 412. 


> Hobson, reference 4, p. 416. 
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separation constants. Since the parameter k = xc is small, (9) and (10) may be treated by perturbation 
methods, taking n= n(n+1)+0(k?). We thus set 


0," =) ¢;P;7"(é), £=cos v, m=(, a 3 en. (11) 
Using the formula 


(n—m+2)(n—m+1) 2n?+2n—1—2m? (n+m)(n+m—1) 


er = pr cae _P m oa anenennaanan a 


(2n+1)(2n+3) “= (2n—1)(2n+3) ' (2n—1)(2n+1) 








we readily find that, aside from constant factors 
(n—m+2)(n—m-+1) (n+m)(n+m—1) 
7 42 “Pe, | +018) 
2(2n+1)(2n+3)? 2(2n —1)?(2n+1) 


= Pita ny2P nyo tB"n—2P™ n-2 + O(k"*). 








(-) ~” _ P.™ a e 
(12) 


An approximate solution of Eq. (10) can now be found by consideration of the unseparated equation 
involving y= O// as dependent variable. This equation is satisfied by 


y= sin” v sinh” n e—* cos § cosh 7 = K(d, n), 
so that we have 


.n(¢) = { K(é, )O."(8)dé,  ¢=cosh 2. 
vi 


We thus have to evaluate integrals of the type 
a” =1"(5? — yn f (€2—1)™/2e-*SEP™(E)dE. 
1 


When ¢ is not real and between —1 and +1, 


(n+m)! 1 q"-™ 
= (¢?-—1)-™/? — (€?—1)" 
(n—m)!2"n! ian 


—m 





Introducing this relation into Eq. (14) and integrating m times by parts gives 


(n+m)!(¢?—1)™? po 
." = 4" f e— "EP, (E)dé. 
(n—m)! (k¢)™ 1 





It can be shown that 


@ 


inth pry} 
{ e-*StP (t)dt= (<) HI n4\(tkE), 
Jy (Ro)iX2 


where J7,,,; is the Hankel function of the first kind,® and enters here because it is the only type of 
Bessel function which vanishes for infinite positive imaginary argument. 

To complete the calculation of the solution of Eq. (10) it is necessary to extend the computation of 
the coefficient a@”,42 in Eq. (12), since the fact that H,.,(tko)/H™ 4;(ik¢) =O(1/k*) necessitates 
the inclusion of more terms to obtain the desired accuracy. The result obtained is, omitting unessen- 
tial constant factors, 


(¢2—1)™/? i 
H,."(¢) = ppt [EO elie +E (— 1)4a™,, 2H sain ike) | +relative O(k*), 
é m+} j=1 


* The definition of H,,,; used here is that given in Jahnke-Emde, Funktionentafeln, 2nd ed., Berlin, 1933, p. 202. In 
€valuating the potential about a spherical ion the point charges of which are not located at its center, Kirkwood (J. 
Chem. Phys. 2, 599 (1934)) was led to the functions eK (xr)r~"— where the K,, are polynomials. These functions are 
identical, aside from constant factors, with the functions r- 1H, 44 (ixr). 
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(n+-m+1)(n+m-+2)-++(n+m+2)) 

2+4+++(27)(2n+1)(2n+3)2(2n+5)2- + +(2n+2j+1)2(2n+2j+3)+++(2n+4j—1) 
2 jk? (4m? —1) 

f (2n —1)(2n+3)?(2n+4j+3) 


a”... 2)= Ri 








| relative O(R*). (15) 
The general expression for the potential outside the ion is taken in the form 
Yo= >> Dd An™(cosh 7)O,"(cos 3)[A nm COS ME+Bnum sin me |. 
n=0 m=0 


The coefficients A,» and Brn may be calculated from the boundary conditions that 


Yi=Ye and VYyi=Vye2 at n=N7o. 
We thus find that 


rd o* rd O”...2 d i” .0) fd QO" 242 d TT” ns2 
dt P,.™ d¢ P*,_- d¢ P,”™ d¢ FP” ne d¢ P,™ 


_— m Fe 
An +n—2, m4 


d H,™ d I me ad I ni” ‘ d I or 














dé Pm 5 dt Pr, 9 dé P,.™ de p*.,, 12 
+relative O(k*), Bym=relative O(k*). (18) 




















In order to simplify the formulae somewhat we will specify that the charge ez; is located on the 
Z axis at a distance p from the origin, so that cosh n’=p/c, cos 8’ =1 or cosh n’=1, cos 8’ =p/c, 
according to whether p>c or p<c. There is no significant loss in generality in this restriction since we 
are mainly interested in the extreme case of a long thin ion in which the charge could not be far 
removed from the major axis. In this case the potentials are independent of the coordinate ¢, so that 
m is to be set equal to zero in all the above formulae. 

Now if c=3X10-* cm and we restrict ourselves to the consideration of solutions in which «!=0.1 
we may write that 

1 


1-3-5+++(2n—1) /2\3 1 2\'1—k¢ 
11 ,\(488) = (=) @orm, mGes)~=(-) 


7 . 1? 
isi " Na? (ks)! 





(19) 
(n+1)(n+2)-++-(n+27) 


27. 
a 


Qn, j;=k? 
OT Dede e+ (29)(Qm-+1) (2+ 3)*(2M-+5)* + (Qn 2741)? 27+3)- ++ (Qn+4j—1) 





The relative error in these approximations is less than 0.5 percent if ¢<1.2. We thus obtain the simple 
result that . 
A,(S)=k-"'0,(6), Hol S) = —1+k-'0o(5), (20) 


so that the expression for ~2(fo, 3), omitting certain negligible terms, is 


> =, | 


€21K p 
v2(So, d) a a ef a > Qn+t)P,(* 
D D anil | 


ont) 


k? n(n—1) p d Qn-2 

ere) ——-P y_»2 gin n— 2—] P's 2 "gad 
+ Gn—1) (‘)|2 sitet di P, Jer 
ke (n+1)(n+2) 
2 (2n+3)? 


Onse 





P : +(¢?—1)(P,)? : 
wo(*)| On g?—1)(P,) — P. 


| Pa(cos 3). 
f=fo 
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It can be shown by using the asymptotic 


values 
rT 4 en (nth)n 
Q,(¢) - (“) Age, 


n/ (2sinh n)} 


1 e(rt))n 
P,(¢) = 





, m>I1 
(n7)} (2 sinh n)} 
valid for ¢=cosh 7 >1, that the relative error in- 
troduced into each term of the sum in (21) by 
omitting the k® terms vanishes with increasing n. 
Moreover these errors are certainly monotonic 
functions of m for values of {9 and p/c sufficiently 
near unity and for greater than 5 or 6. We can 
thus obtain an estimate of the total error in- 
volved in neglecting these terms by examining the 
first few. Taking {9=1.1 and p/c=1.05, we find 
that the error is certainly considerably less than 
1 percent. This value of {9 corresponds to an ion 
having a major axis about 23 times as long as its 
minor axis. Thus even in quite extreme cases we 
may write with sufficient accuracy, 


Valo, 8) = —eeyx/D+(e2,/D) 5 (2n+1) 


n=0 


XP alp/c)Qn(Fo)Pa(cos 8); (22) 


the first term is the potential due to the ionic 
atmosphere and the second term is that due to 
the charge ¢z;, so that we obtain again the simple 
Debye-Hiickel result. 

To evaluate the orientation factor 2 appearing 
in Eq. (2) we assume that the irreversible chem- 
ical reaction between an ion of type 1 and one of 
type 2 (the latter being spherically symmetrical) 
can take place only when the charge of the ion 2 
is located on the spheroid defined by 7= 70, and 
the line joining the centers of the ions is con- 
tained within a solid angle Aw determined by the 
finite reactive portion of the surface of the ion 1. 
This picture should serve as a simplified repre- 
sentation of a reaction between a comparatively 
simple inorganic ion and a complicated organic 
ion, the reaction involving one particular group 
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present in the latter. Then Q is the fraction of 
the ions of type 2 with charges located on the 
spheroid which have the orientation necessary 
for reaction. If we assume a Boltzmann distribu- 
tion’ we have 


"a ¢!” 
f f exp ( - eZopo/kT)dS 
oe’ “e’ 


ff exp (— eSope/kT)dS 





Q(no) _ ’ (23) 


where the limits 3’, 3’, ¢’, ¢’’, correspond to Aw, 
dS is the element of surface, and k is Boltzmann's 
constant. 

It is at once apparent that a potential of the 
form of (22) will lead to an Q which is indepen- 
dent of x. Furthermore, if we use 2 in place of ¢; 
in Eq. (2) we obtain, after introducing appro- 
priate numerical values, 


d log k’’/d(p)'=2120, (24) 


again in agreement with Brénsted’s formula. 

As might be expected, exactly the same result 
is obtained in the consideration of the simpler case 
of a spherical ion whose charge is not located at 
its center. It seems quite safe to draw the con- 
clusion that no reasonable type of dissymmetry 
could lead to purely electrostatic orientation 
effects that would be observable at concentra- 
tions within the range where the Debye-Hiickel 
theory may be applied with any degree of confi- 
dence. The observations of LaMer and Kamner' 
on the 6-brompropionate-thiosulfate reaction can 
only be reconciled with the theory of Brénsted 
by assuming that the log k’’ vs. u* curve has a 
maximum just below the last experimental point. 

In conclusion, it is a pleasure for the author 
to express his gratitude to Professor Lars On- 
sager for many helpful suggestions given in con- 
nection with the treatment of the differential 
equations in this paper. 


7 See Christiansen, reference 2, p. 43. 
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Statistical Mechanics of Fluid Mixtures* 


Joun G. KirKwoop, Department of Chemistry, Cornell University 
(Received March 1, 1935) 


Expressions for the chemical potentials of the com- 
ponents of gas mixtures and liquid solutions are obtained in 
terms of relatively simple integrals in the configuration 
spaces of molecular pairs. The molecular pair distribution 
functions appearing in these integrals are investigated in 
some detail, in their dependence upon the composition and 


density of the fluid. The equation of state of a real gas 
mixture is discussed, and an approximate molecular pair 
distribution function, typical of dense fluids, is calculated. 
Applications of the method to the theory of solutions wiil 
be the subject of a later article. 





I. 


LTHOUGH statistical mechanics has met 

with considerable success in the inter- 
pretation of the properties of rarefied gases as 
well as those of crystalline solids, serious prac- 
tical difficulties have retarded its application to 
dense gases and to liquids. In the present article, 
we shall describe a statistical treatment of fluid 
mixtures, which is formally applicable both to 
real gases and to liquids. It has the advantage of 
yielding expressions which lend themselves more 
readily than the usual statistical formulas, to 
approximations suitable for condensed systems. 
While we make no pretension to a resolution of 
the difficulties presented by the liquid state, the 
more serious of which resist more than a super- 
ficial treatment, we believe that the method to 
be described throws interesting light upon the 
problem, and that it may prove useful in opening 
up new lines of attack, particularly in the field of 
liquid solutions. 

We must first ask whether classical statistical 
mechanics is competent to describe the depend- 
ence of the thermodynamic functions of fluids 
upon intermolecular forces. At ordinary tem- 
peratures, there is good reason to believe that 
the classical approximation is adequate for most 
gases and liquids consisting of chemically satu- 
rated molecules. When certain conditions involv- 
ing the molecular masses, moments of inertia, 
and the intermolecular forces are satisfied, the 
partition function of a system of molecules is 
closely approximated by the product of a func- 
tion of the temperature alone, and the Gibbs 


* The material of Section V of this article was presented 
at the eighty-ninth meeting of the American Chemical 
Society held in New York City during the week of April 
22, 1935. 
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phase integral in the phase space associated with 
the low frequency translational and rotational 
motion of the molecules.! The validity of the 
assumptions involved in this approximation is of 
course more questionable for liquids than for 
gases, particularly when there is strong coupling 
between the rotational motion of neighboring 
molecules as in polar liquids. We shall provi- 
sionally assume the adequacy of the classical 
approximation, remembering that our results are 
applicable only when this is true. 

When the potential of the intermolecular force 
is of special form, for example a sum of terms 
depending upon the relative coordinates of 
molecular pairs, the explicit use of unwieldy 
integrals in a many dimensional phase space may 
be avoided. We shall show how the statistical 
formulas for certain important thermodynamic 
functions of a fluid mixture may be recast in 
terms of comparatively simple integrals in the 
configuration space associated with the relative 
motion of a molecular pair. The molecular pair 
distribution functions, which appear in these 
integrals are, to be sure, formally defined as 
ratios of many dimensional phase integrals. 
However, by certain transformations, they may 
be expressed in a form which leads to reasonable 
approximations both in condensed and rarefied 
systems. The approximation for rarefied systems 
leads to a novel approach to the problem of the 
equation of state of a real gas, while the approx- 
imation for condensed system promises to be of 
use in the theory of liquid solutions. On the 
other hand, the approximate distribution func- 
tions which have thus far been obtained do not 


1E. Wigner, Phys. Rev. 40, 749 (1932); Uhlenbeck and 


Gropper, Phys. Rev. 41, 79 (1932); J. G. Kirkwood, Phys. 
Rev. 44, 31 (1933); J. Chem. Phys. 1, 597 (1933). 
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appear to be sufficiently exact for the calculation 
of the equations of state of liquids. 

The method of the present treatment is based 
upon the same principle as the charging processes 
employed in the theory of electrolyte solutions. 
As Onsager has clearly pointed out, parameters 
which appear in the potential of the inter- 
molecular forces have essentially the same status 
as the parameters of external force,? which 
determine the mechanical interaction between 
a thermodynamic system and its surroundings, 
and may be manipulated in the same manner. 
The application of this principle is by no means 
restricted to intermolecular forces of the Coulomb 
type, in the potential of which ionic charges 
enter as natural parameters. Arbitrary param- 
eters may be introduced into the potential of 
intermolecular force in any system of molecules, 
by means of which the coupling between the 
molecules may be continuously varied in any 
convenient manner. This device allows one to 
express the chemical potential of any component 
of a fluid mixture in a particularly simple form. 
It also proves to be very useful in investigating 
the dependence of molecular pair distribution 
functions upon the composition and density of 
the fluid. 


II. 


The following analysis is formally applicable 
to all homogeneous fluids, whether real gas mix- 
tures or liquid solutions. Differences arise only 
in the method of approximating distribution 
functions. We therefore consider a homogeneous 
fluid of » substances, in thermodynamic equi- 
librium. By suitably restricting the ranges of 
composition and pressure so as to exclude coex- 
istent phases containing any of its components, 
we may suppose that the fluid completely fills 
a vessel of volume v, which confines the motion 
of the system to a finite region of configuration 
space.* We denote the absolute temperature by 

*L. Onsager, Chem. Rev. 13, 73 (1933); J. W. Gibbs, 
Elementary Principles in Statistical Mechanics, Collected 
Works II, Chap. VII (Longmans). 

*We postulate the absence of coexistent phases con- 
taining any of the components of the fluid, since we wish 
to represent the system by a petit ensemble, each example 
of which contains the same number of molecules. Actually 
this is a trivial restriction, for Gibbs has shown that 
although an open system is represented by a grand 
ensemble, its equilibrium properties do not differ sensibly 


from those of a closed system represented by a petit 
ensemble, each example of which is made up of the average 


T, the numbers of molecules of the several com- 
ponents by Nj, ---, N,, and the total number 


of molecules of all species, > N., by N. Ac- 
sal 


cording to classical statistical mechanics, we 
represent the system by a canonically distributed 
petit ensemble, and obtain the usual expression 
for the Helmholtz free energy, Fy. 


di 3(T)%s 
a? -| — ~ lex 
s=1 N, ! 


Bu =|. ‘ fe 8Vndw,-+-dwy, (1) 


where 8 is equal to 1/kT and dw; is a differential 
element of the configuration space associated 
with the low frequency rotational and trans- 
lational degrees of freedom of the molecule 7. 
The integration in the sub-space of each molecule 
is extended over a region bounded by the volume 
of the vessel containing the fluid. In general 


[eros (2) 


where o; is equal to 4m for diatomic molecules 
which have only two low frequency rotational 
degrees of freedom and to 82? for polyatomic 
molecules which in general have three. The 
potential of the intermolecular forces, Vy, is a 
function of the configuration coordinates of all 
N molecules. The functions f,(7), depending 
upon the temperature alone, are products of 
certain momentum integrals associated with the 
low frequency degrees of freedom of a molecule 
of type s and the partition function associated 
with its high frequency vibrational degrees of 
freedom. For our present purposes we do not 
need to characterize these functions more closely. 
The chemical potential per molecule of a com- 
ponent i of the fluid is thermodynamically 
defined by the relation 
ui=(OF/ON;,)T, ¢, Ni, coel¥ pw fo (3) 
It is statistically defined by Gibbs‘ in terms of 
the grand ensemble. For our purposes, it is more 


convenient to define y; in terms of the petit 
ensemble upon which Eq. (1) is based. This may 


number of molecules of a system of the grand ensemble. 


In case coexistent phases exist, we simply regard v as the 

volume of the homogeneous part in which we are interested. 
‘J. W. Gibbs, Elementary Principles in Statistical 

Mechanics, Collected Works II, p. 187, Longmans. 
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be done in a very simple manner. We first observe 
that the treatment of N as a continuous variable 
in the thermodynamic definition (3) is simply a 
convenient mathematical device. A_ physical 
definition of chemical potential requires the use 
of finite differences, AF and AN;,, in the ratio 
AF/AN;,. For values of AN; which do not sensibly 
affect the composition of the fluid, this ratio will 
not differ sensibly from the derivative on the 
right-hand side of Eq. (3). We are therefore 
justified in defining y; as the change in free energy 
produced by the subtraction or addition of a 
single molecule of type i from the macroscopic 
system, while maintaining the temperature, 
volume and the numbers of molecules of the 
other species constant. 


= F(M,, = re Ni, -+-N,, T, v) 
— F(N,, --+, Ni-1,---N,, T,v). (4) 


The fact that the removal of a single molecule 
is scarcely a physically realizable process does 
not alter the value of the definition (4) when 
used in conjunction with the petit ensemble, each 
system of which is supposed to comprise exactly 
the same number of molecules. 

From Eggs. (1) and (4), we obtain the following 
expression for y;: 


Ki= 


f(T) Z(Ni, 
-#T log| 


N; Z(N,,°: 





we 
. to 
>, N,—1, ---N,, T,v) 


Further simplification depends upon the form of 
the potential of the intermolecular forces, Vy. 
We shall treat here only systems, for which Vy 
may be expressed as a sum of terms Vj), each 
depending only upon the relative coordinates 
of a molecular pair (R/).* 


* It is to be remarked that the potential of Eq. (6) in no 
way excludes the simultaneous interaction of groups of 
more than two molecules. It simply states that in such a 
group, the mutual potential energy of any pair is inde- 
pendent of the presence of the other molecules. For example 
the electrostatic energy of a system of point charges is of 
this type. London has shown that the potential of the 
attractive van der Waals forces between chemically 
saturated molecules satisfies the requirements of Eq. (6). 
[F. London, Zeits. f. physik. Chemie 11, 222 (1930).] 

No such clear cut statement can be made regarding the 
repulsive forces, which operate at small distances, de- 
termining the effective collision diameters of molecules. 
However, it is plausible to assume that the potentials of 
these forces are approximately additive, and even if they 
are not exactly so, it is only necessary for our purposes that 
their narrow range of action is not sensibly displaced. 
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Vx=>d View (6) 
kel 
We now define a fictitious potential of the fol- 
lowing type 
Vy (a, ++ Aw) =D MA ed, (7) 
k<l 
where \i-*-*Ay are arbitrary parameters. By 
varying these parameters between zero and 
unity, we can continuously vary the coupling 
between the N molecules from zero to its full 
value. In particular, we notice that when a 
value unity is assigned to each of the N param- 
eters \, we have 


Vy(1, ee 1)= Vy, (8) 


where Jy is given by Eq. (6) and represents the 
actual coupling between the molecules. If, on the 
other hand, a value unity is assigned to N-1 
of the parameters \ and a value zero to a single 
\;, we have 


Vv, oe 0, -*-1)=Vy-i, (9) 


where Vy_; is the potential in a system of 
N-—1 molecules formed by the removal of a 
single molecule 7. Corresponding to Vy(Ax, ++ Ay), 
we define a phase integral, 


Zy(A1, +**Aw) 


= foo ferervor aden: daw. (10) 


When each of the parameters, A, is assigned the 
value unity, we have the following identity : 


Zn(I, -+-1)=Z(M, ---, Ni, -++N,). (11) 


When N—1 of the parameters \ are equal to 
unity, while a single \; corresponding, let us say, 


In general, we should expect terms depending upon the 
relative coordinates of three or more molecules to be 
important only if the total intermolecular field is strong 
enough to produce a large perturbation in the electronic 
structure of a molecule in the fluid. Theoretical con- 
siderations suggest that not only in gases but also in 
liquids, this is not the case. Empirical arguments against 
any appreciable distortion of the electronic structure of 
molecules in most liquids are furnished by spectroscopic an 
optical evidence. For example, the molecular refractivity 
of most liquids is substantially the same as that of their 
vapor. 

There are to be sure some types of intermolecular 
potential which do not satisfy Eq. (6), for example those 
of chemical forces of the homeopolar type. Another 
example is furnished by the energy of polarization of 4 
molecule by a group of ions, 
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to a molecule of type 7 is equal to zero, we may 
write 


Zy(1, tee, ), eee 


1) 
=f. +f eA s-rdey --dwx. (12) 


Now since Vy_; contains no terms depending 
upon the coordinates of the molecule 7, the 
integration over its configuration space in Eq. 
(12) may be carried out immediately. Performing 
this integration and making use of Eqs. (1) and 
(2), we obtain 


Zy(A, tee, Q, ee -1) 


=0Z(Ni, ---,Ni—1,---N,). (13) 
From Eqs. (11) and (13) we may write 

Z(N1, --+, Ni, «++ N,) 

Z(Mi, --+, Ni—1, «+ -N,) 

Zx(1, eos, 1, +++) 


0 iv ° 
Zy(A, ee, 0, ++) 








Further we have the identity, 
Zw(1, -+-, 1, +1) 

log 
Zy(1, «++, 0, +1) 





{ 0 log Zv(Q\0) 


is (15) 


0 On; 


To simplify notation, we write Zy(A,) for 
Zy(1,+++,X;, +1). In the future this conven- 
tion will be adopted for all functions of the 
parameters A, only those differing from unity 
. being explicitly indicated. Returning to Eq. (10), 
we observe that 


Zv(n)= f - . | fears day dwn, 


N 
Vi=d Vin. 
k=1 


(16) 


Differentiation of (16) with respect to ; yields 
0 log Zy(di)/Or =< B Vi(A,), 
f- ° f Vie FIV nit iV dy, - ° -dwy 


Sons fe~PlYn-t Vida, + + «dwn 


—_—_—- 


Vi(X,) = (17) 
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Combining Eqs. (5), (14), (15) and (17), we 
obtain for the chemical potential : 


sl 
m= kT log Ni/ot [ Vid)dd+e(7), 
0 
| (18) 
¢i(T) = —kT log o: f(T). 
The mean potential energy V;(A;) can be ex- 
pressed as follows: 


N 
V (A;) = a Vin(rs), 
=1 


k 


fee fVige 8 Yn at iVilday,+ + -dwy 





Vir(A3) =- 4 (19) 


foes fe-BlYn-t Milde + den 
Since V;, depends only on the coordinates of the 
molecules 7 and k, we may write, 


1 
Vidi) = ff V p08 ¥ OD d wid, (a) 
0 50 }.U" 
(20) 


0; ).U" ee + fe-BIVn-1t iV ildyy,- ° ‘-dwy_s 


f- . + fe-BIVn-1t MV il dy) - ° -dwy 


(b) 


The function W;*(\;), as we shall presently show, 
is the potential of the mean force acting between 
the molecules i and k. At present we merely 
observe that in the absence of external fields of 
force W;*(X;) can depend only upon the relative 
coordinate of the molecule pair 7 and k, except in 
a region of negligible volume near the surface of 
the fluid. 

We may therefore refer the position of the 
center of gravity and the orientation of the 
molecule k to a coordinate system fixed in 
molecule 7. Both Vi, and W;*(A;) will depend 
upon these relative coordinates alone. Moreover, 
the Jacobian of this transformation is unity. If 
we then integrate over the configuration space of 
molecule i, the expression for Vix(A;), Eq. (20), 
becomes : 


eBWkOD = 





Fant = (1/0) f Vige~P¥ FN O0dw,, (21) 


where dw; is understood to refer to the coor- 
dinates of molecule k& relative to a system of 
axes fixed in molecule 7. Introducing Eq. (21) 
into Eq. (18), and remembering that Vjx(A;) 
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is the same function for all molecules k of the 
same type, we may write: 


: oa 
w= kT| log N;/v+ > — Bulrar,| +¢(T), 
k=1 VY 0 
(22) 
Bi(d;) = (1/oukT) f V pe 8 ¥ FO wy, 


where the summation extends over all molecular 
types k present in the fluid, and the function B;,, 
refers to a single pair of molecules of types 7 and 
k, respectively. Each B;, is a function of tem- 
perature, volume and composition, since W;*(A;) 
is a function of all of these variables. The nota- 
tion W;*(A;) means of course that all the param- 
eters \ have the value unity except A; corre- 
sponding to a single molecule of type 7. 

The mean energy E of the fluid may be ex- 
pressed quite simply, if we assume classical 
equipartition with respect to the low frequency 
rotational and translational degrees of freedom 


E=E,(T)+ Vy, (23) 


where Vy is the mean potential energy associated 
with the intermolecular forces and E,(T), a 
function of the temperature alone, is sum of the 
internal energies of the molecules and their 
mean translational and rotational kinetic ener- 
gies. Writing Vy in the form of Eq. (6), and 
remembering that V;, depends only on the 
relative coordinates of molecules k and 71, we may 
write* 


b= 43D (NiN;/2) Vim 
irk 
=1 


Vw=2 Vi 
i<k 


(24) 


Va = (1/ow) f Vipe-8¥ Fd e,. 


We note that W;* appearing in Eq. (24) is 
W;*(1), the value of W;*(A;) when molecule 7 
has its full coupling with the other molecules of 
the system. Finally we may write, 


E=E,(T)+3kT & (NiN;/v)Bix. (25) 

isk 

=! 
The quantity B;, in Eq. (25) is equal to B;,(1) 
of Eq. (22) when ; is unity. Eo(T) may be 


* Eq. (24) has been given by Hildebrand and Wood, 
J. Chem. Phys. 1, 817 (1933). 
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explicitly calculated if the internal partition 
functions of the molecules are known. We shall 
not be interested in this point here. 


Ill. 


The use of Eqs. (22) and (25) for purposes of 
calculation depends upon an investigation of the 
functions W,*(A;) which determine the prob- 
ability distribution of the pair of molecules i 
and k relative to each other. It may readily be 
shown that W;*(A;) is the potential of the mean 
force acting between the molecules 7 and , 
averaged over all configurations of the remaining 
N-—2 molecules of the system. The force acting 
at the center of gravity of molecule z for a fixed 
configuration of the other molecules, is 


F,Q;) = —ViVy(Qai), 
where V; is the gradient operator associated with 
the coordinates of the center of gravity of mole- 


cule 7. If we take a mean value with the pair of 
molecules 7 and & held fixed, we obtain 


(26) 


ik 
“(A,) 
f- ‘ fV iV (dje8Y 4 O0de@): ‘ -dwn_s 
_ : a » 
foes fe-BYNODd a, + -dwy—s 

Differentiation of both sides of Eq. (20b), which 
defines W;*(\;) and comparison with Eq. (27) 
yields 





1, ————_— 
F;(\;) = —ViW* (Ai). (28) 


The function W;*(\;) is therefore the potential 
of the mean force acting at the center of gravity 
of molecule 7 when is held fixed at some point 
in its vicinity. Similarly, it may be shown that 
W;*(\,) is the potential of the mean generalized 
forces associated with the relative orientation of 
i and k. If the relative orientation is specified 
by the Eulerian angles 3, ¢ and y, the mean 
torque tending to produce rotation through the 
angle # is —dW;*(d;)/dd8. These relations have 
been used by Einstein, Schmoluchowski and 
others,® who showed that the simple Boltzmann 
distribution is valid provided the potential of the 
mean force is employed in the exponential factor. 

It is interesting to remark that the potential 
W;*(1), corresponding to full coupling between 
molecule i and the other molecules, is susceptible 


6 For a discussion of this point, see L. Onsager, Chem. 
Rev. 13, 73 (1933). 
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of experimental determination, at least in its 
dependence on the relative coordinates of the 
molecular centers of gravity. According to the 
theory of x-ray scattering by liquids, developed 
by Prins and Debye,’ the angular distribution of 
intensity of the scattered radiation is a sum of 
mean values of functions depending only upon 
the relative coordinates of pairs of molecules. 
Mean values of such functions are determined 
by the probability distribution function e-6”:* 
in a system in thermodynamic equilibrium. By a 
Fourier integral inversion of the theoretical 
scattering formula, (e~*”:*—1) may be com- 
puted from the observed intensity measured as 
a function of the scattering angle. By Eq. (25), 
we see that a knowledge of W;*(1) suffices for 
the calculation of the mean energy of the fluid 
from the potentials V;,. However, it does not 
suffice for the calculations of chemical potentials 
or the Gibbs free energy, ¢, which is equal to 


> Nw. For the latter purpose we must know 


i=l 

W;*(\;) for all values of \; between zero and 
unity. This information is not furnished by 
x-ray scattering experiments. Thus, while they 
yield valuable information concerning the sta- 
tistical structure of liquids, they do not provide 
sufficient data for the calculation of entropy and 
free energy. 

We must consequently depend upon theoret- 
ical information about the potentials of mean 
force W;*(A;). Since the many dimensional phase 
integrals of Eq. (20b), which defines W;*(\,), 
are inconvenient to handle, we follow an indirect 
method of attack, leading to reasonable approx- 
imations both for condensed and _ rarefied 
systems. Logarithmic differentiation of both 
sides of Eq. (20b) yields 





aW.H(d)/As=" Vi)-Vid), — (29) 


where 
th frees f Vie PVN Mild, - - -deoy_o 


Vids) = -, (30) 
feos fe-BIVn-tMiVildyy,- + -dwy_s 





th 

Vi(\,) is the mean value of V; for a fixed con- 
figuration of the molecular pair i and k, while 
Vi(A;) is the mean value of the same function 
taken over all configurations of all molecules of 


‘Zernicke and Prins, Zeits. f. Physik 41, 184 (1927); 
Debye and Mencke, Physik. Zeits. 31, 797 (1930). 
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the system. By Eq. (10), we may write 


tk 


(31) 


—_— 
Vi(A;) = Viet de Vis(Ai), 
Zk 


since the term V;,. is constant for fixed configura- 
tions of molecules 7 and & and is unaffected by a 
mean value operation in which these are held 
fixed. Because V;, depends only on the relative 
coordinates of 7 and s, we have 


"Vesti =(1/o.t) f ViseP® du, (a) 
(32) 


of: . - fe-Blvn IthiVildgy, o- -dwy_3 


Ee BWiks(Qi) = " . 
fees fe-BIVN-1F MV id, + « 


where the integration in the numerator of (32b) 
is to extend over the coordinates of all molecules 
except 7, k and s and in the denominator over 
all but those of 7 and k. It will be seen at once 
that W;,.*(A;) is the potential of the mean force 
acting upon molecule s, averaged over the con- 
figurations of all molecules except i, k and s. 
We note that W;,*(\;) is not invariant under a 
permutation of upper and lower indices, although 
W;*(A;) is. It is possible to define a potential 
of mean force which is symmetric in the coor- 
dinates of 7, k and s, differing from W;,*(d;) by 
a term independent of the coordinates of s, but 
the latter form is not the most convenient for 
the present purposes. The second term of Eq. 
~ 


(29) is the sum >> V;,(A;), each member of which 
s=1 

is given by an expression of the form (21). 

Moreover, each of the mean values “V;, and Vi; 

are identical molecules of the same type s, so 

that we may write, 


OW; *(A;)/Od i= Vix 


7; eT BW ik Od —eBWi*Qd |dw,, 


(33) 


where the sum is extended over all types s. Each 
of the integrals refers to a single group of three 
molecules, 7, k, and a single molecule of type s. 
Properly a term, (v0,)—f VigenB¥* it Oday, 
where k’ is a single molecule of the same type 
as k, should be omitted from the sum in Eq. (33). 
Since, however, this term bears a ratio of the 
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order of 1/N, to terms retained, it can be ignored 
altogether. This is not true of I’;,, which is of 
dominant importance when & is situated in the 
immediate neighborhood of 7. 

Before integrating Eq. (33) we shall find it 
convenient to investigate the behavior of 
W,*(\,;) in the neighborhood of \;=0. First we 
note that 1’;*(0) must be independent of the 
coordinates of 7, since the mean force on 7 
vanishes when \; is equal to zero. Moreover, it 
must also be independent of the coordinates of 
k except in the surface region, since the mean 
force on k must also vanish when \; is equal to 
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must be a constant. Making use of Eq. (20b), 
we write 


oo," f ° fe BVN dq, ee -dwy 9 


e~ OW FO) a ————— : (34) 
f- ‘ -fe BVN idq,-++dwy 
Now the integral in the denominator of (34) is 
obviously independent of the coordinate of 7 and 
k, since it includes integrations over these coor- 
dinates. Since W,*(0) is independent of the 
coordinates of 7 and k, we must conclude that 
the integral in the numerator is also independent 
of them. The factor ojo,v? may therefore be 
written as an integration over the configuration 



















zero. From these facts we conclude that W;*(0) spaces of molecules 7 and k, 


vioie* | e fe BVN dw): . duy-2= [- . fe BVN Ww: . ‘dwy. 


From Egs. (34) and (35) we conclude that e~®”"* has the value unity and that W;*4(0) vanishes. 
Bearing this fact in mind, we integrate Eq. (33) between zero and }; with the result : 





we 
“sa 


v di 
Wi (Ai) =AVine tL (N-/2) (1, of [ Vilererato0 e200 Wo drs (36 
s=1 0 
From Eq. (36), we observe that W;*(\;) approaches \;V;, as the densities of the several components 
of the fluid approach zero. This approximation is also obtained by neglecting the influence of the 
fixed molecule & on the distribution of a third molecule s in the vicinity of molecule 7, that is by 
ignoring the difference between W;,.*(A;) and W;*(\;). In the theory of the equation of state of a 
rarefied gas, this is essentially the approximation which is made. To obtain a more accurate expression 
for W;*(\;) we must know something about the potential W;,.*(A;) involving the group of three mole- 
cules i, k and s. Proceeding from Eq. (32b), which defines W;,*(\;), we obtain by differentiation 
with respect to A; 


OW i.8(\;) A:=AVietDV(M tor) [ Vale-#¥m'00 —¢-#a'0Yla, (37 
l=1 


ows: - - fe-BlVn i+hiVildg,- . -dwy_s 


Sans fe PVN AP Vida,- + -dwy_s 


e- BW iks!(Xi) 





Examination of (32b) shows at once that W;,*(0) is merely W;,* in a system containing N,-! 
molecules of type 7. It is apparent that this function can differ from W,°* in the original system on!) 
by a quantity of completely negligible order, O(1/N;). Thus when Eq. (37) is integrated with 
respect to \;, there is obtained 


y hi 
a ciaditiiiitiae  ute on | J V isle BW ito!) — 9 BWik!D den dd. 38 
!=1 0 


This equation is entirely analogous to Eq. (36). Since W,* is given by an equation of the form ©! 
(36), a first approximation to W;,*(A;) is evidently V.,+A;:V; . Similarly a first approximation © 
W;*(\i) is A:Vis. If these approximations are employed in the integrals of Eq. (36), we obtain 4 
second approximation to W;*(\;), namely 


WF (A,) =i Vi,—RT>D(N,/2) (1 (on) [ (e BriVis— 1) (e-8V ee — 1)dw,. 30 


s=l 





y 
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We should expect this approximation to be adequate in gases of moderate density. The method of 
successive approximation employed to obtain Eq. (39), can of course be extended without great 
difficulty, by transforming the potentials W;,,{(\,), involving groups of four molecules, in exactly 
the same manner as W;*(A;) and W;,.*(4,). Presumably any desired degree of approximation could be 
attained in this fashion, but the expressions become very cumbersome. It should also be remarked 
that the method of approximation just outlined can only be employed when no long range forces, 
such as the Coulomb forces, act between the molecules, for in such cases the integral of Eq. (39) 
diverges. It converges whenever the individual potentials decrease asymptotically at least as rapidly 
as ru‘. 

For an exact calculation of W;*(A;), it would be necessary to solve a set of N—1 simultaneous 
integral equations of the type of Eqs. (36) and (38), involving potentials W*;,...,, in which the 
number of fixed molecules 7, k, «--/, ranges from one to N—1. Moreover, when the number of fixed 
molecules becomes comparable with the total number of molecules, NV, the equations become more 
complicated than (36) and (38), due to growth of terms, which are of zero order in the latter equa- 
tions. Since this is not a practicable procedure, we must resort to some means of approximation. 
Eq. (39) represents a possible method of approximation in which W;*(\;) is developed in powers of 
the densities of the components of the fluid. However, we should like to find an approximation more 
suitable in condensed systems than Eq. (39) is likely to be. An obvious and simple approximation 
may be obtained by assuming superposition in the potentials of mean force, 


Win"(di) = Wi? + Wi"(a,). (40) 


Although the additivity expressed in Eq. (40) holds for the direct action of i and k on s through the 
terms V;.,+A; Vis, it cannot hold exactly for their indirect action upon s through their effect on the 
distribution of the remaining N—3 molecules. The magnitude of this discrepancy will be discussed 
later. If we introduce the approximation (40) into Eq. (36), we obtain the following integral equation: 


y hi 
WA) =A VietDTV(N,/v) 1 fo.) ) J Vaeomoroole BWW — 1 dea, di. (41) 
s=l Jo 


The potentials W;*(A,;) will have the same form for any pair of molecules of the same type. From the 
v types of molecules present in the fluid, v(v+1)/2 types of molecular pairs can be formed. There 
will, therefore, be v(v+1)/2 types of functions W;*(\,;), each satisfying an equation like (41). Simul- 
taneous solution of these integral equations will yield values of the functions W;*(\;) consistent 
with the approximation expressed in Eq. (40). Of course for a pure fluid, having only one com- 
ponent, only one integral equation has to be solved. 


N 
W(qi2, %) =AV (Giz) + (No, Mo) { J Veqnieramcos Te F¥ (a3, DV — | \dwadh, (42) 
0 


where 12, g13 and go3 are the relative coordinates of a set of three molecules, N is Avogadro’s number, 
\f the molecular weight of the substance, and p its density. We note that the approximation (39) 
is obtainable from Eq. (41) as the first step in its solution by successive substitution. In general, it 
will not be feasible to seek analytic solutions of Eqs. (41) or (42). However, by graphical or mechan- 
ical quadrature, it should be possible to determine self-consistent functions satisfying these equa- 
tions. We remark that this process provides a theoretical means of calculating the x-ray scattering 
curves of liquids, since they are determined by W;*(1). In some cases the potentials of intermolecular 
force can be calculated from atomic structure. Although this can be done in relatively few cases, 
quantum mechanics prescribes an approximate functional form for the V;.. With this functional form, 
the numerical values of V;, may be deduced from the temperature dependence of the second virial 
coefficient of the substance in the gaseous state. 

While Eq. (40) represents an approximation which is intuitively attractive, it is of course not 
exact. Although a numerical estimate of the error involved is difficult, we can say something quite 
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definite about the nature of the approximation. To facilitate discussion, we shall introduce the 
symbol ¢;,*° to denote the deviation from additivity 


€ix° = Wi8(Ai) —L Wi + Wii) J. (43) 


It is possible to express €;;* in terms of a third order fluctuation in the potential energy of molecule 7. 
Without entering into the details of the calculation, we state the following result. 


Ms pl pi : 
1h s—— iks—— bihainnne 
cut= ef ff (Vi-._—s Vi)(Vie-—— (Vie) (Vs— Vs) ddd di, (44) 
0 Yo Yo 


where the fluctuation is understood to be a function of the variables \;, Ax, As, mean values being 
taken with the distribution function e~8”N+: «>». A more useful expression for €;,* is obtained 
simply by subtracting the equation for W;*(\;) corresponding to (36) from Eq. (38). The result is, 





= 


hi 
ee ee _aw.l ee _aw..lrn: ” 
é.2= (Ni/eai) f fi ap [eT BW tks i) 4 e—BWi (AD) — @-BWik!(D — e-8Wis' OD) dey dj. (45) 
1 0 


l 


ll 


An exact calculation of ¢;;.* from Eq. (45) is of course impossible without a knowledge of the potential 
Wix-<'(\;). However, we should expect to obtain at least a rough estimate of €;,° by using the addi- 
tivity approximation for the potentials W;;..'(A;), Wi.'(A;), and W;,'(\;). Introduction of this approx- 
imation yields 


v hi 
€ix* = 2) (N 1/201) f f Vise 8 0 (eB! — 1) (eB s! — 1) dad di. (46) 
l=1 0 


If we rely on Eq. (46) for the order of magnitude of €;,*, we see that the integrand will be very small 
except for configurations of molecule such that W;' and W,' are simultaneously large. If s and & 
are situated at some distance from each other so that their spheres of influence on the distribution 
of a third molecule /, do not overlap appreciably, ¢;,° will be small. Returning to Eq. (36) we remark 
that for the calculation of W;*(A;), it is only necessary to know W;,*(A;) accurately in a region in the 
vicinity Of molecule 7, where V;,; is sensibly different from zero. If the molecules 7 and k are separated 
by an appreciable distance, s will be distant from k in the greater part of this region, and we may 
conclude that ¢;,* is small. Thus we may expect Eq. (41) to provide a good approximation to 
W;*(A;) at least when the molecules 7 and k are separated by an appreciable distance, and that 
W;*(\;) will approach the solution of Eq. (41) asymptotically as the distance between i and k is 
increased. This point is of some importance in the case of long range intermolecular forces, such as 
those between electrolyte ions. At small distances between 7 and k the dominant term in W;*(\;) 
becomes V;,, and although the terms depending on the density may still be large, it is not necessary 
to know them with such great accuracy. Once solutions of Eq. (41) have been obtained, they may 
be used in Eq. (46) to obtain a rough estimate of ¢;,*. If this turns out to be small, the solutions of 
Eq. (41) can be accepted with some degree of assurance as reasonable approximations to the poten- 
tials W;*(\;). It would be desirable to have a more rigid criterion for judging the adequacy of these 
approximations, but unfortunately this cannot be obtained without introducing distribution functions 
involving larger and larger groups of fixed molecules. Unless some approximation such as super- 
position is introduced at some early stage in this sequence of distribution functions, we are driven 
back to the many dimensional phase integrals with which we started. 

Having obtained approximate expressions for W;*(\;) by solution of Eqs. (41), we can calculate 
the chemical potential of any component of the fluid by means of Eq. (22), as well as the total 
energy of the fluid by Eq. (25). From the chemical potentials, we may obtain the Gibbs free energy, 
¢, of the system. 

¢=E+pv—TS=>Niwi=kT Ni log Ni/v+ x 
=1 


i=1 i, k=1 


NN; al v a 
| Bunda} +E Niei(T) (47) 
v do iml 
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In principle one can obtain the compressibility, «x, of the fluid by means of the thermodynamic 
formula, 

1/x= —v(dp/dv) p= — (06/00) r, ny, “yy (48) 
However the B;,(A;) depend upon v through the W;*(\;). Even though Eq. (41) may provide good 
approximations to the W;*(A;), it is by no means certain that the derivatives of these solutions with 
respect to v will be good approximations to dW;*(A;)/dv. Thus while Eqs. (22), (35) and (47) may 
express the dependence of the u;, E and ¢ upon the composition of the fluid with some accuracy, 
the compressibility calculated from Eqs. (47) and (48) may be a much poorer approximation. Quite 
formally, we may write 


ap ie N; » NN; f' OB (Xj) 
( ) 6M 4 fF | Bu) —0 at. (49) 
Ov T | i=1 vy? 1, k=1 v' 0 Ov 


In case the fluid under consideration is a gas, we know that as the volume is increased to infinity, 
the pressure decreases continuously to zero. In this case we may integrate Eq. (49) between the 
limits v and infinity. After a simple partial integration, the following equation of state is obtained. 


po > Xe i ° Bi-(X3) 
— =] a > —f M Bald + v? { avin (50) 
eh i, k=1 v 0 e p3 j 

n;=N;/N; xi=ni/>D ny; d=v/)> ny, 


@o 


where N is Avogadro’s number, 7; the number of moles of component 7 present in the gas and x; 
the mole fraction of that component. Eq. (50) may be made the basis for a statistical calculation of 
the equation of state of a real gas mixture. When the potentials W;*(A;) in the B;,(A;) are approx- 
imated simply by \;V;, the usual expressions for the second virial coefficients are obtained. 


IV. 


The present method provides a particularly simple means of determining the equation of state 
of a real gas mixture. We shall derive expressions for the second and third coefficients in the virial 
expansion of the equation of state, indicating how the method may be extended to higher coefficients. 
Results are obtained which are in agreement with those of Ursell, who has treated the same problem 
from a somewhat different point of view.® 

If the function W,*(A;) is approximated by Eq. (39), we obtain the following expression for B;).(\;). 


1 v N,z 
Bii(di) = [ve exp —BVit>d 


o.RT s=1 Us 


fa —eBriVis) (1 —e# day) dan. (51) 


Substitution of this expression in Eqs. (22) and (50) would yield the chemical potentials of the 
several components and the equation of state of the gas mixture. We should expect the equation of 
state, although approximate, to be much more satisfactory than a virial expansion. However, the 
integrals involved cannot be expressed in terms of elementary functions, and we shall not concern 
ourselves with their simplification here. We note that Eq. (39) is exact to terms in 1/v?, so that the 
first two terms of an expansion of the right-hand side of Eq. (51) are exact. The expansion is 


1/1 


kT 


Ck 
» a | . 
+> —— f f V ipe-BMVik(1 —@-BMV ie) (1 —e-8V denne, t+O(1/v8) $. (52) 
s=1 VU O30; 


Fae Proc. Camb. Phil. Soc. 23, 685 (1927); R. H. Fowler, Statistical Mechanics, Cambridge University Press (1929), 
p. 1/5, 
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Substitution of this expansion in Eq. (22), integration with respect to \; and a rearrangement of 


terms yields® 


v N. 
w= kT) log Ni/v +E — | 1-76") de, 


k=1 Vo k 


+ NN, | 
+32 — f fa-e BVik)({ —e- BV is) (1 —e 8Y's)\dwj,dws,+O(1/v*) >. (53) 
0:0; | 


k, s=1 ve 
When the approximate expressions for the y; given by Eq. (53) are introduced into Eqs. (47) and 
(48), we obtain 


ap »N; +» NiNz 1 
T v ok 


ov i=1 gy i, k=1 


NiN-iN; | 
‘*i me — f fa-e BVik) (1 —e~8V is) (1 —e- PY 's)\dwdw,+O(1, [P)¢. 


i, k, s=1 vt OO s 
Integration between the volume limits ~ and v leads to the following equation of state 


pi » Mike yr MGX 
—=1+ > —B;'+ F a Cine + O(1/89) 


RT i, k=1 DP i,k, s=1 9p? 


N 
By®=— | (1-78 *)da, 


2e% 


Cit — f fae BVik) (1 —e-BV is) (1 —e- FY *s)\dwipdas. 


301.0 
If the method of approximating W;*(A;) embodied in Eq. (39) is carried one stage farther, one more 
term may be retained in the expansion of B;,(A;) and the fourth virial coefficient may be calculated. 
The process can be continued indefinitely, but of course the algebraic operations become tedious, 
and the final expressions unwieldy. The method consists simply of continued substitution and 


* After integration with respect to \;, the terms in yu; involving 1/2? have the form 


N 
Viz —8( Vint; r. 
—e BV ks B(ViktVie) _ 4) _ (o-BVik 
2 Fan SS ' ‘ a Viet Vio (e 1) (e 1) } dw,.dws, 


k,s=1 

k+s+1 
where the original summation indices, referring to the individual molecules, have been restored for clearness. For every 
term in the sum with the integrand 


(e BWV iktVis) _ 1) — (ge BVik_ —p BV ks 
ea a (e »}a eBVke), 


corresponding to any specified pair of molecules & and s, there is also a term involving the same pair, having the form 


e BV iRtVis) _ = —BVis__ —_ —BV ks) 
{ote 1) (PM —1) | Le PY), 


The sum of this pair of terms is readily seen to be 
(1 —e BV ik) (1 —e BV is) (1 —e BV ks) 


so that the double sum may be rearranged as follows: 
1s a ff he BMH) 1 0AM) (1 PY iB) dernds 
ee VoLos —_ 
k+s+$t 
This sum contains (N—1)(N—2) terms of similar magnitude. It differs from the sum over molecular types appearing i 
Eq. (55) by (3N—2) terms of the same magnitude. The aggregate of these terms bears a ratio of the order of 1, \ to 
either sum and is consequently entirely negligible. 
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expansion. If one wishes a virial expansion including terms of the order of 1/v", one begins by 


n—1 


approximating the W*;,...,, involving m fixed molecules, by sums of the form A; V;,+20 Vis. These 
l=] 


expressions, substituted in an equation similar to (38), yield an approximate W’;,..., involving (n—1) 
fixed molecules. Continued substitution finally leads to an expression for W;*(\;), which when ex- 
panded in powers of 1/v, is exact to terms involving 1/v". It is only necessary to have the expression 
for W*;,.... involving (7—1) molecules, exact to terms in 1/7, since in the final expansion, it is mul- 


tiplied by a factor 1/v"~'. 


V. 


It is interesting to investigate the form of the 
molecular pair distribution function, e~®”' by 
means of the approximation, Eq. (39). Although 
Eq. (39) can scarcely be accurate at liquid 
densities, we might expect it to be adequate up 
to moderately high densities. Moreover, even 
where inexact, it should serve to illustrate the 
form of the distribution function. 

We shall restrict our attention to a pure fluid, 
consisting of spherically symmetrical molecules. 
For such a fluid, Eq. (39) reduces to the fol- 
lowing expression when \; is equal to unity 


W(r) = V(r) —(NRT/v) 


x fc —e~BV(ris)) (4 — e~ BV (725) \dys, 


r= r—fo| ; rig= |T1—Ts3 > Yas= |¥e—T 3}. (56) 


Since all molecules are of the same type, W;* 
and |’; have the same functional form for every 
molecular pair, and we may omit the indices 7 
and k. We shall suppose the molecules to be 
rigid spheres of diameter }b, so that when r is 
less than b, V(r) becomes positively infinite. 
Attractive forces outside of the sphere a should 
properly be represented by a potential of the 
form, —y/r®. However to simplify calculations, 
we shall represent this potential by a trough of 
depth wu» and width, a—d. Between the con- 
centric spheres of radii 6 and a, V(r) is thus 
equal to —uo. The potential V(r), upon which 
the following calculations are based is repre- 
sented by the following step function. 
V(r) = @, O<r<cb 

=-—u, b&rga (57) 

=0, a<r. 


With this potential, the integral on the left-hand 
side of Eq. (56) may be expressed in terms of the 


volumes of intersection of spheres of radii a and 


b. 
Wr) /kT = Vi(r)/RT —(Noo/v)J(r), 
Wy = 4rb3/ 3 ‘ 


1 
cd } wip t+ 2( war —_ wp) (1 — eFuo) 
Wh 


+ (Waa — 24+ wr») (1 —eFuo)?} 


J(r)= 


where w),, is volume of intersection of two spheres 
of equal radii b, w,, that of two spheres of equal 
radii a, and w,, that of two spheres of unequal 
radii a and b, the distance between the centers 
of each pair of spheres being r. A simple cal- 
culation yields 

Whb 

—=1-—37r/4b+ 7¢(r/d)'; 


Wo 


O0<r<¢2b 


Waa 
— =(a/b)*(1 —3r/4a+ 6(r/a)*]; 


Wo 


O<r<2a 


0, r>2a 
(59) 
Pif1—3(1—(a/r?))(r/b) 
+ 6(1 — (a/r*)3(r./b)?]+(a/b)* 
x (1-$(1+(a/r?))(r/a) 
+ 4(1+(a/r?))3(r/a)*]} ; 
(a?—b?)'<rgat+b 
a=a’>—}*. 


=0, r>at+d; 


By means of Eqs. (58) and (59), W(r)/kT may 
be computed as a function of the ratios, 7/bd, 
a/b, Nwo/v and uo/kT. The probability dis- 
tribution function, which we shall designate by 
G(r) may then be calculated from the relation 
G(r)=e—-¥/kT, We note that Nw is equal to 
eight times the total volume of the molecules, 
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Fic. 1. Molecular pair distribution function based upon 
Eq. (58). 


their radii being 6/2. For most liquids the ratio 
Nw,/v would have a value between 4 and 6. 
Moreover, we should expect u»/kT to be at 
least unity in liquids. The ratio a/b is somewhat 
arbitrary. However, in order to make the rec- 
tangular trough approximate the actual poten- 
tial of the attractive forces, which diminishes 
rapidly with increasing distance, a/b cannot 
differ much from unity. In order to approach 
conditions in the liquid state, we assign the 
following values to the constants in question: 
Nwo/v=2; uo/RT=1; a/b=1.25. The above 
value of Nwo/v corresponds to a density about 
a third of ordinary liquid densities. With the 
above values of the constants, G(r) has been 
calculated with the aid of Eqs. (58) and (59). 
It is plotted as a function of r/b in Fig. 1, and 
tabulated in Table I. Actually there is a discon- 


TABLE I. 








G(r) | r/b 


0.00 
0.00 
2.48 
2.23 
0.81 


G(r) 


0.86 
0.98 
1.20 
1.25 
1.03 
1.00 














tinuity in the curve at r/b=1.25, arising from 
the artificial nature of the step function poten- 
tial. In the figure, this discontinuity has been 
eliminated by drawing a smooth curve through 
points calculated at r/b intervals of 0.2. It will 
be observed that the calculated G(r) bears a 
close resemblance to the molecular distribution 
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functions in liquids derived from measurements 
of x-ray scattering, exhibiting the characteristic 
maxima and minima of the empirically deter- 
mined functions. While we realize that at the 
temperature and density of the above calcula- 
tion, the approximation (39) is being pushed to 
the limit of its applicability, we believe that it 
adequately illustrates the mechanism by which 
the intermolecular forces operate to establish 
molecular pair distribution functions in dense 
fluids. 

It is interesting to observe that the attractive 
van der Waals forces play a dominant role in the 
production of the secondary maximum and 
minimum in the G(r) of Fig.'1. On the basis of 
Eq. (39), the distribution function in a fluid of 
rigid spheres with no attractive forces would 
exhibit a single maximum at r=b. The absence 
of secondary maxima and minima for rigid 
spheres is of course the fault of Eq. (39). If the 
approximation were carried one stage farther, 
they would begin to appear. However, the point 
to be emphasized is that probability distribution 
functions derived from the packing of rigid 
spheres can scarcely be applied to actual fluids, 
where the attractive forces play a dominant role 
in determining the form of the distribution 
function. 

A further limitation of the approximation (39) 
is that it fails to produce more than a single 
pair of secondary maxima and minima. For a 
more exact calculation, it is suggested that the 
distribution function derived from Eg. (39) be 
used as a starting point for the self-consistent 
solution of the more exact integral Eqs. (41) and 
(42). Thus, it will be observed that tertiary 
maxima and minima would doubtless arise from 
the secondary maximum and minimum of the 
approximate W(r, \) based on Eq. (39), if the 
latter function were substituted in the integral 
of Eq. (42). The computational difficulties in- 
volved in the self-consistent solution of Eq. (42) 
are unfortunately rather serious, since even for 
spherically symmetrical molecules, a_ triple 
numerical or mechanical integration is required 
at each stage of the process. However, these dif- 
ficulties are not insurmountable. 

The formulation of a theory of solutions on the 
basis of the present ideas offers intriguing p0s- 
sibilities. Several methods of approach are open, 
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the direct use of our expressions for the chemical 
potentials being the most obvious. An alternative 
and more phenomenological treatment applicable 
to dilute solutions has also been undertaken. 
These investigations will be described in a future 
article. 


Note added in proof: 

(a) A remark about the normalization of the molecular 
pair distribution functions is perhaps worth while. In 
defining the potential of mean force, W;*(A;), ar arbitrary 
constant is at our disposal. This constant has been so 
chosen that 


(1/von)feP¥* inde, =1 


a relation which may be verified by reference to Eq. (20b). 
The approximate forms for W;*(\;) of course do not exactly 
satisfy this condition. 

(b) At sufficiently high temperatures, the exponentials 
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of Eq. (41) may be expanded and the following approxi- 
mate linear integral equation is obtained 

hi N, 1 

W#(\) =i Via -— Do —— | VisWi'(1)dor, 

kT Ss Vv od;* 
where regions in which W;* and W,* become infinite due 
to repulsive forces, are excluded from the integration. If 
we apply the equation to the ions of an electrolyte solution, 
the solvent being treated as a continuous medium of 
dielectric constant D, we know that Vix is equal e:e./Drix 
where e; and e are the ionic charges. For this potential, 
it may be verified without difficulty that the above set of 
linear integral equations has the following solutions 


€:€i- 4r N, 
Ws) =r eres t= as 
ss « DT hs 


if molecular size is neglected. These solutions, together 
with Eq. (22) lead to the Debye-Hiickel limiting law for 
the chemical potentials of the ions. 
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Some Critical Remarks Concerning Bailli Nilsen’s Paper: ‘On the Theory of Free 
Radicals and Organo-Alkali Compounds' ”’ * 


Eceit A. HyLterAas, Chr. Michelsens Institutt, Bergen, Norway 
(Received February 27, 1935) 


ATELY, a paper by Bailli Nilsen on the 
above title appeared in this journal, the 
main subject of the article being a computation 
of the electron affinities of free aromatic radicals. 
The author first mentions two possible ways 
of a quantum-mechanical calculation of the 
electron affinities of such free radicals,—the 
“first method” certainly alluding to some ordi- 
nary first order perturbation method for an atom, 
or, a molecule when adding-an extra electron on 
an unoccupied orbit—whereas the ‘‘second 
method,”’ as he says, is closely connected with 
Pauling and Wheland’s theory for free radicals. 
If now, for instance, we apply such an ordinary 
perturbation method to the simplest of all atoms, 
the hydrogen atom, we will find a very poor ap- 
proximate value of the electron affinity—even 
the sign being wrong—viz., the negative value 
—0.25 Rh or —3.385 electron volts—whereas an 
exact calculation provides a positive value 0.0506 
Rh or 0.685 electron volt. 


i. Chem. Phys. 3, 15 (1935). 
Publikasjoner fra Chr. Michelsens Institutt, Nr. 50. 


It is true, therefore, that there are many diffi- 
culties in determining the energy difference 
caused by putting an extra electron into an atom 
or a molecule. We also know that these difficulties 
are closely connected with considerable changes 
of the electronic configuration of the atom, be- 
cause the electron orbits are deformed by the 
mutual interaction of the electrons. These 
changes can only be accounted for by the aid of 
complicated variational methods. 

The author of the above paper now seems to be 
of the opinion that such difficulties may be con- 
quered simply by dropping them out of sight and 
by making use of the essentially different ‘‘sec- 
ond method”’ of Pauling and Wheland, which is 
especially suited for quite other purposes, v7z., 
the computation of binding forces between atoms 
of given electronic configurations. 

We will, however, quite overlook the consider- 
able errors introduced in this way by neglecting 
energy changes in the individual atoms them- 
selves, and only analyze the author’s application 
of the “second method.” 





314 BGEIL A. 

In Pauling and Wheland’s paper electron or- 
bits in individual atoms are designated by letters 
a, b, c, «++ and exchange integrals between such 
orbits by ab, ac, bc, etc. The exchange integrals 
are, therefore, mathematically well-defined quan- 
tities depending only upon the geometrical 
properties of the supposed electron orbits, which 
may, for instance, be expressed by some definite 
kind of eigenfunctions. The values of the integrals 
have nothing at all to do with the binding elec- 
trons occupying those orbits, they only enter into 
the secular equation of the problem on account 
of the presence of such electrons. 

Now Pauling and Wheland in dealing with the 
seven-electron problem of the benzyl radical, 
formally add an eighth electron, called ‘‘phantom 
electron,’ and a corresponding ‘phantom orbit” 
h. Equating to zero all exchange integrals be- 
tween real orbits and the ‘phantom orbit’”’ they 
get a number of singlet terms corresponding 
exactly to the number of doublet terms of the real 
problem and with the same energy values. 

This procedure may be considered as a purely 
formal trick, but we may also imagine that it has 
a definite physical meaning. The “phantom 
orbit’ may well be interpreted as a real electron 
orbit in a second molecule at some distance 
which allows the exchange integrals mentioned to 
be put approximately equal to zero, thus giving 
qualitatively a picture of a faint, although real 
and active valence force pointing from the free 
radical under consideration to some second dis- 
tant atom or molecule. 

In Bailli Nilsen’s paper, however, there is no 
clear distinction between the binding electrons 
themselves and the possible orbits of the elec- 
trons. Sometimes he speaks of exchange integrals 
between orbits, sometimes between electrons, 
and he never tries to give a clear definition of the 
“phantom orbit” introduced. On the contrary he 
only speaks of the “‘partner’’ of the “charging 
electron” in that or that “canonical structure.” 
Following Pauling and Wheland he assigns defi- 
nite values to the exchange integrals ab, ac, ---, 
bc, bd ---, between the 7 original orbits, a, --- g, 
either the common approximate value a or the 
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approximate value zero according to whether the 
orbits belong to neighboring atoms or not. In 
this case he quite correctly takes no account of 
the direction of the valence forces in various 


“canonical structures.”” In the case of the 


‘phantom orbit,’’ however, he puts the exchange 
integrals ah, bh, ch, etc., sometimes equal to a 
and sometimes equal to zero—in the very same 
secular equation. 

Thus in his matrix element H;; we have to add 
integrals of this kind to the complete expression 


— }(ah+bh+ch+dh+eh+fh) + gh 


and correspondingly in his matrix element H;; 
we get the expression 


s[ah+ch—2 bh—3(dh+eh+fh+gh) }. 


In the first matrix element H;; he now puts 
gh=a, all other integrals equal to zero, whereas 
in the second element Hs; he assigns the value a 
to ah and ch now putting all remaining integrals, 
including the above gh, equal to zero. 

Having completely overlooked the necessity of 
a “materialization”’ of the ‘‘phantom orbit’ he 
can freely fancy that his “phantom electron” is 
sometimes here and sometimes there, in the 
“‘neighborhood”’ of this or that ‘‘partner.”’ 

If it had originally been the author’s opinion 
that the ‘‘charging electron” is in some way 
continually jumping from one atom of the mole- 
cule to another, he should, or course, have made 
use of the available electron orbits, one per atom. 
Then there would be a different number of orbits 
and electrons, and the special formal method of 
“superposition patterns’ put forth by Pauling 
and Wheland could not have been applied at all in 
the given form. Also in this case, of course, the 
difficulties of determining the perturbation 
energy of the individual atoms remain to the 
same degree. 

The conclusion I want to give is, therefore, 
that in my opinion the whole treatment of the 
electron affinities of free aromatic radicals as 
given by Bailli Nilsen is unclear and misleading— 
and if not thoroughly revised—it must be con- 
sidered as meaningless. 
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Remarks on the Theory of Aromatic Free Radicals 


Linus PAULING AND G. W, WHELAND, Gates Chemical Laboratory, California Institute of Technology 
(Received March 5, 1935) 


N the quantum-mechanical treatment of the 

benzene molecule (and of other aromatic 
molecules)! the approximation is usually made of 
neglecting the pairs of electrons which are con- 
sidered to form bonds in the plane of the ring. 
The problem then becomes one of discussing the 
wave function which can be constructed for the 
remaining electrons by use of the p, orbitals, one 
for each carbon atom. This discussion can be 
made either by the molecular orbital method or 
by the valence bond method, the results being in 
general agreement. The treatment of a free 
radical such as phenylmethy] is closely similar, 
seven p, orbitals and seven electrons being con- 
sidered for this molecule. It is found convenient 
in applying the valence bond method? to intro- 
duce a “phantom orbital’’ and an additional 
electron, a procedure which simplifies the treat- 
ment without changing the results. 

The quantum-mechanical treatment of the 
electron affinity of a free radical involves the 
similar treatment of the negative ion, the wave 
function for the ion being formed by the use of 
the same p, orbitals as for the radical, but with 
one more electron (eight electrons and seven p. 
orbitals for the phenylmethyl negative ion). The 
discussion of this problem by the molecular orbi- 
tal method has been given by Wheland,' who 
showed that the electron affinity is approximately 
equal to the difference in coulomb energy of the 
ion and radical, the difference in resonance en- 
ergy being small in all the free radicals. This 
result is not in disagreement with experiment. 
The valence bond method can also be applied to 
this problem. For example, we might consider for 
the phenylmethyl radical (as an approximation) 


'E. Huckel, Zeits. f. Physik 70, 204; 72, 310 (1931); 76, 
628 (1932); 83, 632 (1933); L. Pauling and G. W. Wheland, 
J. Chem. Phys. 1, 362 (1933); L. Pauling and J. Sherman, 
ibid. 1, 679 (1933); G. W. Wheland, ibid. 2, 474 (1934). 

*L. Pauling, J. Chem. Phys. 1, 280 (1933). Cf. also L. 
Pauling and G. W. Wheland, reference 1. 


only the five unexcited structures 
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in which the dot represents the odd electron (con- 
sidered bonded to a phantom orbital in our 
graphical treatment), and for the ion the five un- 
excited structures 
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in which the pair of dots represent two electrons 
(with opposed spins) occupying the same p. 
orbital. It is found, however, that the secular 
equations for the ions involve several integrals 
which are presumably not small enough to be 
neglected and for which we have not been able to 
estimate reliable values; for this reason we have 
not carried out the calculation. It seems prob- 
able, from the form of the canonical structures of 
radical and ion, that the valence bond method 
would give essentially the same result as the 
molecular orbital method, namely, that the dif- 
ference in resonance energy of radical and ion 
is small. 

The attempted application of the valence bond 
method in the calculation of the electron affini- 
ties of free radicals made by Nilsen? is criticized 
by Hylleraas in the preceding paper.‘ We agree 
with Hylleraas that Nilsen has not correctly eval- 
uated the resonance energy, and in addition we 
point out that he has neglected to mention the 
changed coulomb energy, which in our opinion is 
the major contributor to the electron affinity. 


3B. Nilsen, J. Chem. Phys. 3, 15 (1935). 
4E. A. Hylleraas, J. Chem. Phys. 3, 313 (1935). 
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The Isotope Ratio in Petroleum 


The natural separation of isotopes of hydrogen in 
benzene and kerosene found by M. Dole! was called in 
question by the author himself because of later experiments 
with benzene (under modified conditions).? The results of 
his initial experiments were explained by him as a conse- 
quence of fractional condensation of water vapors from 
the products of combustion. 

I burned the motor-petroleum from Machach-Kala in a 
calorimetric bomb, so that the possibility of fractional 
condensation was excluded. All the water obtained was 
subjected to an elaborate purification. The interferometric 
analysis gave 1.5 parts deuterium per 5000 parts of 
hydrogen for the concentration of deuterium in petroleum 
if, as usually, the mole fraction of the former in ordinary 
distilled water is taken as 1 : 5000. If the tank oxygen 
used for the combustion has an increased concentration 
of heavy oxygen isotopes (as it was found by Smith*), 
then according to Luten‘ this would decrease the value 
of An, i.e., the ratio 1.5 : 5000 for concentration of 
deuterium in investigated motor-petroleum is decreased. 

I wish to express my sincere gratitude to Professor 
Dr. A. E. Brodski for suggesting this work and for his 
valuable advice. 

N. S. FILippova 

Institute of Chemical Technology, 

Dniepropetrovsk, U.S.S.R., 
February 25, 1935. 
1M. Dole, J. Chem. Phys. 2, 337 (1934). 
2M. Dole, J. Chem. Phys. 2, 548 (1934). 


3 E. R. Smith, J. Chem. Phys. 2, 298 (1934). 
4D. B. Luten, Phys. Rev. 45, 161 (1934). 


The Intensities of the Vibration-Rotation Bands of HCl 


The writer’s attention has recently been called to 
Professor Mulliken’s Note on Electric Moments and Infrared 
Spectra and to the subsequent note of correction.' In his 
second communication Mulliken takes cognizance of the 
recent intensity measurements of Bartholomé on the 0—1 
bands of HCl, HBr and HI, and of the serious discrepancy 
between the results obtained by Bartholomé and those 
previously obtained by Bourgin.? Mulliken accepts at its 
face value Bartholomé’s argument in favor of the superi- 
ority of his method and combines the Bartholomé result 
for the 0—1 band of HCl with that of Dunham? for the 
0—2 band to derive a pair of alternative values for the 
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coefficient of the quadratic term of the expansion of the 
electric moment p in powers of the vibrational displace- 
ment. 

Since the measurements of Bourgin and Dunham were 
both carried out by the same method under the general 
supervision of the writer, some comment on the validity of 
their conclusions is, perhaps, in order. The method consists 
in the observation of the areas A under the absorption 
curves (measured in percent) of the individual band lines 
for a variety of tube-lengths x. The initial slope of this 
curve (dA /dx)» is the integrated absorption coefficient « 
for the corresponding line. The criticism’s of the method 
made by Bartholomé are (a) that there may be some 
surface absorption due to adsorbed vapor on the end 
plates of the absorption tube, and (b) that there must be 
a very large probable error in the determination of the 
initial slope of the A, x curve in view of the marked 
curvature in the neighborhood of the origin. In the writer's 
opinion these criticisms are wholly inadequate to account 
for the discrepancy between Bourgin’s result and Bar- 
tholomé’s. 

According to Bartholomé the true intensity of the 0-1 
band should be only one-fourth of that observed by 
Bourgin. This discrepancy factor is subject to a correction 
even if we grant for the moment the validity of Bar- 
tholomé’s method and results. In his second paper Bourgin 
re-estimated the initial slope of the A, x curve by several 
methods and corrected numerical errors in the previously 
published computations of the transition probabilities and 
effective electric moment of the molecule. The most 
probable value of the absolute integrated absorption 
coefficient obtained by averaging the results of four 
methods given in the second paper is 41 X10" instead of 
the value 49X10!" given in his first paper. The corre- 
sponding value of the Einstein transition probability for 
the first line of the R branch is‘ 8.6 X10". This transition 
probability, which we denote by By:*', is to be distinguished 
from the average transition probability for the entire band 
(Bu) observed by Bartholomé. These two quantities would 
be the same according to the summation rules if the 
frequency differences between the different band lines were 
negligible. Assuming the theoretical formula for the relative 
intensities of the band lines derived by the writer and 
confirmed within the limits of experimental error by 
Bourgin, calculation indicates that at 20°C Bo, *+! =0.96Bu. 
From Bartholomé’s value of Bo we obtain (Bo:*)parth. 
= 1.8310", or approximately 21 percent of Bourgin’s 
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value. Thus the use of Bourgin’s revised figures together 
with the application of the above correction factor slightly 
increases the discrepancy previously noted. 

Let us now return to the criticisms of Bourgin’s method. 
As regards adsorbed vapor, one may note that the vapor 
pressure of HCI at room temperature is about 45 atmos- 
pheres. Since the pressure of the gas actually used was 
one atmosphere, one may reasonably assume that the 
adsorbed layer would not be much more than one molecule 
deep. Taking the cross section of the molecule standing on 
end to be 107'* sq. cm, one readily calculates that the 
number of adsorbed molecules per square centimeter on 
the front and back quartz surfaces is equal to the 
number per sq. cm in a layer of gas less than a hundredth 
of a millimeter thick. Even if the adsorbed molecules were 
free to rotate normally and contribute as much to the 
absorption as free molecules, their contribution would 
have been too small to detect experimentally. Moreover, 
the base-line readings were taken without removing the 
adsorbed gas by heating, so that most of its contribution 
to the absorption must have been eliminated by sub- 
traction. 

Finally an examination of the A, x curves (cf. Fig. 4 in 
the first paper and Fig. 1 of the second paper) should make 
it evident that Bourgin’s estimate that the extrapolated 
value of the initial slope is reliable within 20 percent, is a 
fair one. These estimates were made with the aid of two- 
parameter semitheoretical curves which were fitted to the 
complete absorption curve and gave due weight to the 
observations made with long tubes. In fact we can derive 
from the A, x curves a lower bound to the value of Bo*! 
which is several times as large as that given by Bartholomé. 
From the fundamental formula for the area, viz., A(x) 
= {(1—e-**)dv, it follows that dA /dx = fu(v)e**dv. Hence 
dA /dx and all the higher derivatives of A(x) are mono- 
tonically decreasing functions of x. It follows that the 
initial slope must be greater than the slope of any chord 
connecting two points on the curve. Joining the origin with 
the points for the two shortest tubes and using the data 
for the third line of the R branch (the strongest in the band) 
one obtains the lower bounds 7.0510" and 5.8410" 
for By **. The other lines yield similar bounds. Thus it is 
quite impossible to reconcile Bourgin’s data with Bar- 
tholomé’s low value of the transition probability. 

Let us next consider the possibility of error in Bar- 
tholomé’s experiment. The observations consisted in 
measuring the absorption with a relatively wide spec- 
trometer slit and an absorption tube containing one 
atmosphere of HCI mixed with air at total pressures up to 
60 atmospheres. The theory of the experiment is that at 
sufficiently high pressures the lines will be broadened 
enough to destroy the rotational structure of the band and 
give a smooth curve for the absorption coefficient plotted 
against frequency. Under these circumstances with a slit- 
width small enough to cover one line at a time, but not 
necessarily small enough to resolve rotational structure, 
the radiation reaching the radiometer should be sensibly 
homogeneous as regards its absorption coefficient. By 
means of the exponential absorption law it should then be 
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possible to determine the smeared-out absorption coeffi- 
cient curve from observations made with a single tube- 
length. It is assumed that the transition probability is 
sensibly independent of the pressure and on that point 
the evidence is reasonably conclusive. 

It is rather difficult and dangerous to criticize an 
experiment from as meager an account of the experimental 
details and data as that given by Bartholomé, but one 
weak point isapparent. He gives his slit width as equiva- 
lent to 800A. If we assume that the entrance and exit 
slits of the spectrometer were of equal width, as is cus- 
tomary in such work, the extreme range of wavelengths 
reaching the thermopile would be 1600A. As the distance 
between the points of maximum intensity for the R and P 
branches is about 1400A, it is quite clear that even when 
the lines were broadened so as to eliminate the rotational 
fine structure, the radiation reaching the microradiometer 
must have been very inhomogeneous as regards its ab- 
sorption coefficient. In the R branch for example, the 4th 
line has an intensity five times as great as the 8th line 
which differs in wavelength by only 750A. Thus it is not 
permissible to extrapolate to zero tube-length by means 
of the exponential absorption law and an attempt to 
determine the transition probability by means of such an 
extrapolation is sure to underestimate seriously its true 
value. With higher resolving power Bartholomé’s method 
should be good. 

Epwin C. KEMBLE 

Harvard University, 

Research Laboratory of Physics, 
March 27, 1935. 
1R. S. Mulliken, J. Chem. Phys. 2, 400, 712 (1934). 
2D. G. Bourgin, Phys. Rev. 29, 794 (1927); 32, 237 (1928). 
3J. L. Dunham, Phys. Rev. 34, 438 (1929). 


4The value given in Bourgin’s first paper was wrong by a factor 
of two. 


Active Hydrogen 


Four years ago I expressed my doubts in what at that 
time was considered to be proofs for the existence of 
triatomic hydrogen.! It appeared that a final decision had 
been reached when Conrad in his mass-spectrographical 
work? found an uncharged molecule corresponding to the 
mass (3). Unfortunately the significance of his investigation 
was greatly reduced’ after deuterium was discovered. A 
strict distinction should be made between HD and H# 
even in mass-spectrograph analysis. This has not been 
done in every case.‘ Hence positive ray analysis as 
applied to this problem must take account of the two facts, 
(1) that the presence of HD may obliterate an unequivocal 
decision, (2) that only such particles’ can be considered as 
a proof of H* which were uncharged for some definite 
length of time and received their charge later on by 
electron transfer. 

The H* molecule may have properties identical with 
those of the H-atom. One would have to assume this, 
reading the publications of Grubb and co-workers* if their 
argumentations were conclusive in every respect. Only 
their most weighty argument of their most recent publica- 
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tion may be considered. Either the decay reaction occurs 
independently of the surface or their wall reaction is 
erroneous. Both cannot be assumed at the same time. 
Conceding for a moment that the wall reaction be 
negligible, a homogeneous gas reaction of the first order is 
improbable because it will be difficult to suggest a mech- 
anism which would fit such an assumption. A mono- 
molecular decay should therefore be explained as a wall 
reaction. It seems possible’ to oppress this wall reaction 
even more than Steiner and Wicke*® were able to do. One 
may argue that active hydrogen in as small a concentration 
as 0.001 to 0.025 percent® would have exactly the same 
properties if explained on the basis of the well-known 
atomic modification. This argumentation may be followed 
up in every instance. 

Wrede’s orifice method'® and Keussler’s study of the 
absorption spectrum of H-atoms!" would be sufficient 
proof to substantiate the atomic nature of that hydrogen 
which is obtained by Wood’s method, unless the bulk of 
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other investigations which were based on this assumption 
is taken as proof in itself. A decision in favor of H* appears 
exceedingly difficult by the use of merely chemical methods. 
GEORG R. SCHULTZE 
Universitat Berlin (Germany), 
Physikalisch-chemisches Institut, 
April 15, 1935. 
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